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PREFACE 


Since  July  1,  1975,  Lawrence  Livermore  Laboratory  (LLL)  has  been  participating  in  the  High 
Altitude  Pollution  Program  sponsored  by  the  U.S.  Department  of  Transportation’s  Federal  Aviation  Ad¬ 
ministration.  This  report  describes  the  major  accomplishments  and  significant  findings  during  the  fiscal  year 
ending  September  30,  1979,  for  work  performed  at  LLL  under  Reimbursable  Agreement  DOT-FA79WAI- 
034.  Two  major  research  areas  are  covered  by  this  agreement:  (1)  numerical  modeling  of  the  atmospheric 
response  to  stratospheric  perturbations,  and  (2)  the  processing,  archiving,  and  analysis  of  satellite  ozone  data. 
Each  of  these  research  areas  has  been  divided  into  a  number  of  subtasks,  and  the  successful  accomplishment 
of  these  subtasks  has  required  contributions  and  cooperation  from  many  participants.  The  work  reported  here 
should  be  considered  the  collective  effort  of  all  those  listed  below. 

Scientific  Administration 

Joseph  B.  Knox,  Division  Leader 
Frederick  M.  Luther,  Principal  Investigator 
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Annual  report  of 
Lawrence  Livermore  Laboratory 
to  the  FAA  on  the  high  altitude 
pollution  program — 1979 


1.  INTRODUCTION 


This  annual  report  documents  the  progress 
made  on  research  in  support  of  the  High  Altitude 
Pollution  Program  (HAPP)  during  FY-1979.  Work 
performed  during  previous  years  at  Lawrence  Liver¬ 
more  Laboratory  (LLL)  is  described  in  our  earlier 
annual  reports  (Luther  el  al .,  1976,  1977,  1978),  A 
review'  of  progress  made  during  previous  years  of  ef¬ 
fort  on  HAPP  is  also  included  in  this  year’s  report. 

During  the  past  decade,  a  great  deal  of  effort 
has  been  devoted  to  understanding  the  potential  en¬ 
vironmental  impact  of  high-altitude  aircraft  flights. 
In  1970,  the  Climatic  Impact  Assessment  Program 
(ClAP)  was  created  to  evaluate  the  impact  of  the 
future  operation  of  supersonic  transport  aircraft 
(SST’s)  in  the  stratosphere.  This  research  program 
was  sponsored  by  the  U  S.  Department  of 
Transportation  and  was  completed  in  1975  with  the 
publication  of  a  series  of  monographs  that  assessed 
the  potential  chemical,  climatic,  biological,  and 
economic  effects  of  stratospheric  aircraft  emissions. 
Additional  assessments  were  prepared  by  the 
Climatic  Impact  Committee  of  the  National 
Academy  of  Sciences  (National  Research  Council, 
1975),  by  the  British  (COMESA,  1976).  and  by  the 
French  (COVOS.  1976). 

The  High  Altitude  Pollution  Program  was  ini¬ 
tiated  in  1975  by  the  Federal  Aviation  Administra¬ 
tion  to  extend  the  investigations  carried  out  during 
CIAP  so  as  to  ensure  that  stratospheric  aircraft 
emissions  would  not  result  in  unacceptable  environ¬ 
mental  effects.  The  environmental  issues  that  were 
raised  relative  to  SST's  pertain  to  all  types  of  air¬ 
craft  operating  in  the  stratosphere.  Consequently, 
the  flight  regimes  of  interest  to  HAPP  include  sub¬ 
sonic  .  supersonic,  and  hypersonic  aircraft.  The  air¬ 
craft  engine  effluents  of  primary  concern  in  terms  of 
impact  on  o/one  are  oxides  of  nitrogen  (NO  and 
NOs)  and  ssater  vapor.  Perturbations  other  than 
aircraft  emissions  have  also  been  proposed  as  hav¬ 
ing  a  potentially  significant  impact  on  o/one  (e.g.. 


Cl  and  CIO  from  chlorofluoromethanes).  The 
HAPP  research  effort  also  considers  these  other 
perturbations  in  order  to  address  the  overall 
problem  of  the  combined  effects  of  stratospheric 
I  urbations. 

Lawrence  Livermore  Laboratory  has  par¬ 
ticipated  in  HAPP  since  July  1975  under  in¬ 
teragency  agreements  DOT-TSC-76-1,  DOT- 
FA76WA 1-653,  and  DOT-FA79WAI-034.  Primary 
research  emphasis  at  LLL  is  on  numerical  modeling 
of  the  atmospheric  response  to  stratospheric  pertur¬ 
bations.  The  modeling  effort  at  LLL  covers  four 
major  research  areas:  photochemical  kinetics, 
coupled  kinetics  and  transport,  radiative  transfer, 
and  meteorological  analysis. 

Photochemical  kinetics  modeling  spans  gas- 
phase  reaction  processes  and  heterogeneous  (i.e., 
gas-particle)  reaction  processes,  constant  and  diur- 
nally  (and  seasonally)  varying  photochemical  reac¬ 
tion  rates,  a  wide  range  of  ambient  temperatures 
and  pressures,  neutral  and  ionized  species,  processes 
with  both  accurately  and  poorly  known  rates  of 
reaction,  and  species  of  both  acknowledged  and 
questionable  importance  in  the  stratosphere.  Con¬ 
siderable  attention  has  been  directed  toward  evalua¬ 
tion  of  the  sensitivity  of  reaction  mechanisms  to 
deficiencies  in  our  knowledge  of  reaction  rates, 
quantum  yields,  and  reaction  mechanisms. 

The  transport-kinetics  models  are  intended  for 
time-dependent  calculations  of  the  response  of  the 
atmosphere  to  stratospheric  perturbations.  These 
models  incorporate  interactive  radiative  and 
chemical  processes  coupled  with  transport,  which 
must  be  parameterized  in  one-  and  two-dimensional 
models.  These  models  are  designed  to  facilitate  per¬ 
turbation  and  sensitivity  studies. 

Radiative  transfer  modeling  includes:  applica¬ 
tion  of  highly  detailed  models  to  assessment  of  the 
effect  on  the  radiation  budget  of  perturbations  to 
the  stratospheric  composition,  and  development  of 
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simplified  radiative  transfer  models  for  inclusion  in 
the  transport-kinetics  models.  Considerable  atten¬ 
tion  has  been  directed  toward  including  interactive 
radiative  processes  in  the  transport-kinetics  models 
in  order  to  assess  possible  feedback  mechanisms. 

Meteorological  analysis  provides  guidance  and 
support  for  the  modeling  effort  through  in-depth 
studies  of  atmospheric  processes  and  phenomena.  A 
primary  responsibility  of  this  research  is  to  develop 
test  situations  to  validate  various  aspects  of  the 
numerical  models. 

In  August  1976,  LLL’s  participation  in  HAPP 
was  extended  to  include  a  feasibility  study  to  deter¬ 
mine  whether  good  quality  total  ozone  data  could 
be  derived  from  infrared  radiance  measurements  by 
a  multi-channel  filter  radiometer  (MFR)  sensor 
carried  aboard  a  series  of  satellites  operated  by  the 
U.S.  Air  Force.  These  satellites  are  a  part  of  the 
Defense  Meteorological  Satellite  Program  Block  5D 
series.  The  first  four  satellites  in  this  series  carry 
MFR  sensors  from  which  total  ozone  data  may  be 
derived.  The  MFR  sensor  is  the  first  cross-track 
scanning  sensor  capable  of  measuring  ozone,  and  it 
has  a  higher  spatial  resolution  and  lower  noise  level 
than  previous  instruments.  The  first  satellite  began 
transmitting  MFR  measurements  in  March  1977, 
and  MFR  data  will  be  received  from  these  satellites 
into  the  early  I980’s.  These  MFR  data  should  be 
useful  for  determining  the  temporal  and  spatial 
variability  of  ozone  over  this  period  of  time. 

The  feasibility  study  was  completed  in  May 
1978,  with  the  successful  processing  of  20  days  of 
data  taken  during  1977.  The  quality  of  the  derived 
total  ozone  data  was  demonstrated  by  comparison 
with  corresponding  ozone  data  obtained  at  selected 
stations  in  the  world  surface  network  of  Dobson 
spectrophotometer  observatories.  A  description  of 
the  methodology  used  and  global  maps  of  the  total 
ozone  data  for  the  20  days  are  contained  in  the 
report  by  Lovill  el  al.  (1978). 

Following  completion  of  the  feasibility  study, 
LLL's  effort  has  been  directed  toward  archiving  and 
processing  the  MFR  data.  Before  routine  data 
processing  could  begin,  an  automated  data-basing 
system  had  to  be  developed.  Each  MFR  sensor  is 


capable  of  taking  up  to  67,500  observations  (each 
observation  consisting  of  16  radiance  measure¬ 
ments)  per  day,  and  as  many  as  three  satellites  were 
taking  data  concurrently,  so  several  million  pieces 
of  data  are  included  in  the  data  base  for  each  day  of 
coverage.  Consequently,  the  archiving  and  process¬ 
ing  of  the  data  would  be  difficult  without  an 
automated  data  base.  We  have  also  been  evaluating 
various  aspects  of  the  ozone  retrieval  methodology 
and  making  refinements.  Thus,  the  methodology 
used  for  regular  data  processing  wul  be  improved 
over  that  used  in  the  feasibility  study. 

The  satellite  data  processing  task  has  three  ma¬ 
jor  areas  of  effort:  (1)  radiative  transfer  and  ozone 
retrieval  methodology,  (2)  data  processing,  archiv¬ 
ing,  and  distribution,  and  (3)  data  analysis,  inter¬ 
pretation,  and  quality  assurance.  Work  on  these 
tasks  is  supported  in  part  by  the  National 
Aeronautics  and  Space  Administration  (NASA). 

The  first  area  of  effort  involves  detailed 
transmittance  and  simulated  radiance  calculations, 
comparisons  between  models  and  laboratory 
measurements,  and  ozone  retrieval  model  develop¬ 
ment  using  regression  analysis.  Sensitivity  studies, 
error  analysis,  and  special  studies  with  satellite 
radiances  are  also  included. 

The  second  area  of  effort  involves  MFR  data 
calibration,  optimization  of  data  processing 
methodology,  design  and  implementation  of  an 
automated  data  base,  data  processing,  mapping  of 
ozone  data,  and  data  distribution.  In  addition, 
satellite  overflight  times  are  predicted  and  dis¬ 
tributed  to  over  36  Dobson  observatories  each 
month  so  that  special  Dobson  spectrophotometer 
measurements  can  be  taken  for  comparison  with  the 
MFR  ozone  data. 

The  third  area  of  effort  involves  development 
of  the  Dobson  and  ozonesonde  data  bases,  develop¬ 
ment  of  quality-assurance  techniques  and  statistical 
analyses,  comparison  of  the  MFR  ozone  data  with 
other  data  sources,  and  analyses  of  the  data  on 
various  temporal  and  spatial  scales.  The  data  will  be 
analyzed  for  both  anthropogenic  and  non- 
anthropogenie  effects. 
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2.  ATMOSPHERIC  MODELING 


A  fundamental  tool  in  the  LLL  effort  in  at¬ 
mospheric  modeling  has  been  the  one-dimensional 
transport-kinetics  model.  This  model  includes  as 
complete  a  set  of  the  important  chemical  and  pho¬ 
tochemical  reactions  as  is  feasible  and  is  designed 
for  time-dependent  perturbation  and  sensitivity 
studies.  This  model  currently  includes  39  chemical 
species  and  134  chemical  and  photochemical  reac¬ 
tions.  Species  concentrations  are  computed  at  44 
levels  in  the  atmosphere,  extending  from  the  ground 
to  55  km.  Vertical  transport  is  parameterized  using 
a  one-dimensional  diffusion  formulation  that 
describes  hemispheric-average  net  vertical 
transport.  The  model  can  include  temperature 
coupling  between  changes  in  composition  and  reac¬ 
tion  rate  coefficients  and  hydrostatic  adjustment  of 
pressure  and  density  as  user  options,  but  these  op¬ 
tions  are  not  utilized  in  our  standard  model  calcula¬ 
tions.  The  one-dimensional  transport-kinetics 
model  is  described  in  more  detail  in  Appendix  A. 

In  addition  to  the  one-dimensional  transport- 
kinetics  model,  we  have  used  radiative  transfer 
codes  of  varying  complexity  for  several  studies.  We 
are  also  in  the  process  of  developing  a  two- 
dimensional  transport-kinetics  model. 

Since  the  beginning  of  the  HAPP  program, 
LLL  has  completed  the  following  major  research 
studies  and  activities: 

•  Sensitivity  of  NO*  catalytic  ozone  destruc¬ 
tion  to  uncertainties  in  rate  coefficients, 

•  Sensitivity  of  ozone  reduction  from 
chlorofluoromethanes  to  parameter  uncertainties. 

•  Effects  of  changes  in  O,  and  NO?  on  at¬ 
mospheric  solar  absorption. 

•  Effect  of  NO  photolysis  on  NOv  mixing 
ratios. 

•  Analysis  of  chemical  rate  data. 

•  Effect  of  multiple  scattering  on  species  con¬ 
centrations  and  model  sensitivity. 

•  Effect  of  changes  in  stratospheric  water 
vapor  on  ozone  reduction  estimates. 

•  Effect  of  atmospheric  nuclear  testing  dur¬ 
ing  the  1 950's  and  1960's  on  ozone. 

•  Effect  of  agriculture  on  stratospheric 
ozone. 

•  Calculation  of  global  tropospheric  OH 
distributions. 

•  Tropospheric  lifetimes  and  potential  effects 
on  ozone  of  halocarbons. 


•  Analysis  of  global  budgets  of  halocarbons. 

•  Effects  of  changes  in  total  o/one  and 
receiver  orientation  on  received  erythema  dose. 

•  Effect  of  possible  variations  in  solar  ul¬ 
traviolet  flux  on  stratospheric  ozone. 

•  Effects  of  a  solar  eclipse  on  stratospheric 
chemistry. 

•  Potential  effects  of  space  shuttle  emissions. 

•  Potential  effects  of  solar  power  satellite 
launch  vehicles. 

•  Effects  of  a  massive  pulse  injection  of  NO* 
into  the  stratosphere. 

•  Effects  of  stratospheric  perturbations  on 
the  earth’s  radiation  budget  and  the  climatic  im¬ 
plications. 

In  addition  to  the  listed  special  studies,  we  have 
continued  to  work  on  model  improvements  and 
comparisons  between  model  calculations  and  obser¬ 
vations.  These  comparisons  are  important  to  the 
verification  of  various  aspects  of  the  model.  We 
have  compared  computed  species  profiles  with  ob¬ 
servations  for  the  ambient  atmosphere  as  well  as 
simulations  of  natural  and  man-made  perturba¬ 
tions.  Several  basic  types  of  data  have  proven  to  be 
useful:  height  and  latitudinal  distributions  of  in¬ 
dividual  species  and  groups  of  species,  partitioning 
of  related  species,  and  local  diurnal  and  seasonal 
variations  of  individual  species.  Interpretation  of 
the  results  from  such  comparisons  is  not  always 
simple  because  of  the  significant  level  of  local 
variability  of  species  concentrations.  It  is  the  overall 
reasonableness  of  the  models  and  their  predictions 
that  provide  a  high  level  of  confidence  in  their  basic 
correctness.  Of  course,  this  is  often  a  matter  of  judg¬ 
ment. 

The  one-dimensional  transport-kinetics  model 
has  been  used  to  assess  the  potential  chemical  ef¬ 
fects  of  several  man-made  perturbations  to  the 
stratosphere.  Whenever  significant  changes  occur  to 
the  model  chemistry  or  to  the  treatment  of  physical 
processes,  the  assessment  studies  are  repeated  in  or¬ 
der  to  assess  the  effect  of  the  changes  on  the  model 
sensitivity. 

The  research  activities  that  were  completed 
during  the  past  year  are  described  in  the  following 
sections.  One  of  the  major  accomplishments  during 
the  past  year  was  the  preparation  of  a  detailed 
assessment  report  on  the  potential  effects  of  aircraft 
emissions  (Luther  et  ol ..  1979).  The  reader  should 
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refer  to  that  report  for  additional  material  and  a 
more  detailed  discussion  of  the  topics  presented 
here. 

2.1  MODEL  SIMULATIONS  OF 
STRATOSPHERIC  OBSERV  ATIONS 

Stratospheric  Composition 

Comparison  with  observational  data  on  trace 
species  concentrations  is  an  important  aspect  of 
validating  the  performance  of  numerical  models  of 
the  stratosphere.  Although  comparisons  with  obser¬ 
vations  are  a  necessary  part  of  model  validation, 
these  comparisons  alone  are  not  sufficient  to 
validate  the  performance  of  the  model  because  we 
know  from  past  experience  (Duewer  el  al.,  1977) 
that  models  with  significantly  different  chemistry 
and  sensitivity  can  predict  very  similar  ambient 
species  profiles. 

Two  somew  hat  different  sets  of  chemical  reac¬ 
tion  rate  coefficients  and  photolysis  cross  sections 
were  used  in  the  calculations  presented  in  this 
report.  Most  of  the  calculations  were  performed 
prior  to  the  NASA  Workshop  at  Harpers  Ferry, 
West  Virginia  (June  1979),  and  utilized  chemical 
rates  based  primarily  on  JPL  (1979).  After  the 
NASA  Workshop,  we  repeated  several  of  the  per¬ 
turbation  calculations  using  chemical  rates  based  on 
the  recommendations  of  the  chemistry  panel  at  the 
NASA  Workshop.  These  two  chemistries  are 
described  in  Appendix  A  as  1979a  and  1979b 
chemistries. 

In  this  section,  we  compare  the  results  of  our 
current  one-dimensional  model  using  1979b 
chemistry  with  measured  trace  species  concentra¬ 
tions  Much  of  the  measurement  data  used  in  the 
comparisons  was  derived  from  LLL-supplied  infor¬ 
mation  communicated  at  the  NASA  Harpers  Ferry 
Workshop.  As  such,  some  is  preliminary  in  nature. 
The  evaluation  of  measurements  is  strongly  influ¬ 
enced  by  that  of  the  NASA  panels  as  reflected  at  the 
Workshop.  Our  discussion  in  this  section  closely 
parallels  some  of  that  in  the  NASA  document, 
reflecting,  in  part,  our  participation  in  writing  parts 
of  that  document 

<)(  JP>.  There  is  good  agreement  between 
model-calculated  distributions  and  the  few  available 
measurements.  However,  the  OUP)  concentration  is 
in  a  near  photoequilibrium  with  Oi,  and  the  ozone 
column  used  in  the  calculations  is  calculated  for 


equinox  conditions  w'hereas  the  measurements  are 
winter  measurements.  Because  Oi  was  not 
measured  at  the  same  time  that  O  was,  there  is  some 
danger  that  the  apparent  good  agreement  is  for¬ 
tuitous. 

Oj.  Ozone  profiles  vary  significantly  with 
latitude  and  season.  Total  column  ozone  ranges 
from  about  200  m  •  atm  •  cm  in  the  tropics  to 
about  400  m  •  atm  •  cpi  near  the  poles.  The  altitude 
of  the  peak  concentration  is  lower  near  the  poles 
than  at  the  equator.  The  model-derived  O3  profile 
resembles  observations  of  O3  taken  at  middle 
latitudes. 

OH,  H02.  Within  the  plausible  error  in  the 
measurements,  there  appears  to  be  reasonable 
agreement  between  theoretical  calculations  and  the 
few  available  measurements  of  OH  and  H02  in  the 
stratosphere.  Data  are  available  only  between  28 
and  28  km  for  H02.  The  point-to-point  variability 
of  the  measurements  for  both  species  exceeds 
measurement  precision.  This  suggests  that  the  dif¬ 
ferent  air  masses  measured  may  have  significantly 
different  histories.  It  would  be  desirable  if,  in  the 
future,  concurrent  measurements  of  various  species 
(OH,  HOs,  O3.  H20,  etc.)  would  be  made.  If  such 
concurrent  measurements  become  available,  it  may¬ 
be  possible  to  explain  local  variability  ,  as  well  as  to 
usefully  examine  ratios  (e.g.,  HO;/HO)  in  relation 
to  those  expected  theoretically  .  Comparison  of  the 
total  OH  column  calculated  by  current  one- 
dimensional  models  with  the  observations  of  Bur¬ 
nett  and  Burnett  (1979)  does  not  indicate  significant 
differences  in  average  column  content. 

NO.  When  measurements  taken  at  middle 
latitudes  and  near  noon  are  compared  to  calculated 
noon  profiles,  the  calculated  profiles  are  signifi¬ 
cantly  larger  than  the  observations  above  10  km. 
The  maximum  difference  occurs  in  the  upper 
stratosphere  where  NO  is  the  dominant  NOx 
species.  One  contribution  to  the  difference  between 
models  and  observations  is  the  overprediction  of  the 
total  NOx  content  by  current  models.  The  apparent 
NO  maximum  near  40  km  lends  some  credibility  to 
the  hy  pothesis  that  the  treatment  of  NO  photolysis 
may  underestimate  the  actual  upper  stratospheric 
NO  sink. 

Recent  measurements  of  NO  (discussed  by  H. 
Schiff  and  others  at  the  N.A.T.O.  Advanced  Study 
Institute  on  Atmospheric  Ozone  held  in  Portugal  in 
October  1979)  find  higher  concentrations  in  the  up¬ 
per  stratosphere  than  previous  measurements. 
These  measurements  agree  better  with  model 
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calculations.  However,  the  reason  for  disagreement 
with  previous  observations  (whether  it  is  due  to 
natural  variability  or  measurement  errors)  must  be 
resolved. 

HNO,.  Current  models  predict  more  HNO^ 
than  is  observed  above  about  25  km.  In  part  this 
may  simply  reflect  the  higher-than-observed  NOv 
concentrations  currently  calculated.  However,  this 
does  not  appear  to  be  sufficient  to  explain  the  entire 
difference,  nor  is  it  likely  that  model  overprediction 
of  OH  between  25  and  35  km  is  responsible  for  the 
difference,  since  model  OH  and  observations  are  in 
good  agreement  between  28-35  km. 

Total  NOv  Our  theoretical  model  presently 
calculates  approximately  a  factor  of  2  more  NOn 
(NO  +  NO;  +  HNO.,)  than  has  been  measured. 
However,  it  should  be  noted  that  in  situ  measure¬ 
ments  of  all  three  species  are  not  available  with 
readily  comparable  techniques  to  permit  an  unam¬ 
biguous  comparison  between  theory  and  measure¬ 
ments. 

Part  of  the  reason  for  the  difference  between 
model  and  observation  may  be  that  the  mode! 
predicts  too  much  N:0  in  the  upper  stratosphere, 
since  reaction  of  N:0  with  0( ' D)  is  the  primary 
source  for  NOv  in  the  stratosphere.  More  measure¬ 
ments  of  NsO  in  the  upper  stratosphere  would  be 
necessary  to  confirm  this  hypothesis.  The  N:0  con¬ 
tent  in  one-dimensional  models  is  sensitive  to  the 
choice  of  vertical  transport  coefficient  used  in  the 
calculation. 

Other  possible  factors  that  might  cause  the 
models  to  overpredicl  NOv  would  include:  (I)  A 
lower  rate  coefficient  for  0(!D)  +  NsO  ->  2NO 
than  is  currently  recommended.  The  data  for  this 
process  are  not  in  good  agreement.  However,  if 
anything,  they  suggest  a  faster  rate  coefficient.  (2) 
O('D)  may  be  overestimated.  This  could  happen  if 
the  quantum  yield  for  OCD)  production  from  O* 
photolysis  has  a  stronger  temperature  dependence 
than  is  used  in  the  models.  (3)  The  rate  at  which  NO 
photolysis  occurs  may  be  underestimated  in  the 
models.  Frederick  and  Hudson  (1 979) estimate  their 
results  to  be  uncertain  by  a  factor  of  3.  Adoption  of 
their  treatment  decreased  the  photolysis  rate  by  a 
factor  of  about  3  from  the  previous  value  in  the 
model  which  resulted  in  an  increase  of  the  model- 
calculated  NOx  by  about  40'i. 

HNO,/N02.  The  ratio  of  UNO,  NO:  cur¬ 
rently  calculated  by  models  is  much  higher  than  is 
observed  hy  I  vans  el  a/.  (1976)  but  is  within  the 
error  bounds  ot  other  measurements.  In  none  of 


these  measurements  were  both  species  measured  by 
similar  techniques.  Concurrent  measurements  taken 
at  a  given  point  with  the  same  technique  are  needed 
to  establish  whether  a  discrepancy  actually  exists. 
Even  then,  because  of  the  difference  in  chemical 
lifetimes,  it  would  be  difficult  to  determine  the  time 
history  of  both  NO,  and  HNOi  that  led  to  the  ratio 
observed.  If  a  discrepancy  does  exist,  one  possible 
explanation  would  be  that  current  estimates  of  OH 
in  the  lower  stratosphere  (below  approximately 
30  km)  are  too  large.  Another  possibility  would  be 
that  the  sinks  for  HNO}  are  underestimated. 

N2O5.  We  currently  calculate  a  NyOs  concen¬ 
tration  of  approximately  2.5  ppb  near  sunrise.  This 
can  be  compared  to  the  value  of  2  ppb  deduced  for 
NSOs  by  Evans  el  at.  (1978)  based  upon  measure¬ 
ments  of  other  species  a  few  hours  after  sunrise.  It 
should  be  noted  that  the  measurement  was  not 
made  at  the  latitude  and  season  the  model  most 
nearly  approximates.  Until  further  measurements 
are  made,  the  comparison  is  inconclusive. 

Cl,  CIO.  The  single  observed  Cl  profile  (and 
one  additional  Cl  observation  at  a  single  altitude) 
and  calculated  Cl  profiles  agree  within  the  measure¬ 
ment  uncertainty  of  ±359f.  However,  the  measure¬ 
ments  suggest  a  larger  vertical  gradient  in  Cl  mixing 
ratio  than  is  predicted.  More  Cl  measurements  are 
needed  to  determine  whether  any  seasonal  or  spatial 
variations  exist  and.  if  so,  whether  these  are  consis¬ 
tent  with  model  predictions. 

For  CIO.  at  the  mixing  ratio  maximum,  with 
the  exception  of  the  14  July  1977  measurement  by 
Anderson  et  at.  (1979),  the  resonance  fluorescence 
measurements  are  w  ithin  a  factor  of  about  2  of  the 
predicted  CIO  abundances.  The  September  through 
December  CIO  measurements  show  a  significantly 
sharper  than  predicted  gradient  below  the  mixing 
ratio  peak.  Except  for  the  summer  ot  1979,  summer 
measurements  give  more  CIO  than  the  September 
through  December  measurements.  The  14  July  1977 
measurement  gives  six  limes  more  CIO  than  is 
predicted  near  the  mixing  ratio  peak,  and  exceeds 
by  about  a  factor  of  3  the  total  chlorine  mixing  ratio 
believed  to  be  present  in  the  stratosphere.  Predicted 
variations  of  the  Cl  and  CIO  profiles  w  ith  season  are 
very  small  compared  to  the  range  of  observed 
profiles.  A  broader  data  base  is  needed  to  explain 
'he  difference  111  these  results  as  well  as  the  large 
(factor  of  10)  variation  in  the  CIO  profiles  observed 
but  not  predicted  by  one-dimensional  models 

The  differences  between  theoretical  calcula¬ 
tions  and  observations  of  the  CIO  gradient  below 
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34  km  also  require  attention.  At  present,  specula¬ 
tions  for  explaining  these  differences  range  from  the 
possibility  of  a  missing  chemical  reaction  or  an  in¬ 
correct  rate  constant  in  present  theory  to  the 
possibility  of  it  being  an  artifact  of  the  transport 
parameterization  used  in  current  models,  or  an  ar¬ 
tifact  of  the  measurement  techniques. 

Cl/CIO  Ratio.  The  Cl/CIO  ratios  from  the  28 
July  1976  resonance  fluorescence  measurements 
(Anderson  el  al.,  1979)  agree  to  within  measurement 
uncertainty  with  most  calculations.  It  is  difficult  at 
present  to  make  any  conclusions  regarding  this 
comparison  because  of  the  experimental  uncertain¬ 
ties,  the  limited  amount  of  data,  the  apparent 
variability  of  CIO,  and  the  poor  agreement  between 
models  and  observations  for  CIO. 

HCI.  Although  there  may  be  significant  dis¬ 
agreement  (as  large  as  30-40%  when  compared  with 
the  1R  spectra  data,  much  larger  if  compared  to  the 
filter  or  radiometer  data)  above  30  km,  the  uncer¬ 
tainty  in  the  observational  data  in  this  altitude 
range  is  such  as  to  probably  encompass  the  model 
predictions. 

CIONOj.  The  only  observational  data  for 
CIONOt  are  from  one  flight  by  Murcray  (1979). 
There  is  excellent  agreement  if  the  fast  production 
rate  for  CIONCU  is  used  in  the  calculations,  whereas 
there  is  a  factor  of  3  difference  if  the  slow  er  produc¬ 
tion  rate  is  used  (see  Section  2.4  discussion 
regarding  the  rate  of  production  of  CI0N02). 

Until  further  data  are  available,  it  is  difficult  to 
make  any  conclusive  statements  regarding  this 
species.  However,  the  limited  amount  of  data 
available  does  not  indicate  any  gross  discrepancies 
between  theory  and  observations. 

Total  C'lX.  There  now  exist  a  limited  number  of 
measurements  for  Cl,  CIO,  HCI,  and  CIONCU  in  the 
stratosphere.  Except  for  three  measurements  of  CIO 
taken  in  summer  (Anderson,  1979:  Menzies,  1979), 
there  does  not  appear  to  be  any  large  discrepancy 
between  the  total  amount  of  CIX  predicted  and  that 
observed.  However,  the  summer  measurements  of 
CIO  are  so  large  as  to  require  much  more  total  CIX 
than  known  sources  are  capable  of  producing.  This 
implies  the  possibility  of  an  unknown  source  for 
CIX  in  the  atmosphere,  an  error  in  the  observation, 
or  some  unusual  episodic  event 

Model  Simulations  of  Past  Perturbations 

Possible  Long-Term  Solar  Variations  and  Ef¬ 
fects  on  Temperature  and  Ozone.  In  the  stud}  of 


potential  anthropogenic  influences  on  the 
stratospheric  ozone  budget  it  is  necessary  to  first 
understand  the  natural  variations  of  ozone.  In  order 
to  achieve  early  detection  of  anthropogenic  trends 
and  to  establish  the  magnitude  of  such  effects,  we 
must  differentiate  the  natural  stratospheric  ozone 
variations  from  the  man-made  effects  that  are  of 
comparable  time  scales.  The  suggested  correlation 
between  ozone  and  cunspot  number  is  of  primary 
interest.  If  this  relation  is  real,  then  for  the  next  few 
years  the  theoretically  predicted  effect  of  CFM's  on 
ozone  may  not  be  directly  detectable  unless  the  ef¬ 
fect  of  this  ozone-solar  cycle  relation  can  be  quan¬ 
tified  and  removed  from  the  ozone  data.  At  this 
time,  the  most  plausible  mechanism  coupling  the 
sunspot  cycle  and  atmospheric  ozone  is  the  change 
in  solar  flux  between  180  and  340  nnt  from  solar 
maximum  to  solar  minimum  (Heath  and 
Thekaekara,  1977:  Callis  and  Nealy.  1978:  Penner 
and  Chang,  1978).  Theoretically,  a  variation  in  solar 
flux  between  180  and  340  nm  of  approximately  30% 
(maximum  to  minimum)  can  lead  to  local  ozone 
changes  as  large  as  10%  near  35  km,  and  total  ozone 
column  changes  of  approximately  5%  (Penner  and 
Chang,  1978).  This  is  comparable  to  the  predicted 
CFM  effect  of  ozone  decreases  of  up  to  5%  around 
the  year  1990.  Consequently,  the  possible  influence 
of  periodic  or  aperiodic  changes  in  solar  UV  flux  in¬ 
tensities  tends  to  obscure  changes  due  to  the  present 
levels  of  CFM’s,  and  this  may  continue  for  the  next 
decade.  The  same  is  also  true  for  the  local  ozone 
concentrations  in  the  upper  stratosphere. 

The  data  that  support  a  variation  in  solar  UV 
fluxes  during  a  solar  cycle  are  limited  (Heath  and 
Thekaekara,  1977)  and  need  to  be  confirmed  with 
other  independent  measurements.  Simon  (1 978)  has 
pointed  out  the  difficulty  in  the  absolute  calibration 
of  the  measuring  instruments,  especially  those 
aboard  a  satellite. 

Since  the  direct  monitoring  of  solar  UV  (luxes 
and  ozone  are  both  of  limited  use  in  establishing  the 
ozone-solar  cycle  relation,  we  are  left  with  the 
monitoring  of  trace  species  other  than  ozone  as  a 
possible  independent  method  for  validating  the  ef¬ 
fects  of  solar  LIV  flux  variations.  We  examined  the 
possible  variations  of  31  trace  species  that  result 
from  hypothetical  solar  U V  \ an.it ions  with  the  LLL 
time-dependent  model  of  stratospheric  chemistry, 
including  sell-consistently  calcu1  ted  stratospheric 
temperature  and  atmospheric  adiustment  to 
hydrostatic  equilibrium  (Penner  and  (  hang.  1 97V). 
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From  this  study  we  identified  NiO  us  a  most  likelv 
candidate  for  monitoring  to  study  coupling  between 
variations  in  solar  UV  and  atmospheric  composi¬ 
tion.  Furthermore,  there  remains  the  question  of  the 
effect  of  transport-related  variability  that  can  only 
he  resolved  with  a  set  of  upper  stratospheric  data. 

The  Atmospheric  Nuclear  Tests  of  the  1950’s 
and  I960’s.  Past  atmospheric  tests  of  nuclear  devices 
larger  than  about  one  megaton  (TNT)  yield  provide 
a  potentially  significant  source  of  NOx  to  the 
stratosphere.  During  the  late  1950's  and  early 
I960's,  large  numbers  of  such  tests  were  carried  out, 
and  the  NOx  released  to  the  stratosphere  in  1961 
and  1962  should  have  been  comparable  to  the 
amount  that  would  be  released  by  about  2.7  years 
operation  of  a  large  SST  fleet  (i.e.,  one  emitting  ap¬ 


proximately  10x  molecules  cm-V1).  The  paper  by 
Chang,  Duewer  and  Wuebbles  (1979)  gives  an  ac¬ 
count  of  the  predicted  effects  through  early  1978. 
Figure  1  gives  the  computed  o/one  column  change 
as  a  function  of  time.  For  computations  made  with 
1979b  chemistry  (see  Appendix  A),  the  two  curves 
give  the  ozone  depletion  calculated  for  injection 
altitudes  based  on  a  parameterization  by  Foley  and 
Ruderman  (1972)  or  on  the  lower  altitude  observa¬ 
tions  of  Seitz  a  al.  (1968).  The  two  treatments  are 
thought  to  bound  the  plausible  range  of  injection 
altitudes. 

Both  treatments  predict  maximum  ozone 
depletions  of  less  than  2%  for  the  largest  annual- 
average  change.  The  largest  change  is  only  about 
1.25%  il  a  1-2-1  smoothing  function  is  applied  to 
annual  perturbations.  The  computed  changes  are 


HG.  I  Computed  change  in  total  ozone  resulting  from  atmospheric  nuclear  testing  during  1950’s  and  I960’s. 
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clearly  within  the  observed  variability  in  the  ozone  ion-pair  are  used.  However,  the  agreement  between 

record  as  analyzed  by  Angell  and  Korshover  (1978).  observation  and  computation  at  lower  latitudes  was 

Figure  2  gives  the  computed  changes  at  se-  less  satisfactory,  although  still  qualitatively  en- 

lected  altitudes.  The  ozone  concentration  decreases  couraging.  Borucki  et  at.  (1978)  report  similar 

at  all  altitudes  above  25  km.  However,  at  lower  findings,  and  also  find  that  their  model  substantially 

altitudes  the  concentration  of  ozone  was  computed  underpredicts  the  ozone  perturbations  at  altitudes 

to  increase  significantly,  which  accounts  for  the  net  near  30  km. 

increase  in  total  ozone  prior  to  1962  shown  in  f  he  August  1972  PCA  events  seem  to  provide  a 

Fig.  I.  Some  of  these  computed  changes  are  much  useful  test  of  the  short-term  (several  days)  response 

larger  than  the  change  in  the  ozone  column.  of  stratospheric  ozone  to  NO,  increases  above 

However,  the  predicted  ozone  variations  at  specific  about  40  km.  Unfortunately,  the  test  is  not  directly 

altitudes  are  still  only  about  the  same  magnitude  as  applicable  to  lower  altitude,  midlatitude  NO,  per- 

the  variations  in  ozone  concentrations  determined  turbations  occurring  over  long  (several-year)  time 

bv  the  IJmkehr  technique.  periods.  The  nature  of  the  NO,  injection  and  of  the 

Thus,  the  model-predicted  response  to  the  at-  O,  data  do  not  permit  resolution  of  questions  about 

mospheric  nuclear  tests  does  not  lead  to  conflict  ihe  adequacy  of  model  simulations  of  transport 

with  the  ozone  record,  but  it  also  does  not  seem  to  phenomena  or  of  the  chemistry  of  the  lower 

explain  much  of  the  observed  ozone  variation.  The  stratosphere.  Thus.  PCA  events  provide  the  only 

predicted  effect  on  ozone  of  the  nuclear  tests  should  phenomenon  for  w  hich  models  forecast  an  observ- 

be  considered  in  any  attempt  to  model  ozone  varia-  a^e  and  observed  stratospheric  ozone  perturbation 

tions  during  the  1960’s  since  the  predicted  changes  'n  response  to  an  NO,  change,  but  the  nature  of  the 

in  ozone  at  specific  altitudes  are  comparable  to  both  perturbation  and  the  response  differ  from  the 

the  observed  variations  and  to  the  variations  problems  associated  with  SST  operations  to  an  ex- 

calculated  for  such  phenomena  as  the  hypothesized  lenl  lhat  precludes  its  direct  use  as  a  calibration 

variation  of  ultraviolet  light  tied  to  the  solar  cycle,  Poinl  for  SST  predictions  of  column  ozone  changes, 

or  the  predicted  change  in  ozone  from  CFM  Nonetheless,  the  good  agreement  above  40  km  is 

production  through  1978.  encouraging  since  it  suggests  at  least  that  models  are 

adequately  representing  the  short-term  response  of 
Polar  Cap  Absorption  Events.  Crutzen  et  al.  ozone  to  NO,  injections  at  higher  altitudes. 

( 1975)  noted  that  large  polar  cap  absorption  (PCA)  The  Solar  Eclipse  on  February  26,  1979.  In  ad- 

events  should  produce  significant  quantities  of  NO,  dition  to  ozone  observations,  measurements  of 

in  the  middle  and  upper  stratosphere,  especially  at  other  minor  constituents  during  a  solar  eclipse 

high  latitudes.  Indeed  the  PCA  events  of  August  could  provide  validation  of  the  short  time  constant 

1972  were  estimated  to  have  produced  several  times  chemistry  in  atmospheric  models.  Consequently,  ex- 

as  much  NO,  as  the  ambient  content  of  the  at-  periments  for  upcoming  solar  eclipses,  when 

mosphere  in  the  region  above  40  km.  In  a  com-  properly  supported  by  theoretical  analysis,  could 

parison  with  Nimbus  4  ozone  data.  Heath  el  al.  contribute  significantly  to  our  understanding  of  at- 

(1977)  found  that  the  agreement  between  predicted  mospheric  chemistry.  In  fact,  given  the  proper  data 

ozone  change  and  observation  was  quite  good  north  on  trace  species  concentrations  during  an  eclipse, 

of  75°  latitude  and  above  4  mb  (the  model  under-  such  measurements  could  provide  a  direct 

predicted  the  perturbation  by  roughly  30%).  Since  demonstration  that  the  currently  proposed  major 

those  calculations  were  made,  several  important  reactions  of  NO„  HO„  and  CIX  species  are  indeed 

changes  in  model  chemistry  have  occurred.  concurrently  functioning  in  the  stratosphere  in  the 

However,  these  changes  have  had  only  a  modest  ef-  manner  suggested  by  laboratory  chemistry, 

feet  on  the  sensitivity  of  the  model  to  large  injec-  However,  the  measurements  y  ield  little  information 

tions  of  NO,  at  high  altitude  and  high  latitude.  about  slower  processes  that  can  have  large  effects 

There  are  two  recent  recalculations  of  the  effect  of  on  model  sensitivity  to  perturbations, 

the  August  1972  PCA  event  in  the  literature.  Fabian  We  have  examined  theoretically  (Wuebbles  and 

et  al.  (1979)  found  that  their  computed  ozone  Chang.  1979)  the  expected  effect  of  a  solar  eclipse 

change  was  in  excellent  agreement  with  observation  on  stratospheric  minor  constituents.  Primary 

if  more  recent  estimates  of  the  NO  production  per  emphasis  was  given  to  the  total  eclipse  that  occurred 


1956  1958  1960  1962  1964  1966  1968  1970 

Year 

FIG.  2.  C  omputed  change  in  local  ozone  concentration  at  various  altitudes  resulting  from  atmospheric  nuclear 
testing  during  I95fl's  and  I960’s. 
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over  North  America  on  February  26,  1979.  Varia¬ 
tions  similar  to  those  computed  for  this  particular 
case  should  be  expected  for  other  total  eclipses. 
Totality  was  the  longest  (approximately  3  minutes) 
at  50°N  latitude  for  the  February  1979  eclipse. 

While  ozone,  tropospheric  water  vapor,  and 
temperature  were  held  fixed  in  the  model  for  most 
of  our  eclipse  calculations,  the  model  was  also  run 
with  calculated  ozone  to  examine  the  expected 
response  of  ozone  to  an  eclipse.  Our  analysis  has 
shown  that  fixing  the  ozone  distribution  does  not 
significantly  affect  the  temporal  variations 
calculated  for  other  species  during  the  eclipse,  Solar 
flux  variations  during  the  eclipse  were  based  on 
Hunt  (1965). 

Those  species  having  chemical  lifetimes  less 
than  a  few  hours  are  expected  to  vary  significantly 
from  normal  diurnal  behavior  during  a  solar 
eclipse.  Local  concentrations  of  the  species  could  be 
quite  variable,  and  therefore  we  should  focus  on 
relative  changes  rather  than  absolute  magnitudes. 

The  model-calculated  response  of  ozone  during 
an  eclipse  essentially  agreed  with  Hunt  (1965).  A 
significant  increase  in  03  is  to  be  expected  in  the  up¬ 
per  stratosphere  and  in  the  mesosphere  as  a  result  of 
the  conversion  of  0(3P)  to  ozone  through  the  reac¬ 
tion  0(3P)  +  02  +  M  -  03  +  M  accompanied  by 
decreased  photolysis  of  02  and  O3.  The  maximum 
increase  in  03,  found  at  the  end  of  totality,  was 
computed  to  be  15%  and  45%  at  50  and  55  km, 
respectively.  Larger  fractional  changes  are  expected 
in  the  mesosphere.  Since  most  of  the  atmospheric 
ozone  is  at  lower  altitudes  in  the  stratosphere,  an  in¬ 
significant  change  in  the  total  ozone  column  is  ex¬ 
pected.  Significant  changes  were  also  predicted  to 
occur  in  NO,  N02,  Cl,  CIO,  OH,  and  H02  concen¬ 
trations. 

The  results  of  this  study  suggest  that  significant 
and  detectable  variations  are  expected  for  some  of 
the  important  stratospheric  minor  constituents  dur¬ 
ing  a  solar  eclipse.  Such  observations,  particularly 
simultaneous  observations  of  trace  species,  would 
demonstrate  clearly  the  simultaneous  functioning  of 
the  various  important  photochemical  catalytic  cy¬ 
cles  in  the  stratosphere. 

The  model  calculations  were  completed  prior 
to  the  eclipse.  Measurements  during  the  solar 
eclipse  of  February  26,  1979,  were  made  with 
NASA's  aircraft.  This  plane  carried  instruments  to 
measure  NO,  03,  and  temperature  at  20  km  (instru¬ 
ment  of  Starr,  Craig,  and  others  of  NASA-Ames) 
and  to  measure  the  NO  and  N02  columns  above 


20  km  ( by  David  Murcray  of  University  of  Denver). 
Preliminary  results  by  Starr  and  Craig  show  no 
change  in  03  at  20  km  during  the  eclipse  (as  was 
predicted),  and  they  show  excellent  agreement  with 
the  theoretically  expected  change  in  NO.  A  detailed 
comparison  between  this  measurement  and  model 
calculations  is  being  made.  In  any  case,  the  results 
demonstrate  that  the  reactions  NO  +  03  and  N02 
+  he  occur  in  the  stratosphere  at  rates  similar  to 
those  computed  theoretically. 

2.2  POTENTIAL  CHANGES  IN 
OZONE  CAUSED  BY 
AIRCRAFT  EMISSIONS 

Emission  Indexes  and  Fleet  Projections 

We  consider  the  potential  effects  on  at¬ 
mospheric  ozone  of  several  different  aircraft  emis¬ 
sions  scenarios.  These  emission  scenarios  were 
developed  for  three  basically  different  applications: 
(I)  the  projected  1990  fleet,  (2)  a  commercially 
viable  fleet  of  supersonic  transports,  and  (3)  a  com¬ 
mercially  viable  fleet  of  hypersonic  transports.  In  all 
of  the  perturbation  calculations  discussed,  the  “am¬ 
bient”  or  “unperturbed"  atmospheric  conditions 
refer  to  the  model-calculated  initial  state  rather  than 
to  atmospheric  measurements. 

1990  Fleets.  We  developed  emission  profiles  for 
both  high  and  low  projected  1990  aircraft  fleet  sizes. 
In  the  case  of  the  projected  1990  high  fleet,  we  adop¬ 
ted  the  emissions  factors  and  fleet  projections  used 
by  Oliver  el  al.  (1977)  in  their  Table  2.33.  The  pro¬ 
jection  used  by  Oliver  et  al.  (1977)  was  based  on  A. 
D.  Little.  Inc.  (1976)  with  corrections  to  the  mean 
flight  altitude  of  the  projected  SST  fleet  of 
142  Concordes  and  TU  I44's  and  an  emissions  fac¬ 
tor  four  times  larger  than  used  by  A.  D.  Little.  Inc. 
(1976)  for  CF6  engines  (which  Oliver  et  al..  1977. 
treated  as  having  the  same  emissions 
indexes  as  JT9D  engines).  The  NO,  emission  index 
for  SST's  was  assumed  to  be  20  g  N02/kg  fuel. 

In  converting  the  projected  emissions  in 
Table  2.33  of  Oliver  et  al.  (1977)  to  a  format  com¬ 
patible  with  the  one-dimensional  model,  we  treated 
the  model  as  a  Northern  Hemisphere  model, 
summed  emissions  over  all  latitudes  between  zero 
and  90°N  at  a  given  altitude,  and  converted  kg/year 
at  each  altitude  to  molecules  cnr1  s-1  over  a  l-kni- 
thick  layer  centered  at  even-kilometer  altitudes. 
Table  1  gives  the  emission  rates  used  in  our  calcula¬ 
tion. 
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TABLE  I.  Projected  1990  aircraft  emissions  of  N0x 
[Oliver  et  al.  (1977)  high  estimate j;  low  estimate  of 
emission  rates  equals  0.633  times  high  estimate. 


Injection 

altitude 

(km) 

NO,.  injection  rate 
total  fleet 

(molecules  cm'V*) 

NOx  injection  rate 
subsonic  only 
(molecules  cm 

6 

90 

90 

7 

179 

179 

8 

265 

265 

9 

665 

665 

10 

1167 

1167 

tl 

1161 

1 161 

12 

520 

520 

13 

75 

75 

14 

I8a 

12 

15 

18 

— 

16 

33 

— 

17 

43 

— 

18 

29 

— 

19 

8 

— 

Emissions  from  the  SST  fleet  are  included  at  14  km  and 
above. 


We  also  examined  the  effects  of  the  projected 
subsonic  fleet  without  the  SST  contribution  using 
the  emissions  given  in  Table  I  for  the  subsonic  fleet 
only.  We  generated  a  1990  low  estimate  by  multiply¬ 
ing  the  1990  high  estimate  emissions  by  0.633,  the 
ratio  of  the  1990  low  estimate  fuel  usage  to  the  1990 
high  estimate  fuel  usage  given  in  Table  2.25  of 
Oliver  et  al.  (1977). 

SST  Fleets  of  Commercially  Viable  Size,  In 

treating  hypothetical  SST  fleets  of  commercially 
viable  magnitude,  we  have  considered  a  fleet 
emitting  1000  molecules  of  NO,  cm-As-l  and 
1.77  X  I05  molecules  of  HiO  cm- -V1  over  a  l-km- 
thick  layer  in  the  Northern  Hemisphere,  which  is 
equivalent  to  6.2  X  10s  kg  of  NO,  (as  NO:)  and 
4.3  X  I010  kg  of  HiO  per  year.  These  emission  rates 
correspond  to  a  fuel  usage  of  3.5  X  I0ln  kg  yH  by 
SST's,  assuming  the  emission  indexes  are  18  g 
NOj/kg  fuel  and  1.25  kg  HiO/kg  fuel.  The  NOx 
emissions  index  for  currently  realizable  SST  engines 
has  been  estimated  to  be  as  high  as  22  g  NOi/kg 
fuel  or  as  low  as  15.6  g  NOi/kg  fuel.  A  still  un¬ 
resolved  discrepancy  exists  between  spectroscopic 
and  probe  sampling  methods  of  determining  NO,  in 
aircraft  exhaust  (Oliver  et  al.,  1977).  If  the  spec¬ 
troscopic  analyses  are  correct,  the  above-cited  emis¬ 
sion  index  should  be  increased  by  a  factor  of  2  to  3 
for  SST  operations. 


Future  technologies  may  be  capable  of  reduc¬ 
ing  the  NO,  emission  index  by  several  fold.  A  three¬ 
fold  reduction  in  the  NO,  emission  index  is  pro¬ 
jected  for  some  existing  design  concepts  (Popoff  et 
al.,  1978).  We  consider  NO,  emission  indexes  of 
18  g  NOj/kg  fuel  for  current  technology  and  6g 
NOj/'kg  fuel  for  future  technology. 

The  number  of  SST’s  corresponding  to  a  fuel 
usage  of  3.5  X  1 0 ,0  kg  yr'1  is  difficult  to  determine 
on  an  absolute  scale.  The  number  of  SST’s 
corresponding  to  a  given  fuel  usage  is  inversely 
proportional  to  the  expected  hours/day  of  flight 
and  the  fuel  usage  per  hour.  Estimates  of  the  hours 
per  day  of  flight  at  cruise  altitude  have  varied  from 
4.4  to  7.5  hr.  Estimates  of  fuel  usage  have  ranged 
from  16,800  to  19,100  kg/hr  for  the  Concorde;  from 
52,000  to  60,000  kg/hr  for  the  B2707;  and  as  low  as 
44,000  kg/hr  for  a  hypothetical  advanced  design 
similar  to  the  B2707  (Popoff  et  al.,  1978).  Thus, 
more  than  a  two-fold  variation  in  engine  emissions 
from  a  projected  SST  fleet  is  possible  based  on  dif¬ 
ferent  estimates  of  aircraft  operation  and  fuel  con¬ 
sumption  rate.  A  fuel  usage  of  3.5  X  10'°  kg  yr1. 
therefore,  corresponds  to  about  750-1000  Con¬ 
cordes,  or  i80-660  advanced  SST's.  Thus,  our  stan¬ 
dard  emissions  rates  are  compatible  with  a  large 
fleet  of  SST’s. 

Hypersonic  Transport  Emissions.  In  addition  to 
the  effects  of  subsonic  and  supersonic  fleets,  w;e  pro¬ 
jected  a  range  of  emissions  and  environmental  im 
pacts  of  a  possible  hypersonic  transport  (HST)  fleet 
For  this  purpose,  we  adopted  the  emissions  indexes 
given  in  CIAP  Monograph  2  (1975),  Chapter  6  for  a 
fleet  of  research  HST’s  of  gross  takeoff  weight 
2.27  X  10'  kg  operating  at  Mach  8  and  cruise 
altitude  of  “36.68"  km.  According  to  data  in  CIAP 
Monograph  2.  the  hypothesized  HST  operations 
would  consist  of  an  initial  acceleration  to  hyper¬ 
sonic  speeds  (and  cruise  altitude)  using  rocket 
engines,  followed  by  about  five  minutes  of 
SCRAMjet  operation,  followed  by  a  cruise-mode 
during  w  hich  most  emissions  would  consist  of  liquid 
hydrogen  boil-off  for  cooling.  The  SCRAMjet 
would  use  670  kg/flight  of  liquid  Hi;  the  rocket 
engines  would  use  2450  kg/flight  of  liquid  Hi,  and 
14,650  kg/flight  of  liquid  oxygen  (CIAP 
Monograph  2,  p.  6-7);  but  it  would  only  release  a 
total  of  1650  kg/flight  of  water  vapor  and  615 
kg/ flight  of  Hi  (CIAP  Monograph  2,  p.  6-8).  Thus, 
the  fuel  use  and  emissions  estimates  given  in  CIAP 
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Monograph  2  for  rocket  operations  are  incompati¬ 
ble,  with  not  enough  mass  being  emitted.  This  dif¬ 
ficulty  would  be  resolved  by  increasing  the  water- 
vapor  emissions  by  a  factor  of  10.  Cooling  at  cruise 
altitude  would  consume  an  additional  309  kg  of  liq¬ 
uid  hydrogen  per  flight. 

The  emissions  consist  of  HiO,  Hi,  H,  OH,  NO, 
and  O  as  given  in  Table  2.  In  analyzing  the  effects  of 
HST  emissions,  only  cooling  and  SCRAMjet  emis¬ 
sions  were  treated,  and  these  were  based  on  the  data 
presented  in  Table  6-9  of  CIAP  Monograph  2. 
Rocket  emissions  were  neglected  since  much  of  the 
rocket  emissions  would  occur  in  the  troposphere.  In 
order  to  approximate  a  commercially  viable  fleet, 
we  increased  the  emissions  to  correspond  to  an  HST 
fleet  with  1000  flights /day.  Because  even  these  emis¬ 
sions  are  relatively  small,  and  the  CIAP- 
Monograph-2-based  emissions  estimates  may  not 
apply  to  commercial-scale  fleet  operations  (they  are 
based  on  research  flights),  we  also  examined  the  ef¬ 
fects  of  a  10-fold  larger  emission,  perhaps  inter¬ 
pretable  as  longer  flights  or  flights  of  heavier  air¬ 
craft.  The  emissions  of  H20,  Hi,  NO,  and  OH  were 
treated  simply  as  a  source  of  the  molecules  in  ques¬ 
tion.  Emissions  of  H-atoms  (the  code  calculates 
H-atoms  as  an  equilibrium  species)  were  treated  as 
emissions  of  an  equivalent  number  of  HOi 


molecules,  while  the  small  emissions  of  O  atoms 
were  ignored. 

Because  NOx  emissions  are  responsible  for  the 
major  portion  of  the  predicted  change  in  total 
ozone,  it  might  be  attractive  to  use  more  than  the 
stoichiometric  amount  of  OH  in  the  SCRAMjets  in 
order  to  reduce  NOx  emissions.  According  to  CIAP 
Monograph  2,  operating  at  a  stoichiometric  ratio  of 
1.5  instead  of  1  would  reduce  NOx  emissions  by 
more  than  a  factor  of  3  while  increasing  HiO  and 
Hi  emissions  by  factors  of  about  1.1  and  3,  respec¬ 
tively.  The  net  effect  of  these  changes  in  emissions 
should  be  a  smaller  ozone  perturbation. 

We  wish  to  emphasize  that  we  performed  no  in¬ 
dependent  calculations  of  HST  emissions  indexes, 
and  that  the  fleet  projections  are  no  more  than 
crude  parametric  estimates  of  the  approximate  level 
of  activity  that  might  be  associated  with  a  mature, 
commercially  successful  fleet  of  HST's.  If  specific 
aircraft  designs  were  to  be  proposed,  the  emissions 
should  be  independently  projected. 

Aircraft  Assessment  Results 

The  LLL  one-dimensional  transport-kinetics 
model  has  been  used  to  assess  the  potential  chemical 
effects  of  aircraft  engine  emissions  in  the 
troposphere  and  the  stratosphere.  Calculations  were 


TABLE  2.  HST  emissions  on  a  per-flight  basis.  Data  from  CIAP  Monograph  2,  Ch.  6,  assuming  cruising  speed  of 
MACH  8,  cruise  altitude  of  36.68  km,  and  H2:02  ratio  of  1 . 


Emissions  per  flight 

Altitude 

(km) 

Emission  rate  for  1000 
flights/daVj 
(molecules  cm  s  r 

(kg) 

(moles) 

SCRAMJet 

hto 

5705.0 

317.000 

36.68 

8504 

”2 

136.7 

68.350 

36.68 

1820 

H 

13.1 

13.100 

36.68 

350 

OH 

132.3 

7,800 

36.68 

210 

NO 

137.4 

4,500 

36.68 

123 

O  34.4 

Rocket  (not  included  in  calculation) 

2,150 

36.68 

58 

H2Ob 

16,500 

917,000 

0-36  mostly 
tropospheric 

— 

»2 

615 

307,500 

0-36 

— 

Cooling  in  flight 

H2 

308.5 

154,200 

36 

4130 

g 

Emissions  in  a  I  -km-thick  layer. 

**H20  emissions  from  rocks  have  been  increased  by  a  factor  of  10  to  permit  mass  balance. 
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made  using  projected  fleet  sizes  for  subsonic  and 
supersonic  aircraft  as  well  as  for  a  range  of  emission 
rates  at  particular  injection  altitudes.  In  order  to 
determine  the  effect  of  uncertainties  in  various  input 
parameters  on  the  assessment  calculations,  several 
sensitivity  studies  were  also  performed  and  are 
reported  in  the  following  section, 

Subsonic  and  Supersonic  Aircraft  Fleets.  The 
altitude  of  injection  has  a  significant  effect  on  the 
computed  change  in  ozone  because  of  the  increase 
in  residence  time  with  altitude  and  because  of  the 
variation  with  altitude  of  the  dominant  chemical 
reactions  and  cycles.  Table  3  shows  the  effect  on 
total  ozone  at  steady  state  of  the  same  NOx  injection 
rate  at  different  injection  altitudes.  An  injection  rate 
of  1000  molecules  cm-^s-1  over  a  1-km-thick  layer 
was  used  for  NO,,  and  there  was  no  injection  of 

h2o. 

Injections  of  NO,  over  the  altitude  range  of  7 
to  20  km  resulted  in  a  net  increase  in  total  ozone  in 
each  case.  The  change  in  total  ozone  increased  with 
increasing  injection  altitude  for  the  troposphere  and 
lower  stratosphere.  Throughout  this  region  HO, 
chemistry  is  the  dominant  chemical  destruction 
process  for  ozone.  When  HO,  chemistry  is  more  ef¬ 
ficient  than  NO,  catalytic  destruction  of  ozone,  in¬ 
jections  of  NO,  lead  to  a  net  increase  in  odd  oxygen 
production  through  the  reaction  sequences: 

NO  +  N02  -*  N02  +  OH 
NO,  h-  NO  +  O 

... _  0) 

net  H02  -  OH  +  O 


TABLE  3.  Change  in  total  ozone  resulting  from  NO 

l  _ ■  * 

injection  of  1000  molecules  cm  s  distributed  over 
a  1-km-thick  layer  centered  at  injection  altitude. 
Calculations  were  made  using  1979a  chemistry  (Ap¬ 
pendix  A). 


Injection  altitude 
(km) 

Change  in  total  ozone 

m 

7 

0.25 

9 

0.40 

II 

0.64 

13 

0.84 

17 

1.34 

20 

1.31 

N02  +  OH  “  HN03 
OH  +  HN03  ->  H20  +  N03 
no3  h-£  N02  +  O 


net  20H-H20  +  0 

The  reaction  sequence  in  Eq.  (1)  acts  to  shift  HO, 
more  toward  OH.  thereby  increasing  the  conversion 
of  OH  and  HO;  to  H;0  by  the  reaction  HO;  + 
OH  ->  HjO  +  O;.  Reaction  sequence  in  Eq.  (2)  acts 
to  reduce  HO,  as  OH  and  HO;  are  converted  to 
H 20  and  HNO3.  Conversion  to  H;0,  and  to  a  lesser 
extent  HNO3,  acts  as  a  sink  for  HO,,  since  the 
chemical  lifetimes  of  these  reservoir  species  can  be 
long  compared  to  the  time  for  transport  into  the 
troposphere.  Injections  of  NO,  act  to  increase  the 
rates  of  the  above  reactions.  The  net  effect  of  reac¬ 
tion  sequences  ( 1 )  and  (2)  is  a  reduction  of  HO,  and 
an  increase  in  odd  oxygen  production,  both  of 
which  contribute  to  an  increase  in  ozone.  Of  course, 
the  odd-hydrogen  species  may  have  been  generated 
through  reaction  of  O('D)  with  H;0.  If  this  is  the 
case,  the  odd  oxygen  produced  via  (1)  or  (2)  only 
regenerates  that  used  in  the  initiation  of  a  catalytic 
cycle.  It  is  important,  nonetheless,  because  these 
sequences  compete  for  HO,  with  odd-oxygen 
destroying  sequences. 

Increases  in  NO,  concentration  in  the  roper 
stratosphere  lead  to  a  reduction  in  ozone  concentra¬ 
tion  in  this  region  due  to  the  greater  importance  of 
the  NO,  catalyzed  ozone  destruction  at  these 
altitudes.  The  20-km  injection  of  NO,  resulted  in  a 
smaller  net  increase  in  total  ozone  than  did  the  17- 
km  injection  case,  since  more  NO,  reached  the  up¬ 
per  stratosphere  causing  a  greater  reduction  in 
ozone  concentration  at  high  altitudes  relative  to  the 
increase  at  lower  altitudes.  This  effect  was  found  to 
be  sensitive  to  the  shape  of  the  K,  profile  used, 
however. 

The  standard  model  uses  fixed-value  boundary 
conditions  for  all  species  (except  CCI4)  at  the 
ground.  Using  flux  boundary  conditions  for  the 
species  N;0,  CH4,  CH3CI,  and  HCI.  resulted  in 
slightly  larger  increases  in  total  ozone  (approx¬ 
imately  0. 1  jn  absolute  change)  when  compared  to 
the  fixed  boundary  condition  case.  Consequently, 
the  choice  of  boundary  condition  has  only  a  very 
small  effect  on  the  assessment  results.  Fixed-value 


boundary  conditions  arc  used  in  the  model  calcula¬ 
tions  unless  otherwise  noted. 

The  model  was  tested  for  interference  effects  by 
comparing  the  effect  of  simultaneous  injections  at 
the  various  injection  altitudes  given  in  Table  3  with 
the  sum  of  the  changes  in  ozone  for  individual  injec¬ 
tion  altitudes.  Table  4  shows  that  there  is  a  destruc¬ 
tive  interference.  That  is,  the  simultaneous  injec¬ 
tions  resulted  in  less  change  in  total  ozone  than  the 
sum  of  the  changes  caused  by  individual  injections. 
This  implies  that  summing  the  results  of  separate 
calculations  for  subsonic  and  supersonic  fleets  will 
not  give  the  same  answer  as  considering  the  effect  of 
the  two  fleets  simultaneously. 

The  effect  of  NOx  emissions  by  subsonic  and 
supersonic  fleets  projected  for  1990  are  given  in 
Table  5.  Using  the  high  estimate  of  the  fleet  sizes, 
the  effect  of  subsonic  and  supersonic  aircraft  com¬ 
bined  is  estimated  to  be  an  increase  in  total  ozone  of 
2.01%.  This  number  represents  the  steady  state 
change  in  ozone  due  to  constant  NOx  emission  rates 
at  the  rates  estimated  for  1990.  The  subsonic  fleet 
alone  (injection  altitudes  up  to  14  km)  is  estimated 
to  cause  an  increase  in  total  ozone  of  1.86%.  The 
small  effect  of  the  supersonic  fleet  is  due  primarily 
to  the  small  fleet  size  projected  for  1990.  Using  the 
lower  estimate  for  fleet  emissions,  the  corres¬ 
ponding  numbers  for  AO,  are  1.39%  for  the  com¬ 
bined  fleets  and  1.29%  for  the  subsonic  fleet  alone. 

It  should  be  noted  that  the  effects  of  the  ex¬ 
isting  1979  subsonic  fleet  are  of  the  order  of  a  0.5% 
increase  in  total  ozone.  This  effect  is  not  negligible 
when  compared  to  the  estimated  present-day  effects 
of  CFM’s.  Thus,  an  accurate  assessment  of  the  ef¬ 
fects  of  subsonic  aircraft  will  be  important  in  inter¬ 
preting  ozone  data  for  trends  due  to  other 
anthropogenic  influences. 

The  changes  in  ozone  concentration  for  the 
high  and  low  estimates  of  subsonic  and  supersonic 
fleet  emissions  are  shown  in  absolute  concentration 
in  Fig.  3.  The  largest  absolute  increase  in  ozone 
concentration  occurs  near  12  km,  and  the  region  of 
increasing  ozone  extends  up  to  about  26  km.  Be- 


TABLE  4.  Synergism  test  using  fixed  boundary  con¬ 
ditions  and  1979a  chemistry. 

Caw 

A03(%) 

Simultaneous  injections 

3.68 

Sum  of  individual  injections 

4.78 

TABLE  5.  Computed  change  in  total  ozone  at  steady 
state  using  NOx  emission  estimates  for  1 990  subsonic 
and  supersonic  fleets  (1979b  chemistry). 


Change  in  total  ozone  (%) 

Subsonic  and 

Subsonic 

Case 

supersonic 

only 

Oliver  et  al.  ( 1977)  high 

1990  fleet  estimate 

2.01 

1.86 

Low  1990  fleet  estimate 

1.39 

1.29 

tween  26  and  39  km  there  is  a  small  decrease  in 
ozone  concentration. 

It  should  be  recognized  that  the  bulk  of  the 
ozone  change  calculated  for  the  subsonic  fleet  oc¬ 
curs  in  the  upper  troposphere.  As  a  result,  it  is  sen¬ 
sitive  to  the  treatment  of  such  poorly  understood 
phenomena  as  wet  and  dry  removal  processes  and 
surface  boundary  conditions  for  rarely  measured 
species.  The  tropospheric  pressures  differ  from  the 
pressure  conditions  used  in  most  direct  measure¬ 
ments  of  chemical  rate  coefficients  (to  a  greater 
degree  than  stratosphere  pressures).  Heterogeneous 
reactions,  which  might  have  a  significant  effect  on 
the  calculations,  are  not  included,  and  the  one¬ 
dimensional  treatment  is  not  as  good  an  approxima¬ 
tion  for  the  troposphere  as  it  is  for  the  stratosphere 
because  of  the  strong  latitudinal  and  longitudinal 
gradients  in  many  trace  species.  Thus,  the  results  for 
the  1990  fleet  estimates  are  suggestive  but  are  not 
definitive. 

Assessments  of  potential  changes  in  ozone  due 
to  future  large  fleets  of  supersonic  transports  have 
focused  on  injection  altitudes  of  17  and  20  km.  In 
updating  these  assessments,  we  have  assumed  an 
emission  index  of  18  g/kg  fuel  for  NOx  and  1250 
g/kg  fuel  for  H,0.  These  emission  indexes  are  based 
on  current  engine  technology.  It  may  be  possible 
through  future  advances  in  technology  that  the  NOx 
emission  index  could  be  reduced  to  6  kg  (one-third 
of  present  value).  Calculations  of  AO,  given  a  con¬ 
stant  fuel  consumption  rate  at  cruise  altitude 
(3.5  X  I0in  kg/yr  in  a  hemispheric  shell)  are  shown 
in  Table  6  for  different  emission  indexes.  The  NOx 
emission  rate  of  1000  molecules  ctrrV1,  which 
roughly  estimates  the  emissions  from  a  commer¬ 
cially  viable  fleet,  was  chosen  as  a  reference  case. 
The  equivalent  annual  hemispheric  injection  rates 
of  NOx  and  H,0  are  given  in  Table  7. 
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FIG.  3.  Change  in  ozone  concentration  expressed  as  molecules  cm'3  ( 1979b  chemistry)  computed  for  1990  high 
and  low  fleet  estimates. 


TABLE  6.  Change  in  total  ozone  due  to  NOx  and  HjO  emissions  distributed  over  a  1-km-thick  layer  centered  at 
injection  altitude. 


Injection  altitude 
(km) 

Injection  rate  (molecules  cm'^s'S 

Change  in  total  ozone  (%) 

N°x 

h2o 

1979a  Chemistry 

1979b  Chemistry 

17 

1000 

0 

1.34 

1.25 

17 

1000 

177,000 

1. 18 

— 

17 

333 

177,000 

0.32 

— 

20 

1000 

0 

1.31 

1.14 

20 

1000 

177,000 

0.91 

_ 

20 

333 

177,000 

0.07 

— 

For  each  of  the  perturbations  considered,  there 
was  an  increase  in  total  ozone.  When  the  water 
vapor  injection  is  included  with  the  NOx  injection, 
there  is  less  of  an  increase  in  total  ozone  because  of 
the  additional  ozone  destruction  caused  by  HOx 
which  is  produced  from  dissociation  of  FGO.  Even 
for  the  advanced  technology  case  (reduced  NOx 
emission),  the  increase  in  total  ozone  due  to  the 
NOx  injection  is  greater  than  the  reduction  due  to 
the  H20  injection,  the  net  effect  being  a  very  small 
increase  in  total  ozone. 


TABLE  7.  Equivalent  annua)  hemispheric  injection 
rates  for  various  NOx  and  HjO  emissions. 

Injection  rate 

Species 

—3  -I  a 

(molecules  cm  s  )  (kg/yr) 

NOx 

333  2,1  x  108 

NOx 

1.000  6.2  X  108 

HjO 

177,000  4.3  XI010 

hemispheric  injection  assuming  uniform  injection  over  a 
1-km-thick  layer. 
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Hypersonic  Transport  Fleet.  The  hypothesized 
hypersonic  transport  lleet  making  1000  flights  per 
day  was  calculated  to  cause  a  -0.218%  change  in  the 
ozone  column.  Emissions  from  a  lleet  ten  times 
larger  were  calculated  to  cause  a  -2.06%  change  in 
the  ozone  column.  The  NO,  emissions  from  the  lat¬ 
ter  fleet  by  themselves  caused  an  ozone  change  of 
-2.13%. 

Figure  4  shows  the  local  percentage  changes  in 
ozone  vs  altitude  for  the  10,000  (light  per  day  HST 
emissions  and  the  NO,  emissions  alone  from  10,000 
flights  per  da>  of  HST's.  (Because  AO;  is  roughly 


linear  in  emissions  rate  over  the  range  considered, 
large  fleets  were  chosen  to  avoid  comparing  very 
small  local  changes  that  might  contain  significant 
numerical  noise.)  As  is  evident,  the  NO,  emissions 
dominate  the  ozone  changes  below  approximately 
40  km,  while  HiO  and  H;  emissions  are  responsible 
for  most  of  the  ozone  change  at  higher  altitudes. 
Because  most  of  the  ozone  column  is  below  40  km, 
the  integral  column  change  is  largely  a  result  of  the 
NO,  emissions. 

Effect  of  V  ariations  in  Background  (  IX.  The  ef 

feet  of  aircraft  engine  emissions  on  stratospheric 


-20  -15  -10  -5  0  5  10  15 


Change  in  ozone  concentration  —  % 

FIG.  4.  Change  in  ozone  concentration  caused  by  HST  emissions  for  a  fleet  with  10,000  flights  per  day  ( 107% 
chemistry ). 
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composition  is  coupled  with  other  stratospheric  per¬ 
turbations.  As  chlorofluoromethanes  (CFM’s)  con¬ 
tinue  to  be  released  to  the  atmosphere,  there  will  be 
a  gradual  increase  in  the  background  concentration 
of  C1X  (Cl  +  CIO  +  HC1).  Model  calculations  in¬ 
dicate  a  present  day  stratospheric  CIX  concentra¬ 
tion  of  more  than  1  ppbv.  Assuming  CFM  releases 
at  the  1976  release  rate  indefinitely  into  the  future, 
the  background  CIX  concentration  is  projected  to 
increase  to  2  ppbv  around  1990  and  to  4  ppbv 
around  2020.  To  assess  the  effect  of  changes  in 
background  CIX  concentration  on  the  model  sen¬ 
sitivity,  several  NOx  perturbation  calculations  were 
repeated  using  a  range  of  CIX  concentrations.  The 
results  are  shown  in  Table  8.  In  each  case  the  CIX 
concentration  was  held  constant.  So  the  change  in 
total  ozone  is  the  steady  state  value  for  a  constant 
injection  rate  (i.e.,  the  AO,  value  does  not  pertain  to 
any  particular  year)  relative  to  a  natural  atmosphere 
with  the  same  CIX  concentration. 

*  Table  8  shows  that  the  change  (increase)  in 

total  ozone  resulting  from  an  NOx  injection  in¬ 
creases  as  the  CIX  concentration  increases.  The 

/  reason  for  the  increase  in  AO,  is  related  to  the 

mechanism  by  which  an  NOx  injection  leads  to  an 
increase  in  total  ozone.  The  increase  in  total  ozone 
comes  about  because  for  these  magnitudes  of  NOx 
injections,  NOx  interferes  with  the  more  efficient 
ozone  destruction  by  HOx  and  CIX.  Increasing  the 
background  CIX  concentration  means  that  there  is 
more  ozone  destruction  by  CIX  in  the  ambient  case, 
so  there  is  more  that  can  be  interfered  with  by  the 
injected  NOx. 

Effects  of  I'ncertainties  in  Chemical  Rate  Coef¬ 
ficients.  Probably  the  best  method  for  assessing  the 
sensitivity  of  model  predictions  to  errors  in  the  in¬ 
puts  that  describe  fairly  well  understood  quantities 
(e.g.,  rate  constants  for  which  experimental  preci¬ 
sion  is  the  dominant  source  of  error)  is  a  Monte 


Carlo  calculation  (e.g.,  Stolarski  el  al.,  1978).  As 
vet,  no  such  calculation  has  been  carried  out  for  a 
perturbation  resulting  from  aircraft  operations.  In  I 

part  this  reflects  the  hypothesis  that  prediction  error  ' 

is  dominated  by  systematic  error  in  poorly  un¬ 
derstood  reactions. 

Even  a  full  sensitivity  analysis  (like  that  of  j 

Butler,  1979)  has  not  yet  been  performed.  The  most 
recent  sensitivity  study  on  the  NOx  system  was  a 
partial  sensitivity  analysis  carried  out  in  1975  by 
Duewer  el  al.  (1976.  1977).  The  major  conclusions 
of  that  study  were  that  large  prediction  errors  were 
possible,  and  that  a  few  reactions  dominated  the 
error.  Subsequent  events  have  shown  these  conclu¬ 
sions  to  be  valid.  The  large  changes  in  model  sen¬ 
sitivity  to  NOx  injection  that  have  occurred  have 
been  dominated  by  changes  in  the  rate  coefficients 
for  the  five  reactions  that  were  identified  as  plausi¬ 
ble  sources  of  error  in  the  1975  and  1977  reports.  In 
the  intervening  years,  many  rate  coefficients  have 
been  determined  or  redetermined  with  substantially 
improved  accuracy.  For  the  majority  of  the  reac¬ 
tions  in  the  model,  experimental  precision  is 
probably  a  reasonably  good  estimate  (as  well  as  a 
lower  bound)  to  the  actual  error  in  the  measured 
rate,  and  model  predictions  are  only  weakly  sen¬ 
sitive  to  modest  errors  in  most  rate  coefficients. 

Nonetheless,  the  composite  error  in  prediction 
caused  by  such  errors  could  be  significant, 
although,  since  it  is  unlikely  that  random  errors  will 
reinforce  each  other,  this  is  not  too  likely  to  be 
realized. 

Table  9  gives  the  change  in  O,  calculated  for  , 

several  perturbations  using  the  structure  of  our  ; 

current  model  but  using  the  rate  coefficients  recom-  ; 

mended  in  NBS  866  (1975)  and  NBS  513  (1978)  as  j 

well  as  the  current  results.  In  calculating  the  j 

response  for  1975  and  1977  chemistries,  we  excluded  > 

species  (such  as  HOCI  or  HN04)  that  the  LLL  j 


TABLE  8.  Effect  of  variations  in  background  CIX  on  computed  change  in  total  ozone  (1979a  chemistry). 


Injection  altitude 
(km) 

NO  injection  rate 

x  _3  _( 

(molecules  cm  s  ) 

Change  in  total  ozone  (7c) 

1.14  ppb  CIX 

1.83  ppb  CIX 

3.76  ppb  CIX 

17 

500 

0.65 

0.72 

0.89 

17 

1000 

1.22 

1.34 

1.69 

17 

2000 

2.14 

2.37 

3.03 

20 

500 

0.61 

0.75 

1.12 

20 

1000 

1.04 

1.31 

2.03 

20 

2000 

1.42 

1.91 

3.28 

17 


-*»!•*&*  i* 


TABLE  9.  Effect  of  choice  of  rate  coefficients  on  model  sensitivity  using  structure  of  current  LLL  model. 


Change  in  total  ozone  (%) 


Chemistry 

1950  ambient  (K 

(DU) 

1990  high 
subsonic 

17-km  SST3 

20-km  SST“ 

Cl  M  1978 

CFM 

steady-state 

1975  (NBS  866) 

288 

-0.31 

-3.94 

-9.12 

-0.76 

-7.5 

1977  (NBS  513) 

366 

1.43 

1.20 

-0.52 

-0.98 

-12.7 

1979a  (JPL,  1979c) 

324 

2.01 

2.57 

2.18 

-1.92 

-19.4 

1979b  (NASA,  1979) 

327 

1.86 

2.00 

1.06 

-1.30 

-14.2 

aNOx  injection  rate  of  2000  molecules  over  a  1-km-thick  layer.  Contribution  to  C1X  due  to  CFM’s  is  not  included  in 

calculations. 

bCFM’s  released  at  1976  rate. 

cRate  of  OH  +  CIO  -*  HC1  +  O2  is  set  to  zero. 


model  contemporary  with  the  rate  compendium  in 
question  did  not  include.  In  addition,  we  include  the 
results  of  calculations  made  using  the  1979a 
chemistry  with  the  rate  coefficient  for  the  reaction 
OH  +  CIO->  HCI  +  O;  set  to  zero. 

The  conclusion  to  he  drawn  from  Table  9  is 
that  over  the  last  five  years,  changes  in  evaluated 
recommendations  for  chemical  rate  coefficients 
have  resulted  in  substantial  changes  in  the  model 
predicted  response  to  perturbation,  and  that  even 
recent  recommendations  appearing  a  few  months 
apart  carry  a  significant  amount  of  uncertainty. 
(However,  the  differences  between  model  predic¬ 
tions  are  not  such  as  to  change  any  qualitative  con¬ 
clusions  if  one  restricts  one's  self  to  post- 1976 
models  with  consistent  physics.) 

A  comparison  of  the  results  in  Table  9  wi-.h  the 
results  of  LLL.  calculations  in  earlier  years  given  in 
Fig.  5  for  some  of  the  perturbations  demonstrates 
that  other  changes  in  model  structure  have  had 
quantitatively  significant  effects,  but  would  have 
had  little  effect  on  the  qualitative  conclusions  drawn 
from  model  predictions  of  aircraft  effects.  The 
results  in  Fig.  5  demonstrate  the  combined  effects  of 
the  evolution  of  our  understanding  of  stratospheric 
chemistry  and  evolution  of  the  treatment  of  physical 
phenomena  such  as  multiple  scattering  of  light,  the 
averaging  of  reaction  rates  over  diurnal  cycles,  the 
treatment  of  boundary  conditions,  and  the 
transport  parameterization.  Although  many  dif¬ 
ferent  factors  have  contributed  to  the  variation 
described  in  F  ie.  5,  the  evolution  of  model 
chemistry  has  been  the  most  important  single  fac¬ 
tor. 


The  uncertainty  limits  for  reactions  in  the 
various  evaluations  [NBS  513  (1978),  NBS  866 
(1975),  NASA  1010(1977),  JPL(I979),  NASA  1049 
(1979)]  have  generally  become  smaller  when  in¬ 
dividual  reactions  are  considered.  However,  over 
the  same  period,  new  and  often  highly  uncertain 
reactions  have  been  recognized,  and  it  is  far  from 
obvious  that  the  error  that  might  be  associated  with 
rate  coefficients  actually  included  in  models  has  got¬ 
ten  smaller  over  the  last  decade.  Moreover,  “new" 
reactions  have  been  introduced  to  the  modeling 
community  at  a  fluctuating  but  not  obviously 
decreasing  rate. 

Effect  of  Temperature  Feedback  and  Hydro¬ 
static  Adjustment.  In  the  calculations  up  to  this 
point,  the  i'.S.  Standard  A'mosphere  (1976)  tem¬ 
perature  profile  was  used,  and  temperature  was  not 
allowed  to  change.  Changes  in  the  concentrations  of 
species  that  are  radiatively  important  (either  solar 
or  longwave)  will  affect  the  temperature  profile  by 
changing  radiative  fluxes  and  heating  or  cooling 
rates.  At  a  given  pressure  level,  as  the  temperature 
changes,  the  air  density  will  also  change  as  defined 
by  the  equation  of  state.  Since  chemical  reaction 
rates  are  affected  by  changes  in  temperature  and  air 
density,  it  is  important  to  assess  the  magnitude  of 
these  effects  on  the  assessment  calculation. 

Air  density  can  be  computed  in  a  one¬ 
dimensional  model  by  assuming  hydrostatic 
equilibrium  (i.e.,  the  pressure  at  any  height  is  deter¬ 
mined  by  the  weight  of  the  column  of  air  above  that 
point,  and  the  air  density  is  determined  by  the  equa¬ 
tion  of  state  given  the  temperature).  Expressed 
mathematically , 
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FIG.  5.  Historical  evolution  of  LLL  model  calculations  of  change  in  total  ozone  due  to  an  NOx  injection  of  2000 
molecules  cm'V1  over  a  l-km-thick  layer  centered  at  either  17  km  (open  circles)  or  20  km  (solid  circles).  This  in¬ 
jection  rate  was  a  standard  for  assessment  and  comparison  purposes  and  did  not  pertain  to  any  particular  fleet  es¬ 
timate.  Values  for  1978  are  approximate. 


dP  =  -pg  dz  (3) 

and 

P  =  pRT  (4) 

where  P  is  pressure,  p  is  density,  g  is  the  acceleration 
of  gravity,  z  is  altitude.  R  is  the  gas  constant  and  T 
is  temperature.  The  temperature  is  computed  using 
a  radiative  transfer  model  assuming  radiative 
equilibrium  in  the  stratosphere  (see  Appendix  A).  In 
our  calculations,  the  temperature  profile  is  com¬ 
puted  above  13  km,  and  is  specified  below  this 
altitude.  Changes  in  surface  temperature  may  be  im¬ 
posed,  but  they  are  not  calculated  in  this  version  of 
the  model.  The  model  includes  solar  absorption  and 


longwave  exchange  by  O,,  H;0  and  CO;  along  with 
solar  absorption  by  NO;  (Luther  ei  a!.,  1977). 

We  assume  that  the  change  in  surface  tem¬ 
perature  is  negligible.  Ozone  reductions  of  up  to 
30%  due  to  NOx  injections  were  computed  to  cause 
less  than  a  0. 1  K  change  in  surface  temperature  by 
Ramanathan  et  al.  (1976).  An  increase  in  the 
stratospheric  water  vapor  mixing  ratio  of  2  ppmv  is 
estimated  to  cause  an  increase  in  surface  tem¬ 
perature  of  <0.2  K.  (Grobecker  el  al 1974.  p. 
FI25).  Increasing  the  surface  temperature  in  our 
model  by  0.2  K  causes  a  change  in  total  ozone  of 
-0.06%.  Neglecting  changes  in  surface  temperature 
of  this  magnitude  has  no  significant  effect  on  the 
results. 

In  studying  the  effect  of  temperature  feedback 
and  hydrostatic  adjustment  on  model  sensitivity,  we 
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considered  four  calculation  alternatives:  (1)  tem¬ 
perature  feedback  with  hydrostatic  adjustment.  (2) 
temperature  feedback  without  hydrostatic  adjust¬ 
ment,  (3)  constant  temperature  using  the  ambient 
computed  temperature  and  pressure  profiles,  and 
(4)  constant  temperature  using  the  U.S.  Standard 
Atmosphere  (1976)  temperature  and  pressure 
profiles.  The  results  are  presented  in  Table  10. 

The  constant  temperature  calculation  using  the 
ambient  computed  temperature  profile  (consistent 
with  hydrostatic  adjustment)  is  the  reference  case 
with  which  to  compare  the  effects  of  temperature 
feedback  and  hydrostatic  adjustment.  Table  10 
shows  that  inclusion  of  temperature  feedback  leads 
to  a  larger  increase  in  total  ozone  compared  to  the 
constant  temperature  calculation  (ambient  profile). 
When  hydrostatic  adjustment  is  included,  the  in¬ 
crease  in  total  ozone  is  even  larger.  To  understand 
these  effects,  "e  need  to  look  at  the  change  in  tem¬ 
perature  versus  height. 

For  analysis  purposes,  we  consider  the  results 
for  an  NOx  injection  rate  of  1000  molecules  cm'V1 
and  an  H^O  injection  rate  of  177.000  molecules 
crrr-V  at  20  km.  The  change  in  local  ozone  con¬ 
centration  is  shown  in  Fig.  6  and  the  change  in  the 
temperature  profile  is  shown  in  Fig.  7.  The  tem¬ 
perature  increases  below  20  km  causing  an  increase 
in  the  chemical  destruction  rate  of  ozone.  Thus, 
there  is  less  of  an  increase  in  ozone  concentration  in 
this  region  when  temperature  feedback  is  included. 
Conversely,  above  20  km  the  temperature 
decreases,  thereby  slow  ing  the  rates  of  reaction  and 
leading  to  less  ozone  reduction  in  this  region. 

In  the  case  with  hydrostatic  adjustment,  the  air 
in  the  14-20  km  region  expands  due  to  the  rise  in 


temperature.  As  the  air  expands  it  lifts  the  at¬ 
mosphere  above  so  that  a  given  pressure  level  is 
raised  in  altitude  relative  to  the  ambient  profile. 
Conversely,  where  the  temperature  decreases,  there 
is  contraction. 

The  results  with  temperature  feedback  and 
hydrostatic  adjustment  represent  the  most  complete 
model  calculations  in  terms  of  physical  processes 
and  feedback  mechanisms.  The  changes  in  ozone 
concentration  computed  with  these  processes  in¬ 
cluded  are  shown  in  Ffig.  8  for  various  injection 
rates  of  NO^  and  HsO  at  an  injection  altitude  of 
17  km.  The  results  for  an  NC\  injection  rate  of  1000 
molecules  cnrV 1  with  and  without  a  simultaneous 
HsO  injection  are  nearly  identical  below  26  km.  The 
injection  of  HsO  only  leads  to  noticeable  differences 
in  the  computed  change  in  ozone  concentration 
above  26  km. 

Given  the  same  (HO  injection  rate,  the  results 
for  NOv  injection  rates  of  1000  or  333  molecules 
cnr-V1  are  quite  different.  The  lower  NOv  injection 
rale  causes  much  less  change  in  ozone  concentration 
in  the  middle  and  lower  stratosphere;  the  relative 
magnitudes  of  the  change  in  ozone  concentration 
being  roughly  proportional  to  the  NC\  injection 
rate.  This  result  indicates  that  the  NC\  injection  rale 
is  a  much  more  important  (and  sensitise)  parameter 
than  the  HsO  injection  rate  in  estimating  the  effects 
of  projected  aircraft  fleets. 

Another  important  conclusion  from  Table  10  is 
that  the  choice  of  temperature  profile  in  a  constant 
temperature  calculation  om  have  a  significant  effect 
on  model  sensitivity.  The  difference  between  the 
AO}  values  using  the  ambient  or  the  I  S  Standard 
Atmosphere  (1976)  temperature  profile  is  the  same 


TABLE  10.  Effect  of  temperature  feedback  and  hydrostatic  adjustment  on  model  sensitivity  (1979a  chemistry). 

Injection 

altitude 

(km) 

Injection  rate 
(molecules  cm  ^  s  * ) 

Change  in  total  ozone  (r'<) 

Temperature  feedback 
with  hydrostatic 
adjustment 

Temperature  feedback 
without  hydrostatic 
adjustment 

Constant 

ambient 

temperature3 

l  .S.  standard 
atmosphere 
temperature 

NO 

\ 

H2<3 

17 

1000 

0 

1.45 

1.30 

1.20 

1.34 

17 

1000 

177.000 

1.37 

1.24 

1.04 

1.18 

17 

333 

177.000 

0.43 

0.40 

0.26 

0.32 

20 

1000 

0 

1.36 

1.17 

1.10 

1.31 

20 

1000 

177.000 

1.29 

1.03 

0.70 

0.91 

20 

333 

177,000 

0.33 

0.27 

-0.02 

0.07 

*The  temperatures  and  pressures  were  computed  at  equilibrium  using  temperature  feedback  and  hydrostatic  adjustment;  they  then 
were  kept  fixed  for  the  perturbation  calculation. 
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magnitude  as  the  difference  between  including  and 
not  including  temperature  feedback. 

2.3  POTENTIAL  C  HANGES  IN 
OZONE  RESULTING  FROM 
OTHER  PERTURBATIONS 

Mans  potentially  significant  perturbations  to 
stratospheric  o/onc  have  been  proposed.  In  order  to 
put  the  effects  of  aircraft  engine  emissions  presented 
in  Section  2.2  in  perspective,  we  present  assessments 


for  other  anthropogenic  perturbations.  Just  as  the 
effects  of  aircraft  emissions  were  show  n  to  depend 
on  the  effects  of  other  perturbing  influences,  such  as 
the  concentration  of  CIX  and  the  atmospheric  tem¬ 
perature  structure,  the  changes  in  O*  resulting  from 
the  other  perturbations  considered  here  also  depend 
on  the  interference  effects  of  multiple  perturbations. 

Our  effort  has  been  primarily  directed  toward 
the  assessment  of  predicted  changes  for  individual 
perturbations  since  future  perturbation  scenarios 
are  often  quite  uncertain.  Of  the  perturbation 
scenarios  considered  here,  only  the  increases  in  the 
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Change  in  temperature  —  K 


FIG.  7.  Change  in  temperature  vs  attitude  for  20-km  injection  perturbation  described  in  Fig.  6  ( 1979a 
chemistry). 


concentrations  of  chlorofluoromethanes  and  C02 
are  well  established.  Increases  in  N20  and  increases 
in  stratospheric  HC3CCI3  are  predicted  on  uncer¬ 
tain  knowledge  of  budgets  and,  in  the  case  of 
CH3CCI3,  uncertain  tropospheric  chemistry. 

Chlorofluoromethanes 

It  has  been  firmly  established  that  the 
chlorofluoromethanes,  CFCI3  and  CF2CI2,  have 
been  increasing  in  the  troposphere  and  stratosphere, 
and  the  observed  increase  is  consistent  with  model 
estimates  based  on  the  historical  production  rates. 
The  CFM’s  are  photolyzed  in  the  stratosphere  to 


yield  free  chlorine  which  may  catalytically  destroy 
ozone. 

Quantitative  estimates  of  the  depletion  of 
ozone  have  changed  as  models,  physics,  and 
chemistry  have  improved.  In  1976  the  predicted 
change  in  total  O3  due  to  steady-state  production  of 
CFM’s  at  1973  rates  was  estimated  to  be  -7.5% 
(National  Research  Council,  1976b).  In  1977  a  ma¬ 
jor  change  in  the  model  chemistry  occurred  when 
the  rate  for 

H02  +  NO  -*■  OH  +  N02  (5) 
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H0C1  h-S  OH  +  C1  .  (7) 


Destruction  of  O,  is  more  efficient  at  lower  altitudes 
if  HOCI  is  included  in  the  model.  When  HOC1  is  in¬ 
cluded,  CIX  shifts  to  more  Cl  and  less  CIO.  Also,  a 
new  catalytic  ozone  destruction  cycle  is  created  that 
does  not  depend  on  the  limited  amount  of  O  atom, 

oh  +  o3  -  ho2  +o2 

Cl  +  O,  -  CIO  +  O, 

(8) 

cio  +  ho2  -  hoci  +  o2 

HOCI  Cl  +  OH 
net  203  -►  3O3 

Figure  9  shows  a  number  of  calculated  time 
histories  for  ozone  under  various  assumptions 
regarding  the  future  release  rates  for  CFM’s.  For 
curve  A  it  was  assumed  that  CFM  production  con¬ 
tinued  at  the  1976  rate  until  1982,  then  it  was 


lowered  by  25%.  In  1987,  it  was  lowered  again  by 
25%.  The  eventual  steady-state  ozone  depletion  for 
this  case  was  7.5%.  For  curve  B,  the  1976  produc¬ 
tion  rate  continued  until  1982  and  was  then  cut  by 
25%.  The  ozone  depletion  at  steady  state  was  10.8% 
for  this  case.  For  curve  C  the  1976  production  rate 
was  used  for  the  entire  time,  and  the  steady-state 
ozone  depletion  was  14.0%.  For  curve  D  the  1976 
production  rate  was  assumed  to  continue  through 
1980,  then  the  CFM  production  was  increased  by 
7%/year  up  to  the  year  2000.  The  steady-state  deple¬ 
tion  for  this  case  is  37.9%.  Case  E  is  almost  the  same 
as  case  D  except  that  the  increased  CFM  produc¬ 
tion  begins  in  1978.  The  rate  of  increase  is  such  that 
the  production  rate  doubles  by  1990  and  doubles 
again  by  2000.  The  steady-state  ozone  depletion  is 
38.6%.  The  1979a  chemical  rates  were  used  in  all  of 
these  model  calculations. 

Increase  in  NzO 

Concern  that  human  perturbations  to  the 
nitrogen  cycle  might  lead  to  enhanced  concentra¬ 
tions  for  atmospheric  N20  has  stimulated  interest  in 


FIG.  9.  Change  in  total  ozone  as  a  function  of  time  for  various  CFM  release  rate  scenerios.  Curve  A:  1976 
CFM  release  rate  until  1982,  reduced  by  25%,  and  further  reduced  by  25%  after  1987,  Curve  B:  1976  release  rate 
until  1982  and  reduced  by  25%  thereafter.  Curve  C:  1976  release  rate  constant  throughout  period.  Curve  D:  1976 
release  rate  through  1980  followed  by  release  rate  increase  of  7%  per  year  up  to  year  2000  followed  by  being  con¬ 
stant  thereafter.  Curve  E:  Same  as  D  except  the  increase  in  release  rate  begins  in  1978.  All  calculations  were 
based  upon  1979a  chemical  rates. 
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The  level  of  atmospheric  N20  determines  the 
atmospheric  response  to  chlorine  changes,  and  vice 
versa,  mainly  because  of  the  coupling  of  chemistry 
by  CIONOt.  Figure  10  illustrates  how  the  percent¬ 
age  change  in  total  O3  with  increased  N20  depends 
on  the  level  of  background  CIX  and  on  the  chemical 
rate  constants  in  the  model.  The  dashed  curves  were 
produced  using  chemical  rate  constants  that  were 
used  in  1978  (see  Luther.  1978).  The  solid  line  refers 
to  results  using  the  1979a  rate  coefficients.  The  ma¬ 
jor  differences  between  the  1978  and  1979 
chemistries  as  they  affect  this  perturbation  are:  (1) 
the  NO  photolysis  rate  is  slower  in  the  1979 
chemistries  so  that  the  background  level  of  NOx  is 
approximately  40%  higher,  and  (2)  the  CIONOi  for¬ 
mation  rate  via  CIO  +  N02^  C10N02  is  almost 
four  times  slower  in  the  1979a  chemistry.  Both 


Multiple  of  present  l\l20 


FIG.  10.  Change  in  total  ozone  due  to  an  increase  in  NjO  expressed  in  multiples  of  present  ground  level 
concentration  which  is  325  ppbv.  Results  produced  using  1979a  rate  coefficients  are  indicated  by  the  solid  line. 
Results  using  chemical  rate  data  that  were  used  in  1978  ( Luther,  1978)  are  shown  by  dashed  curves.  Dash-dot 
line  refers  to  results  using  fast  CIONOj  formation  rate  (see  text). 


the  budget  for  this  gas.  Tropospheric  N20  is  the 
major  source  for  stratospheric  NOx  so  perturba¬ 
tions  in  N20  are  expected  to  alter  stratospheric 
chemistry. 

The  fertilizer  source  of  fixed  nitrogen  is  cur¬ 
rently  estimated  to  be  about  4.2  X  10  tN/yr,  which 
converts  to  a  source  of  1.5  X  10h  tN/yr  as  N20 
(estimate  of  Logan  et  al .,  1978).  Combustion  pro¬ 
vides  about  1 .5  X  106  tN /yr  as  N20  directly  (Weiss 
and  Craig,  1976;  Pierotti  and  Rasmussen,  1976). 
These  anthropogenic  sources  are  thus  currently 
significant  relative  to  the  natural  sources  and  could 
grow  in  the  future  and  lead  to  a  doubling  of  N20 
perhaps  by  the  early  part  of  the  next  century  (Logan 
et  a! .,  1978).  This  doubling  time  is  very  uncertain 
due  to  the  lack  of  detailed  understanding  of  the  cy¬ 
cle  for  N->0. 
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changes  tend  to  diminish  the  predicted  O,  increase 
as  NiO  is  increased.  Even  for  1.95  ppbv  CIX,  ozone 
decreases  for  all  N20  perturbations  with  the  1979a 
chemistry.  This  remains  true  even  for  the  fast  rate  of 
CIONOi  formation  as  indicated  by  the  dot-dash 
line  in  Fig.  10,  so  most  of  the  qualitative  change  in 
results  is  due  to  the  change  in  NO  photolysis. 

Doubling  of  C02 

Systematic  measurements  of  C02  since  1958 
(Keeling  et  al„  1976a  and  b)  have  shown  a  rise  in  at¬ 
mospheric  C02  concentrations  that  has  been  at¬ 
tributed  primarily  to  the  use  of  fossil  fuels.  The  C02 
levels  were  315  ppm  in  1958,  320  ppm  in  1965,  and 
334  ppm  in  1978.  By  the  year  2000  atmospheric  C02 
is  expected  to  be  between  365  and  400  ppm. 
Detailed  prediction  of  the  doubling  time  depends  on 
uncertain  knowledge  of  the  budget  for  C02.  Assum¬ 
ing  that  fossil  fuel  usage  continues  to  increase  at 
4.3%/yr  and  that  about  half  of  the  C02  released 
resides  in  the  atmosphere,  the  atmospheric  C02 
concentration  would  double  by  about  2030,  but  es¬ 
timates  using  other  assumptions  vary  from  2015  to 
2070.  It  is  possible  that  the  concentration  could  be 
limited  to  less  than  500  ppm  by  shifting  away  from 
fossil  fuels  and  relying  more  on  solar  and  nuclear 
energy. 

An  increase  in  CO;  is  expected  to  lead  to 
changes  in  the  thermal  structure  of  the  atmosphere. 
In  particular,  a  doubling  of  C02  has  been  estimated 
to  increase  the  global  mean  surface  temperature  by 
1.5  to  3  K  due  to  the  greenhouse  effect  (Schneider, 
1975;  Augustsson  and  Ramanathan,  1977).  The 
temperature  should  decrease  in  the  stratosphere, 
where  the  infrared  opacity  is  smaller  than  in  the 
troposphere. 

We  tested  the  effect  of  tropospheric  tem¬ 
perature  changes  by  increasing  the  specified  tem¬ 
perature  below  14  km  by  2  K  when  C02  was 
doubled  from  320  to  640  ppm.  The  calculated 
change  in  the  ozone  profile  (Fig.  II)  for  this  case 
was  nearly  the  same  as  for  the  case  where  the 
tropospheric  temperature  remained  unchangeo 
(local  ozone  concentrations  were  within  2.6%  at  all 
altitudes).  The  calculated  temperature  charges 
above  14  km  were  also  similar,  differing  by  less  than 
0.7  K.  In  this  test  the  background  atmosphere 
remained  fixed  (no  hydrostatic  adjustment).  When 
hydrostatic  adjustment  was  included,  the  local 
ozone  increase  near  40  km  was  about  15%',  which 
was  larger  than  the  perturbation  for  no  tropo¬ 


spheric  temperature  change  but  smaller  than  the 
calculated  change  when  hydrostatic  adjustment  was 
neglected  entirely.  Without  hydrostatic  adjustment, 
the  stratospheric  temperature  decrease  leads  to  an 
increase  in  03  at  all  altitudes.  The  O3  increase  tends 
to  offset  the  calculated  temperature  decrease  due  to 
C02  above  40  km  by  increasing  solar  absorption  by 
O3.  The  calculated  temperature  decrease  is  a  max¬ 
imum  near  42  km. 

With  hydrostatic  adjustment,  however,  ozone 
decreased  above  45  km  even  though  the  temperature 
decreased.  This  is  due  to  the  decrease  in  background 
air  density.  As  a  result  of  the  ozone  decrease,  the 
temperature  decreased  further  above  45  km. 

The  increase  in  total  O3  for  a  doubling  of  CO-> 
was  5.06%  with  hydrostatic  adjustment  and  4.74% 
without.  The  larger  change  for  the  hydrostatic  case 
reflects  the  behavior  of  O3  near  25  km,  as  shown  in 
Fig.  1 1 .  At  this  level,  the  photolysis  rates  for  02  and 
O3  both  decreased  for  doubled  C02  as  a  result  of 
the  larger  optical  depth  (more  O3  above  25  km).  The 
change  in  optical  depth  is  larger  for  the  case  without 
hydrostatic  adjustment,  so  the  photolysis  rate  for  02 
decreases  more  for  this  case,  causing  a  smaller  in¬ 
crease  in  O3.  When  the  surface  temperature  was  in¬ 
creased  2  K,  the  change  in  total  ozone  was  4.17% 
without  hydrostatic  adjustment. 

Increase  in  CH3CCI3 

The  use  and  release  of  methylchloroform, 
CH3CCI3,  which  is  used  as  a  cleaning  agent,  has 
been  increasing  at  a  steady  rate.  Its  presence  in  the 
atmosphere  has  been  observed  since  1974  (Cox  et 
al.,  1976).  and  it  has  been  suggested  that  its  con¬ 
tinued  use  will  lead  to  a  reduction  of  ozone 
(McConnell  and  Schiff,  1978).  The  sinks  for 
CH3CCI3  are  photolysis  and  reaction  with  OH  in 
the  stratosphere, 

CH3CCI3  +  ha  -*•  products  (9) 

CH3CC13  +0H->CH2CC13  +  H20  .  (10) 

Stratospheric  destruction  of  CH3CCI3  leads  to 
release  of  Cl  atoms  which  are  able  to  destroy  ozone. 
Reaction  (10)  is  also  effective  in  the  troposphere,  so 
the  growth  of  stratospheric  CH3CCI3  is  limited.  The 
effectiveness  of  reaction  (10)  for  removing 
tropospheric  CH3CCI3.  however,  is  fairly  uncertain. 
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FIG.  II.  Change  in  ozone  concentration  due  to  a  doubling  of  CO2  concentration  ( 1979a  chemistry). 


The  major  uncertainties  result  from  inadequate 
knowledge  of  the  tropospheric  OH  distribution  and 
the  rate  of  reaction.  The  rate  for  reaction  (10)  is  un¬ 
certain  since  there  are  several  conflicting  measure¬ 
ments  of  its  rate  at  room  temperature  and  of  its  tem¬ 
perature  dependence  (JPL,  1979).  Comparison  of 
the  budget  for  CH3CCI3  with  available  measure¬ 
ments  leads  to  an  estimate  for  the  tropospheric 
lifetime  of  between  8  and  11  years  (Penner  and 
Chang.  1978).  Model  calculated  lifetimes  are  con¬ 
siderably  shorter,  implying  larger  tropospheric 
destruction  rates,  on  average,  than  conform  to  the 
measurements.  For  example,  using  the  historical 


release  rate  data  from  Neely  and  Plonka  (1978)  in 
our  one-dimensional  model,  we  calculate  an  average 
abundance  near  the  surface  for  January  1978  of  47. 3 
pptv  using  the  1979a  chemistry.  The  rate  coefficient 
for  reaction  (10)  was  2.5  X  10~12  exp(-1450/T)  (see 
JPL,  1979).  We  calculate  66.1  pptv  using  the  1979b 
chemistry  and  a  rate  coefficient  of  5.4  X  10~12 
exp(-1820/T)  for  (10).  These  rates  differ  by  about  a 
factor  of  2  at  temperatures  characteristic  of  the  up¬ 
per  troposphere  and  by  about  50%  near  the  surface. 
Measurements  taken  by  R.  Rasmussen  (private 
communication.  1979)  give  an  average  concentra¬ 
tion  of  101  pptv  for  CH3CCI3. 
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With  CH3CCI3,  CFCI3,  and  CF:CI2  releases  in¬ 
cluded  in  our  model,  the  calculated  reduction  in 
total  o/one  in  1978  was  1 .7%  (relative  to  1950)  using 
1979b  chemistry,  whereas  without  CH^CCI^,  O, 
decreased  by  1.3%.  The  total  stratospheric  chlorine 
burden  increased  by  7%.  Figure  12  shows  the  time 


history  for  O3  depletion  at  constant  1976  produc¬ 
tion  rates  for  CFCI3  and  CFiCU.  The  CH3CCI3 
release  rate  at  the  earth’s  surface  is  constant  at 
1.23  X  I07  molecules  cm'2  s'!  beyond  1978.  At 
steady  state  ozone  decreased  15.2%,  whereas 
without  CH3CCI3  it  decreased  14.2%o.  Of  course,  the 


FIG.  12.  Enlargement  of  period  from  1980  to  2000  for  calculations  presented  in  Fig.  9,  Curve  C  is 
same  scenario  as  Curve  C  except  1979b  model  chemistry  is  used  with  exception  that  CH3CCI3  is  added. 


effects  of  CH^CCI}  could  be  much  greater  if  its  use 
were  to  continue  to  increase.  Because  of  significant 
uncertainties  in  the  model  treatment  of  photo¬ 
chemistry  in  the  lower  troposphere,  these  results 
must  be  considered  preliminary. 

2.4.  EFFECTS  OF  SPECULATIV  E 
REACTIONS  AND  MECHANISMS 

There  are  several  reactions  that  have  been 
suggested  to  be  of  possible  importance  in  the 
stratosphere  but  have  not  been  measured  vet.  Also, 
there  are  a  few  reactions  for  which  anomalous 
pressure,  temperature  dependencies  have  been 
suggested  but  not  demonstrated.  We  have  examined 
the  sensitivity  of  model  predictions  to  a  subset  of 
such  hypothetical  reactions.  The  subset  selected  u  as 
based  on  the  following  criteria:  (I)  the  hypothesized 
reaction  or  mechanism  has  at  least  some  support 
from  laboratory  work,  and  (2)  the  hypothesized 
reaction  appears  to  have  the  potential  for  causing  a 
major  change*  in  predicted  effects  of  SST's  or 
CFM’s.  The  reasons  for  selecting  these  two  criteria 
were  to  reduce  the  number  of  possible  cases  to  be 
considered  and  to  focus  attention  on  mechanisms 
that  appear  to  have  a  maximum  potential  for  alter¬ 
ing  model  predictions.  It  must  be  stressed  that  most 
of  the  mechanisms  considered  here  have  no  more 
than  a  weak  basis  in  actual  measurement.  The  pur¬ 
pose  of  this  section  is  to  suggest  experiments  that 
might  deserve  some  priority.  It  is  not  to  suggest  that 
the  hypothesized  mechanisms  or  their  computed  ef¬ 
fects  are  particularly  likely  to  be  true. 

In  this  section  we  w  ill  discuss  the  major  effects 
on  the  computed  ambient  atmosphere  and  model 
sensitivity  of  the  various  mechanisms  considered. 
Although  the  changes  to  the  ambient  species  con¬ 
centration  profiles  may  tend  to  either  improve  or 
worsen  the  agreement  with  observations  for  various 
species,  we  have  not  established  a  criterion  for 
agreement  (or  disagreement)  upon  which  to 
evaluate  the  likelihood  that  the  mechanism  actually 
occurs  or  to  eliminate  it  from  consideration.  By 


•In  this  context  a  major  change  would  consist  of  a  change  in  the 
computed  steady-state  effect  of  chlorofluorvimethancs  at  1975 
levels  by  at  least  50*?  or  a  change  in  sign  or  two-fold  increase  m 
computed  SST  effects  Not  ail  of  the  cases  studied  actually  have 
such  large  effects  and  some  mechanisms  not  studied  might  have 
larger  effects 


combining  the  description  of  the  effects  on  the  am¬ 
bient  profiles  discussed  in  this  section  with  the  com¬ 
parison  with  observations  discussed  in  Section  2.1, 
Ihe  reader  can  form  his  own  opinion  on  this  matter. 

Chlorine  Nitrate  Formation 

In  terms  of  published  evaluations  of 
mechanistic  data,  perhaps  the  most  plausible 
speculation  is  that  roughly  75%  of  the  reaction  be¬ 
tween  CIO  and  NOz  forms  a  relatively  short-lived 
species  isomeric  with  chlorine  nitrate,  and  the  effec¬ 
tive  rate  of  chlorine  nitrate  formation  is  only  about 
one-fourth  of  the  observed  rate  of  reaction  between 
CIO  and  N02.  In  JPL  (1979)  (and  NASA,  1979)  no 
firm  choice  is  made  between  the  assumption  that  all 
of  the  reactions  between  CIO  and  N02  lead  to 
chlorine  nitrate  and  the  assumption  that  only  about 
25%  do,  although  most  modelers  have  adopted  the 
first  assumption.  The  choice  of  treatment  for  this 
reaction  is  one  of  the  major  differences  between  the 
1979a  and  1979b  chemistries.  As  can  be  seen  front 
Table  1 1,  if  the  slower  rate  of  chlorine  nitrate  for¬ 
mation  is  assumed,  the  calculated  effects  of  an  SST 
fleet  become  more  positive  than  in  the  base  case, 
while  the  computed  effects  of  CFM's  become  more 
negative.  In  each  case,  the  fractional  change  is  of  the 
order  of  a  factor  of  1.3.  Clearly,  the  correct  treat¬ 
ment  of  chlorine  nitrate  formation  is  of  substantial 
importance.  It  is  also  worth  noting  that  if  the 
isomeric  products  have  greater  photolytic  stability 
than  CI0N02,  or  if,  as  impurities,  they  account  for 
some  of  the  observed  absorption  spectrum  of 
CIONOs,  the  effect  of  chlorine  nitrate  isomers  on 
predictions  might  actually  be  reversed.  This  could 
happen  if  the  mean  rate  of  loss  of  some  CINOj 
products  were  to  be  slower  than  that  presently 
assumed  for  ClONOj  so  that  larger  concentrations 
of  CINO,  species  would  be  generated. 

When  the  low  rate  is  assumed  for  chlorine 
nitrate  formation,  a  significant  (>10%)*  change  in 
the  ambient  state  is  computed  for  Cl.  CIO.  HCI, 
HOCI,  and  CINOs,  all  of  w  hich  increase  by  nearly  a 
factor  of  2  at  25-30  km,  while  CINO^  is  reduced  by 
a  factor  of  3  to  4  at  almost  all  altitudes.  Thus,  this 
hypothesis  exacerbates  the  apparent  disagreement 
between  observation  and  calculation  for  the  shape 
of  the  CIO  vertical  profile  and  destroys  the  apparent 


•All  percent  changes  reported  in  this  section  refer  to  fmaiimal 
changes  in  the  quantity  relative  to  the  standard  (1979b  model) 
ambient  conditions  or  model  sensitivity . 
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TABLE  1 1 .  Results  of  sensitivity  calculations  for  speculative  chemical  reactions  and  mechanisms. 


Model  content 

Ambient  ozone 

column 

(Dobson  units) 

Change  in  total  ozone  (%) 

17-km  NO*  20-km  NO  * 

injection3  injection3 

CFM 

production** 

Baseline  Model  ( 1979b)-NASA  1049  (1979)  2  ppb  OX 

322 

1.25 

1.14 

-14.6 

Slower  rate  of  CIO  NO  2  formation 

312 

1.42 

1.54 

-19.3 

Photolysis  of  XONC^  to  XO  +  NOj  products 

324 

1.35 

1.20 

-12.2 

Pressure-  and  temperature-dependent  HO*  disproportion¬ 
ation  reactions 

339 

0.47 

-0.40 

-8.9 

HC1  production  from  CIO  +  OH,  CIO  +  HO2 

334 

0.74 

0.02 

-4.4 

CIO 3  production,  K  j|ibrjum  as  for  CIOj  from  NASA 
1010;  no  other  reaction  of  CK)^ 

322 

1.24 

1.11 

-14.0 

OClOj  production  ^equilibrium  as  for  ClOj  from  NASA 
1010;  subsequent  chemistry  described  in  text 

337 

1.51 

1.40 

-6.91 

0CI°2  production  Kequilibrium  as  for  C102  from  JPL 
(1979);  subsequent  chemistry  described  in  text 

323 

1.27 

1.16 

-14.0 

OCI02  production  ^equilibrium  as  for  C102  from  NASA 
1010;  subsequent  chemistry  described  in  text  and 
P,T-dependent  HO*  disproportionation  reactions 

353 

0.49 

-0.51 

0.73 

aRate  of  emission  is  1000  molecules  (NO)  cm~3s-l  OVei  a  1-km-thick  layer. 
^Constant  release  rate  at  1976  levels. 


agreement  between  calculation  and  the  preliminary 
chlorine  nitrate  measurement  reported  by  Murcray 
(1979). 

Photolysis  of  X0IN02  Species 

Almost  equally  uncertain  are  the  products  of 
the  photolysis  of  species  of  the  form  XONO,,  where 
X  =  OH,  Cl,  NO,.  In  the  base  case,  it  has  been 
assumed  that  all  of  these  reactions  lead  to  X  +  N03. 
This  assumption  is  largely  based  on  the  experimen¬ 
tal  results  of  Chang  el  al.  (1979)  for  chlorine  nitrate. 

Other  products  are  possible,  perhaps  the  most 
distinct  (in  terms  of  model  predictions)  would  be 
XO  +  NO,  (or  for  N203,  2NO,  +  O).  In  a  com¬ 
putation  in  which  all  of  these  photolyses  were 
assumed  to  yield  XO  +  HO,  products  (except  N,Os 
which  was  assumed  to  yield  2NO,  +  O),  the  com¬ 
puted  effects  of  SST  emissions  became  more 
positive  by  5-10%  (of  the  change  computed  for  the 
base  case)  while  the  effects  of  CFM's  became  less 
negative  by  about  1 7%.  Thus,  resolving  the  question 


of  photolysis  products  of  XOCIO  species  is  of 
modest  importance  to  perturbation  calculations. 

In  terms  of  species  concentrations,  only  NO3 
and  N,03  displayed  large  changes  (both  were 
reduced  by  a  factor  of  two  to  three  at  the  region  of 
their  largest  mixing  ratios,  and  by  up  to  90%  at 
some  (lower)  altitudes).  Several  other  species  dis¬ 
played  changes  on  the  order  of  10-15%  at  some 
altitudes. 

Pressure-Dependent  Rates  for  HO* 
Disproportionation  Reactions 

A  significant  pressure,  inverse  temperature  and 
water  vapor  dependence  has  been  reported  for  the 
reaction  HO,  +  HO,  MH-°  H,Q,  +  0,(Cox.  1978; 
Hamilton  and  Lii,  1977;  Cox  and  Burrows,  1979). 
Furthermore,  while  the  measurements  of  the  reac¬ 
tion  HO,  +  OH  at  low  pressures  seem  only  mildly 
inconsistent  with  each  other,  the  indirect  values  in¬ 
ferred  at  higher  pressures  (Hochanadel  et  al.,  1972; 
DeMore  and  Tschuikov-Roux,  1974;  DeMore, 


1979)  are  substantially  (two  to  four  fold)  faster  than 
the  values  measured  at  low  pressures.  In  this  con¬ 
text,  the  poorly  understood  pressure,  temperature, 
and  water  vapor  dependence  of  H02  +  H02 
suggests  the  possibility  that  H02  +  OH  may  also  be 
affected  by  similar  phenomena.  As  can  be  seen  in 
Table  11,  when  we  assumed  that  the  expression 
given  by  Cox  (1978)  for  kH02+M02, 

kH02  +H02  = 

3.25X  108+4X  nr1  °(M)/ [  1  +3. SX 1  Q-‘ 6(M)e~2Q60/T  ] 
8(M)+4.08X1020 

(II) 

and  an  expression  for  koH  +  ttO,  designed  to  yield  a 
low-pressure  value  compatible  with  Chang  and 
Kaufman  (1978),  and  a  high  pressure  value  com¬ 
patible  with  the  larger  rate  constants  reported  at 
high  pressure,  i.e., 

kH02  +  OH  = 

(5  X  10'30  X  M  +  3X  10'1 1 K300/T)3  ,  (12) 

there  were  major  changes  in  model  predictions  for 
both  SST’s  and  CFM’s.  Here  it  should  be 
emphasized  that  the  expression  used  for  the  rate 
constant  for  reaction  of  H02  with  OH  is  completely 
arbitrary.  It  was  designed  to  approximately  fit  high- 
and  low-pressure  data.  Had  a  Lindeman- 
Hinshelwood  expression  been  fitted  to  the  same 
data,  or  had  a  constant  ratio  to  the  H02  +  H02  rate 
been  assumed,  the  effects  would  have  been  different 
(probably  larger)  but  equally  plausible.  The  form 
chosen  reflects  the  postulation  of  a  pressure  in¬ 
dependent  abstraction  mechanism  coupled  with  a 
pressure  dependent  enhancement  reflecting  complex 
formation  (with  a  limiting  value  well  above 
1  X  10'10  at  room  temperature). 

Elucidation  of  the  pressure  and  temperature 
dependencies  of  the  HOx  radical  disproportionation 
reactions  should  be  of  high  priority  to  those  con¬ 
cerned  with  SST  effects  and  of  considerable  impor¬ 
tance  to  those  concerned  with  CFM  production  ef¬ 
fects.  The  assumed  change  in  the  H02  +  OH  reac¬ 
tion  is  substantially  more  important  than  the 
change  in  the  H02  +  H02  reaction  within  the 
stratosphere.  However,  the  H02  +  H02  reaction  is 
of  significant  importance  in  the  troposphere,  and 
the  effects  of  subsonic  aircraft  operations  are  ex¬ 
pected  to  be  sensitive  to  it,  as  would  be  the  effects  of 
CH3CCI,  emissions. 


When  the  above  assumptions  were  made  about 
the  HOx  disproportionation  reactions,  many  species 
changed  by  roughly  a  factor  of  2  between  approx¬ 
imately  10  and  20  km.  Roughly  two-fold  reductions 
(between  approximately  10  and  20  km)  occurred  for 
OH,  H02,  Cl,  CIO,  and  CH20,  while  H202  was 
reduced  by  roughly  the  square  of  the  fractional 
reduction  in  HO-i.  Computed  NO  and  N02  dis¬ 
played  two-fold  increases  over  the  same  altitude 
range.  Species  that  were  only  slightly  changed  in¬ 
cluded  HNO3  (<10%  changes)  and  CH4  (approx¬ 
imately  10%  increase  at  high  altitudes). 

Thus,  changes  in  the  HOx  disproportionation 
rates  akin  to  those  considered  here  would  cause  sub¬ 
stantial  changes  in  computed  profiles  of  short  lived 
species  in  the  region  from  10  to  20  km.  Some  of 
these  changes  tend  to  improve  agreement  with  ob¬ 
servation,  others  tend  to  reduce  it.  Although  CH4 
increases  at  higher  altitudes,  the  increase  is  not  large 
enough  to  demand  a  major  revision  in  K,  even  if 
future  experiments  should  suggest  the  pressure- 
dependent  disproportionation  rates  used  here  to  be 
fair  approximations  of  reality.  For  CH4  to  increase 
substantially  at  higher  altitudes,  it  would  be 
necessary  for  an  increase  in  k0H+HOi  10  extend  to 
substantially  lower  air  densities. 

All  of  the  above  speculations  have  at  least  some 
direct  basis  in  laboratory  measurements.  The  specu¬ 
lations  that  follow  are  even  more  tenuous  than  the 
above. 


HCI  Formation  from  HOx  +  CIO 

The  reaction  between  OH  and  CIO  could  yield 
H02  +  Cl  or  HCI  +  02  at  low  pressures,  or  a  three- 
body  process  might  yield  a  moderately  stable 
species  of  the  form  HOOCI  (a  peroxide)  or  HOCIO 
(an  acid).  Either  of  these  last  species  might  be  ex¬ 
pected  to  be  unstable  with  respect  to  dispropor¬ 
tionation  in  a  condensed  phase,  but  gas  phase 
stability  seems  possible.  In  any  case,  the  production 
of  HCI  +  02  from  HO  +  CIO  would  seem  likely  to 
have  the  most  drastic  effect  on  model  predictions. 
The  overall  reaction  has  a  rate  constant  of 
9,2  X  I0-11  at  room  temperature,  and  the  products 
H02  +  Cl  account  for  at  least  65%  of  the  products 
(l.eu  and  Lin,  1979).  These  products  have  little  ef¬ 
fect  on  model  sensitivity. 

An  investigation  of  the  effect  of  this  reaction 
producing  HCI  assuming  a  rate  constant  of 
2  X  IfH2  cm  \  s  suggested  that  model  sensitivity  to 
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a  Cl  M  perturbation  could  be  reduced  by  approx¬ 
imately  40%;  at  I  X  Kb1-  cm-'/'s  the  sensitivity 
would  be  reduced  by  roughly  20%. 

A  similar  argument  could  be  applied  to  the 
reaction: 

HO,  +CIO  -  HOCI  +  O,  (13) 

HO,  +  CIO  ‘  HC1  +03  .  (14) 

The  only  evidence  suggesting  a  plausible  role  for 
reaction  (14)  is  that  at  temperatures  below  room 
temperature,  the  rate  of  the  overall  reaction  in¬ 
creases  substantially  as  the  temperature  is  reduced 
(C  Howard  et  at.,  NOAA  Aeronomy  Laboratory, 
Boulder,  private  communication,  1979).  This 
suggests  the  possibility  that  a  five-centered  complex 
might  be  formed  and  HC1  +  Oy  e  -ninated.  If  the 
rate  of  HC1  formation  via  this  reaction  were  to  be  as 
large  as  1  X  Kb1-  cm-’/s.  then  the  effect  of  CFM’s 
would  be  reduced  by  approximately  50%,  and,  if  the 
rate  were  to  be  even  larger,  increases  might  be  com¬ 
puted  from  CFM  releases. 

Table  1 1  gives  the  sensitivities  obtained  when 
both  HO  +  CIO  -  HCI  +  Oy  and  HOy  +  CIO  - 
HCI  +  Oy  were  assumed  to  have  rate  constants  of 
I  X  Kb1-  ent-’/s.  As  is  evident,  these  reactions 
might  have  a  major  impact  on  both  CFM  and  SST 
effects,  and  they  might  have  a  smaller  but  signifi¬ 
cant  effect  on  computed  CFM  effects  at  rates  as  low 
as  approximately  I  X  l(b|1cm’  s. 

The  effect  on  the  model  sensitivity  to  an  NO, 
perturbation  is  largely  a  result  of  there  being  a 
lessened  rate  of  CIX  and  HO,  catalytic  destruction 
of  ozone  with  which  NO,  can  interfere.  Also,  null 
cycle  sequences  involving  CIX  and  HO,  are  not  as 
effective  in  competing  with  ozone-destroying 
sequences  for  NO.  The  reduced  sensitivity  to  CFM 
perturbations  results  from  both  an  increase  in  the 
rate  at  w  hich  HCI  is  formed  and  a  transfer  of  CIO 
and  HO,  radicals  out  of  odd-oxygen  destroying 
sequences  and  into  null  sequences. 

There  is  the  possibility  that  the  reaction  HOy  + 
HOy  -  HyO  +  Oy  might  occur  to  some  extent.  An 
analysis  of  available  data  (H.  Johnston,  University 
of  California  at  Berkeley,  private  communication, 
1979)  suggests  an  upper  limit  of  about  5%  of  the 
total  rate  for  this  reaction.  If  the  similar  reaction 
HOy  +  CIO  -•  HCI  +  Oy  had  a  branching  ratio  of 
less  than  0.05.  it  would  have  no  more  than  a  modest 
effect  on  model  sensitivity  (based  on  calculations 


made  in  1978,  the  reaction  HOy  +  HOy-.  HyO  + 
Or  would  need  to  have  a  branching  ratio  of  nearly 
0.1  to  be  of  even  marginal  significance  in  strato¬ 
spheric  perturbation  studies). 

The  model-calculated  ambient  atmosphere 
with  rates  for  both  OH  +  CIO  -*  HCI  +  Oy  and 
HOy  +  CIO-.  HCI  +  03  set  to  Kb1-  cm’  s  con¬ 
tained  about  half  as  much  Cl,  CIO,  ClONOy, 
CINOy  and  HOCI  as  did  the  normal  ambient, 
whereas  HCI  increases  of  about  a  factor  of  2  were 
calculated  near  30  km.  No  other  measured  species 
were  significantly  (>10%)  affected,  although  H  Oy 
and  CH,0 OH  were  reduced  by  about  30%  between 
25  and  35  km. 

Although  model  predictions  are  indeed  sen¬ 
sitive  to  these  reactions,  the  likelihood  that  they  oc¬ 
cur  is  probably  not  very  great.  Thus,  while  an  effort 
should  be  made  to  measure  them,  an  upper  limit  less 
than  approximately  10-1’  cm’  s  would  resolve  most 
of  the  issues  they  raise. 

CIO,  Chemistry 

A  final  mechanism  that  was  studied  was  the 
possibility  that  Oy  and  CIO  add  to  form  a  molecule 
with  a  binding  strength  comparable  to  that  of 
CIOO,  the  complex  formed  from  Cl  and  Oy.  There  is 
evidence  for  some  such  phenomenon  from  the  effect 
of  added  Oy  on  the  quantum  yield  for  ozone  loss  in 
systems  containing  Cly,  Oy  and  Oy  (Wongdontri- 
Stuperer  «/.,  1978;  J.  Birks.  University  of  Colorado, 
private  communication.  1979). 

If  it  is  assumed  that  the  formation  and  decom¬ 
position  rate  constants  for  the  process  CIO  + 
Oy  ->  ClOy  arc  those  given  in  NASA  1010  (1977)  for 
Cl  4-  Oy  ->  ClOy  and  no  subsequent  chemistry  oc¬ 
curs  (only  the  formation  and  thermal  decomposi¬ 
tion  of  the  complex  occur),  then  model  sensitivities 
are  virtually  unaffected.  Approximately  2  to  5%  of 
the  inorganic  chlorine  in  the  lower  stratosphere  is  in 
the  form  of  CIO,,  and  all  computed  changes  in  am¬ 
bient  concentrations  are  small.  Even  though  the 
coupling  between  CIO  and  ClOy  is  very  rapid,  and 
ClOy  is  computed  to  be  larger  than  CIO  between  ap¬ 
proximately  6  and  22  km.  most  of  the  ClOy  ap¬ 
parently  comes  from  the  much  more  abundant  HCI 
and  ClONOy.  However,  if  the  structure  of  the  ClOy 
complex  were  OClOy.  then  it  is  possible  that  the 
following  reactions  would  occur: 

OCI02  O  +  OCIO  (15) 


* 


f 

> 


and 

0C102  +N0  -  0C10  +  N02  •  (16) 

Because  CIO  and  OCIO  have  comparable  heals 
of  formation,  the  endoergicity  of  the  photolysis 
would  approximate  60  kcal  plus  the  binding  energy 
of  the  complex.  This  suggests  an  endoergicity  for 
reaction  (15)  comparable  to  that  for  NO_i 
photolysis.  Thus,  at  least  in  principle,  reaction  (15) 
might  have  an  effective  J  value  as  large  as  about 
10--V.  The  reaction  OCIO  +  NO  -  NO,  +  CIO 
has  a  rate  constant  of  about  3  X  I0-14  cm'/s  at 
room  temperature  (JPL,  1979).  Given  the 
similarities  in  the  reactions,  a  similar  sort  of  rate 
constant  might  be  estimated  for  0C102  +  NO.  If 
one  assumes  that  the  J  value  for  0CI02  photolysis  is 
I0"4s_l.  that  the  rate  of  reaction  with  NO  is 
IX  1(H~  cm-1  s.  and  that  the  formation  and 
decomposition  rates  are  those  given  for  CIO-,  in 
NASA  1010  (1977),  then  one  estimates  a  moderate 
(20%)  increase  in  the  (small)  ozone  increase  com¬ 
puted  from  SST  operations  and  more  than  a  50% 
reduction  in  the  effect  of  CFMs.  If  on  the  other 
hand  one  assumes  the  same  subsequent  chemistry 
but  a  lower  stability  for  the  complex,  say  that  given 
for  CI02  in  JPL  (1979).  then  no  significant  effect  on 
model  sensitivity  is  predicted. 

Thus,  for  the  formation  of  a  complex  between 
02  and  CIO  to  significantly  affect  model  sensitivity, 
the  complex  must  have  a  binding  energy  ap¬ 
proaching  6  kcal.  and  either  photolysis  or  reaction 
with  NO  must  be  reasonably  fast.  However,  if  the 
complex  were  to  be  stable  (bound  by  more  than  8 
kcal)  and  the  subsequent  chemistry  were  reasonably 
fast,  one  might  even  compute  a  decrease  in  ozone 
for  NO,  injections  and  an  increase  in  column  ozone 
from  CFM  increases.  Similarly,  if  the  HO,  reactions 
were  pressure-dependent  and  0CI02  were  to  have 
the  chemistry  discussed  above,  one  would  calculate 
a  decrease  in  column  ozone  for  the  20-km  NO,  in¬ 
jection  and  an  increase  in  column  ozone  from  the 
CFM  scenario  (see  Table  1 1).  There  is,  of  course,  a 
lengthy  chain  of  speculative  assumptions  required 
to  achieve  such  results. 

Even  though  the  OC102  chemistry  discussed 
above  is  completely  speculative,  its  potential  for 
large  effects  may  justify  the  attempt  to  study  it  by 
those  chiefly  interested  in  CFM  effects. 

Finally,  it  should  be  noted  that  the  simplest 
(but  probably  not  the  most  likely)  route  back  to 


computed  effects  of  SST  operations  (like  those  ob¬ 
tained  in  C1AP  Monograph  3  (1975)  and  National 
Research  Council  (1975a)  would  be  for  the  HO,  dis¬ 
proportionation  reactions  to  have  fast  rates  (either 
via  the  pressure/temperature  effects  hypothesized 
above,  or  through  an  error  in  the  post  1976 
measurements)  and  for  the  reaction  of  H02  +  NO 
to  be  slower  than  is  now  thought  to  be  the  case. 
Although  it  seems  unlikely  that  the  grow  ing  body  of 
experimental  results  involving  H02  chemistry 
would  all  be  wrong  (at  least  at  low  pressures  where 
many  of  the  measurements  have  been  direct  and  the 
data  analysis  reasonably  straightforward),  it  cannot 
be  denied  that  the  new  measurements  have 
produced  a  startling  number  of  unexpected  negative 
temperature  dependencies  for  apparently 
bimolecular  reactions  involving  H02  (and  CIO  as 
well).  While  theories  of  reaction  rates  can  no  doubt 
be  created  to  fit  such  data,  the  pre-existing  theories 
do  not  easily  lead  to  the  observed  rate  constants. 
Although  it  seems  very  likely  that  the  resolution  of 
such  problems  will  indeed  come  via  a  modification 
of  theory,  additional  confirmatory  measurements  of 
some  of  the  anomalous  temperature  dependences 
using  independent,  if  perhaps  less  direct,  techniques 
would  still  be  welcomed. 


2.5  STRATOSPHERIC  WATER 
VAPOR 

The  review  of  our  knowledge  and  obser¬ 
vational  data  on  stratospheric  H20  revealed  several 
apparent  paradoxes  which  remain  to  be  resolved. 
These  are  illustrated  by  the  data  summarized  in 
Figs.  13  and  14.  These  data  indicate  upward  and 
poleward  directed  gradients  of  H20  immediately 
downstream  (in  the  Hadley  flow)  from  the  tropical 
tropopause  source  for  air  (and  contained  tracers  in¬ 
cluding  H20)  entering  the  stratosphere.  As  pointed 
out  previously  (Ellsaesser,  1974),  these  appear  to  re¬ 
quire  a  sink  for  H20  in  the  stratosphere.  However, 
possible  sinks  proposed,  such  as  the  Antarctic  win¬ 
ter  freeze-out  and  the  “cold  finger”  of  overshooting 
cumulonimbi,  do  not  appear  adequate.  An  alter¬ 
native  possibility  is  that  observational  data  have 
been  accumulated  mostly  from  the  downwind  edges 
of  continents  where  cumulonimbi  penetrations  are 
most  likely  to  be  concentrated;  thus  the  gradients 
indicated  in  Figs.  13  and  14  may  exist  only  in  such 
areas  and  not  be  representative  of  zonal  averages  as 
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FIG.  13.  Summary  of  observations  of  the  latitudinal  gradient  in  stratospheric  H20.  Data  sources  identified  by 
references. 


they  have  been  interpreted  to  be.  it  valid,  disap¬ 
pearance  of  these  gradients  would  obviate  the  need 
for  an  HiO  sink  in  the  stratosphere. 

An  even  greater  paradox  is  indicated  by  two 
other  observational  phenomena:  ( 1 )  the  frequent  oc¬ 
currence  of  dry  air  with  mixing  ratios  <3  ppmm 
(i.e.,  typical  of  the  lower  stratosphere)  as  much  as 
100  mb  beneath  the  tropical  tropopause  in  clearly 
tropospheric  air  (such  mixing  ratios  are  not  only  un¬ 
typical  of  what  has  until  now  been  our  concept  of 
tropospheric  air,  but  they  are  actually  lower  than 
would  be  produced  by  passage  through  the  overly¬ 
ing  ambient  tropical  tropopause)  and  (2)  on  those 


so-called  “wet"  days  when  soundings  do  conform  to 
the  Brewer-Dobson  theory  by  showing  saturated  air 
at  and  below  the  tropical  tropopause  and/or  when 
cumulonimbus  penetrations  of  the  tropopause  are 
sighted,  the  mixings  ratios  at  and  just  about  the 
tropical  tropopause  are  generally  >6  ppmm. 

These  observations  appear  to  imply  that  the 
aridity  of  the  stratosphere  is  not  controlled  by 
freeze-drying  at  the  tropical  tropopause  and  that 
there  exists  in  the  troposphere  a  mechanism  for  dry¬ 
ing  the  tropical  tropospheric  air  above  200  mb  to 
mixing  ratios  <3  ppmm  typical  of  the  lower 
stratosphere.  As  of  now  the  only  mechanism  which 
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suggests  itself  is  cirrus  clouds.  If  there  exist  seeding 
agents,  possible  including  the  cirrus  ice  crystals 
themselves,  which  at  temperatures  of  -50  to  -80°C 
and  pressures  below  200  mb.  can  condense  vapor 
out  of  what  is  now  regarded  as  unsaturated  air  and 
cause  the  resulting  hydrometeors  to  fall  to  higher 
pressures  and  temperatures  before  re-evaporating, 
we  could  explain  many  currently  puzzling  observa¬ 
tions.  An  additional  strong  argument  for  a 
tropospheric  mechanism  for  drying  the  air  to 
stratospheric  mixing  ratios  is  that  it  would  open  a 
new  pathway  for  significant  transfers  of  air  from  the 
troposphere  to  the  stratosphere— and  that  is 
through  the  tropopause-gap.  Many  have  suggested 
this  pathway  before  based  on  tracers  such  as  0<,  but 
most  have  considered  it  as  insignificant  since  it  ap¬ 
peared  to  offer  no  way  of  maintaining  the  aridity  of 
the  stratosphere.  However,  if  the  air  is  dried  bv  a 
mechanism  operating  in  the  troposphere,  this  objec¬ 
tion  goes  away. 

A  third  paradox  is  due  to  the  H:0  believed  to 
be  added  to  the  stratosphere  by  the  oxidation  there 
of  CH4  entering  with  the  Hadley  circulation  from 
the  troposphere.  This  should  impose  a  downward 
directed  gradient  of  HsO  in  the  middle  and  upper 
stratosphere  with  a  total  difference  in  HA)  mixing 
ratio  somewhat  less  than  2  ppmm  and  an 
equatorward  directed  gradient  in  the  lower 
stratosphere  since  the  return  Hadley  flow  through 
the  polar  tropopauses  presumably  must  carry  more 
HiO  than  is  carried  by  the  entering  Hadley  flux 
through  the  tropical  tropopause.  A  downward 
directed  gradient  is  shown  by  many  soundings,  but 
most  of  these  show  a  total  difference  in  mixing  ratio 
of  more  than  2  ppmm.  Also,  an  equatorward  di¬ 
rected  gradient  of  approximately  ihe  righi 
magnitude  is  shown  by  what  are  considered  to  be 
the  two  most  reliable  seis  of  observations  of 
stratospheric  HA),  the  MRI*  and  the  Mastcn- 
brook  series.  Such  a  gradient  is  ,dso  shown  by  the 
observations  of  McKinnon  and  Moiewood  (|R~0) 


•MRI  is  (he  licronwri  lor  the  British  meir.'f.*i« .i1 
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(see  fig.  I')  However,  for  the  inferences  drawn 
above  to  he  correct  requires  that  the  data  of  most 
investigators  of  stratospheric  HA)  are  incorrect  or 
unrepresentative. 
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3.  SATELLITE  OZONE  DATA  PROCESSING 
ARC  HIVING,  AND  ANALYSIS 


3.1  OVERVIEW 

Measurements  of  o/one  have  been  made  on  a 
regular  basis  for  several  decades.  The  Dobson  spec¬ 
trophotometer  instrument  is  the  most  widely  used 
ground-based  instrument.  This  instrument 
measures  the  differential  attenuation  of  sunlight  in 
adjacent  spectral  bands  in  the  UV  Huggins  hand  of 
ozone  from  which  the  total  column  of  o/one  is 
determined.  The  Dobson  instrument  has  been  in  use 
since  1929,  and  approximate!)  100  o/one  obser¬ 
vatories  have  used  this  instrument.  Since  I960  the 
data  from  many  of  the  Dobson  Mtes  have  been  com¬ 
piled  and  published  by  the  Xtmospheric  Environ¬ 
ment  Service  of  Canada.  Presently,  data  from  61 
Dobson  sites  are  being  published  routinely.  These 
sites  are  located  principally  on  continents  in  the 
Northern  Hemisphere.  The  largest  data  voids  are 
over  the  oceans  and  in  the  Southern  Hemisphere. 

Global  satellite  measurements  of  total  o/one 
were  first  made  in  1969  using  the  infrared  inter¬ 
ferometer  spectrometer  (IRIS).  1  ater.  the  backscat- 
ter  ultraviolet  spectrometer  (Bl'V)  and  the  IRIS 
sensors  were  both  flow  n  on  the  Nimbus  3  (launched 
in  1970)  and  Nimbus  4  (launched  in  1972)  satellites. 
The  longest  continuous  satellite  o/one  data  record 
is  that  of  the  Nimbus  4  BL’V  instrument  that 
operated  from  April  1972  to  July  1977.  The  mul¬ 
tichannel  filter  radiometer  (MI  R)  infrared  sensor 
was  first  flown  in  1977  on  a  Defense  Meteorological 
Satellite  Program  (DMSP)  Block  5D  series  satellite, 
which  is  operated  by  the  I  S.  An  Force.  Four 


satellites  in  this  series  have  been  launched,  and  the 
data  from  these  satellites  are  expected  to  extend  into 
1980.  The  period  covered  by  each  of  the  four  DMSP 
satellites  is  shown  in  Table  12. 

In  addition  to  the  DMSP  satellite,  o/one  data 
are  being  obtained  by  sensors  aboard  the  NASA 
Nimbus  7  and  NOAA  TIROS  N  satellites  that  were 
launched  in  November  1978.  The  TIROS  N  satellite 
carries  a  high  resolution  infrared  sensor  and  the 
Nimbus  7  satellite  carries  a  solar  BUV  sensor  and  a 
TOMS*  instrument.  The  series  of  NASA  and 
NOAA  satellites  is  expected  to  prov  ide  global  ozone 
data  into  the  late  1980‘s,  The  DMSP  MFR  sensors 
were  the  only  satellite  o/one  monitoring  systems 
operating  between  mid- 1977  and  November  1978. 

All  of  the  DMSP  satellites  are  in  polar,  sun- 
synchronous  orbits  with  basically  the  same  orbital 
parameters.  F4  is  somewhat  different  in  that  it  is  in 
a  nighttime  ascending  orbit  (i.e..  F4  ascends  from 
southeast  to  northwest  on  the  night  side),  whereas 
the  others  are  in  daytime  ascending  orbits.  Conse¬ 
quently.  although  F2  and  F4  are  only  about  an  hour 
apart  in  their  overpass  times,  their  flight-track  (and 
site-scan)  geometries  are  different,  thereby  better 
enabling  F4  to  fill  data-void  areas  of  F2  and  F3.  We 
receive  the  MFR  radiance  data  on  magnetic  tape 
from  the  U.S.  Air  Force  Global  Weather  Central, 
and  the  raw  data  are  forwarded  to  the  Environmen¬ 
tal  Data  Service  of  NOAA  where  they  are  archived. 


•Total  O/one  Mapping  Spectrometer 


TABLE  12.  Data  periods  covered  by  MFR  sensors3  on  DMSP  satellites. 


DMSP  satellite 

Data  period 

Local  overpass  time 

FI 

March  1977  -  July  1977 

11*12  a.m.  and  p.m. 

F2 

July  1977  -  Spring  1980  (estimated)*5 

9-10  a.m.  and  p.m. 

F3 

August  1978  -  January  1980 

6-7  a.m.  and  p.m. 

F4 

June  1979  -  January  1980 

10-11  a.m.  and  p.m. 

*Thc  F4  MFR  sensor  may  be  turned  hack  on  at  a  later  time  if  satellite  battery  problems  are  solved  or  if  there  is  a  change  in  the 
power  management  plan. 

^F2  has  already  exceeded  its  designed  lifetime. 


The  MFR  sensor  measures  infrared  radiances 
for  16  channels.  The  channel  characteristics  are 
given  in  Table  13.  Total  ozone  amounts  are  deter¬ 
mined  from  sets  of  radiance  measurements  using  an 
empirical  relationship  that  is  developed  using  linear 
regression  analysis.  Total  ozone  is  currently 
modeled  as  a  linear  combination  of  terms  involving 
functions  of  the  MFR  radiances  for  four  channels 
(1,  2,  7,  and  16)  and  the  secant  of  the  zenith  angle. 
The  coefficients  used  in  the  retrieval  model  are 
determined  by  regression  analysis  using  sets  of 
simulated  radiances  derived  from  detailed  radiative 
transfer  calculations.  Historical  vertical  tem¬ 
perature  and  ozone  concentration  profiles  are  used 
as  input  to  the  calculations,  and  calculations  are 
made  for  various  cloud  amounts,  surface  tem¬ 
peratures,  and  instrument  scan  angles.  Different 
sets  of  regression  coefficients  are  derived  for  each 
latitude  band  (currently,  we  are  using  1 1  bands)  and 
for  different  times  of  the  year. 

Temperature-dependent  molecular  absorption 
coefficients  for  the  9.6-jum  ozone  band  and  the  15- 
Mtn  carbon  dioxide  band  are  determined  from  high- 
resolution,  line-by-line  calculations.  These  coef¬ 
ficients  are  used  to  compute  transmittance  profiles 
for  the  various  channels,  given  the  temperature  and 
optically  active  gas  concentration  profiles. 
Transmittances  due  to  the  water  vapor  continuum 
absorption  in  the  atmospheric  window  are 
calculated  with  an  “e-type”  model.  A  regression 
model  based  on  data  from  line-by-line  calculations 
plus  continuum  effects  is  used  for  the  tropospheric 
water  vapor  transmittance  calculations  in  the  l5-/um 
spectral  region.  A  Goody  random  band  model  is 
employed  to  calculate  the  transmittances  due  to  the 
14-#im  ozone  band.  See  Lovill  et  at.  (1978)  for  more 
detail  regarding  the  retrieval  methodology  as  it  was 
employed  during  the  feasibility  study. 

The  total  column  ozone  data  derived  from  the 
MFR  measurements  are  compared  with  an  indepen¬ 
dent  set  of  ground-based  measurements  of  total 
ozone  in  order  to  evaluate  the  MFR-retrieval 
methodology.  This  set  of  data  consists  of  Dobson 
ozone  measurements  that  are  taken  close  in  time  to 
DMSP  satellite  overpasses.  Currently,  36  Dobson 
observatories  are  providing  special  sets  of  ozone 
measurements  for  comparison  purposes.  We 
provide  each  observatory  predictions  of  the  daily 
overpass  times  for  each  DMSP  satellite.*  The  North 


TABLE  1 3.  Nominal  MFR  channel  characteristics.* 


Channel 

number 

Center 

Half 

width 

(cm-1) 

Species 

NESRb 

(uml 

(cm  *) 

i 

15.0 

668.5 

3.5 

co2 

0.30 

2 

14.8 

676.0 

10.0 

C02 

0.09 

3 

14.4 

695.0 

10.0 

cxi2 

0.10 

4 

14.1 

708.0 

10.0 

C02 

0.11 

5 

13.8 

725.0 

10.0 

C°2 

0.11 

6 

13.4 

747.0 

10.0 

co2 

0.12 

7 

12.0 

835.0 

8.0 

Window 

0.11 

8C 

18.7d 

535.0 

16.0 

h2o 

0.15 

9 

24.5d 

4083 

12.0 

h2o 

0.14 

10 

22. 7d 

441.5 

18.0 

h2o 

0.09 

11 

23.9d 

420.0 

20.0 

h2o 

0.12 

12 

26.7d 

374.0 

12.0 

h2o 

0.18 

13 

25.2d 

397.5 

10.0 

h2o 

0.16 

14 

28.2d 

355.4 

15.0 

h2o 

0.25 

15 

28.3d 

353.5 

11.0 

h2o 

0.33 

16 

9.8 

1022.0 

123 

°3 

0.05 

a  After  Nichols  (1975). 
h 

“NESR  =  Noise  Equivalent  Spectral  Radiance  in  mW/(m  sr 
cm-*). 

cNot  on  Flight  Model  1. 
dNot  used  in  this  investigation. 

American  Air  Defense  Command  (Colorado 
Springs)  satellite  tracking  facility  provides  us  with 
the  orbital  parameters  for  each  of  the  DMSP 
satellites,  which  are  input  to  the  computer  code  that 
predicts  the  satellite  overpass  times. 

In  addition  to  the  comparison  with  Dobson 
measurements,  comparisons  will  also  be  made  be¬ 
tween  the  MFR-derived  ozone  measurements  and 
total  ozone  data  from  TIROS  N  and  Nimbus  7. 


3.2  RADIATIVE  TRANSFER  AND 
RETRIEVAL  METHODOLOGY 

Several  of  the  many  investigations  required  to 
validate  and  document  the  physics,  mathematics, 
and  software  as  applied  in  the  development  of  the 
total  ozone  retrieval  model  implemented  in  early 
1978  have  been  completed  or  are  in  progress.  The 


•Special  Dobson  measurements  are  not  being  taken  al  F3 
overpass  times  because  its  overpasses  are  near  sunrise  and  sunset. 
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assumptions  made  and  shortcuts  taken  in  develop¬ 
ing  the  initial  methodology  and  several  items  point¬ 
ed  out  by  various  reviewers  are  also  being 
evaluated. 

Techniques  and  software  were  developed  to 
compare  line-bv-line  calculated  absorption  spectra 
against  laboratory  spectroscopic  measurements. 
The  laboratory  data  of  Dr.  Grvvnak  and  Dr  Burch 
at  Ford  Aerospace  and  Communications  Corp.  for 
the  15-uni  carbon  dioxide  bands  were  obtained,  and 
software  to  read  and  display  the  data  were 
developed.  A  version  of  the  line-by-line  software 
was  modified  to  produce  calculations  for  the 
homogeneous  paths  used  in  the  laboratory  in¬ 
vestigation  for  direct  high-resolution  comparisons. 
Initial  results  indicated  considerable  disagreement. 
A  spectral  shift  was  obvious  and  the  magnitude  of 
the  absorption  was  somewhat  lower  for  the 
theoretical  calculations.  Discussions  with  Dr. 
McClatchey  (Air  Force  Geophysics  Laboratory) 
and  Dr.  Burch  revealed  that  similar  disagreements 
had  also  been  found  in  their  studies.  Dr. 
McClatchey  felt  that  the  position  of  lines  and  bands 
in  the  spectrum  were  quite  good  in  his  theoretical 
line  parameter  data  that  are  used  as  input  to  our 
line-by-line  calculations.  Our  analysis  indicated  that 
the  spectral  shift  was  probably  due  to  some  instru¬ 
ment  error  in  the  spectroscopic  data.  Further  dis¬ 
cussions  with  Dr.  Burch  lead  to  the  conclusion  that 
the  laboratory  spectrum  was  observed  through  a 
“skewed"  instrument  slit  function. 

Software  was  developed  to  degrade  our  line-by- 
line  calculations  with  a  skewed  slit  function.  Subse¬ 
quent  displays  of  the  data  using  slit  function  widths 
appropriate  for  the  laboratory  data  indicated  that 
this  was  probably  the  cause  for  the  apparent  shift, 
since  good  agreement  was  obtained  for  the  positions 
of  lines  and  bands. 

The  disagreement  in  absorption  magnitude  has 
been  investigated  by  adjusting  two  different  input 
parameters  which  might  be  the  source  of  error.  One 
possible  source  of  error  in  the  laboratory  data  is  in 
establishing  the  amount  of  carbon  dioxide  in  the  op¬ 
tical  path  of  the  spectrometer.  An  adjustment  of  ap¬ 
proximately  +  8.6%  in  the  reported  concentrations 
produces  results  that  are  greatly  improved  Dr. 
Burch  felt  that  an  error  that  large  was  highly  un¬ 
likely.  He  suggested  that  a  similar  effect  would 
result  from  a  theoretical  half-width  error  for  the  ab¬ 
sorption  lines,  and  that  these  parameters  were  not 
well  known. 


Calculations  using  a  +8%  adjustment  in  the 
half-widths  produces  results  that  give  good  agree¬ 
ment  with  the  laboratory  data.  This  was  true  for  the 
data  that  were  observed  through  a  skewed  slit  func¬ 
tion  and  for  later  higher-resolution  data  that  were 
properly  taken. 

The  results  of  our  efforts  were  made  available 
to  Drs.  Burch  and  McClatchey,  and  further  discus¬ 
sions  will  be  held  in  order  to  reach  an  agreement  on 
the  most  probable  cause  for  the  discrepancies.  If 
these  questions  can  be  resolved  within  the  next  few 
months,  new  line-by-line  carbon  dioxide  calcula¬ 
tions  will  be  carried  out  including  the  necessary  ad¬ 
justments  before  regular  data  processing  begins. 

Comparisons  similar  to  those  for  the  carbon 
dioxide  bands  are  also  being  made  for  the  ozone 
bands  in  the  near  future.  The  laboratory  data 
available  for  the  ozone  bands  are  not  of  high  resolu¬ 
tion,  and  the  fine  details  of  the  spectra  cannot  be 
compared.  However,  the  low-resolution  absorption 
magnitudes  will  still  be  useful  for  evaluating  the  ac¬ 
curacy  of  our  calculations.  The  comparison  is  ac¬ 
complished  by  degrading  the  high-resolution  line- 
by-line  calculations  by  numerically  applying  an  ap¬ 
propriate  spectrally-weighted  instrument  slit  func¬ 
tion.  Initial  comparisons  for  ozone  also  indicate  dis¬ 
crepancies  in  results  for  the  9.6-Mm  bands. 

The  homogeneous  path  line-by-line  calcula¬ 
tions  and  laboratory  measurement  comparisons  are 
critical  evaluations,  since  measurements  for  in¬ 
homogeneous  paths  are  very  limited  and  of  little  use 
(because  exact  temperature,  pressure,  and  concen¬ 
tration  specifications  are  not  well  known). 
Transmittances  for  atmospheric  paths,  therefore, 
must  be  determined  from  accurate  homogeneous 
calculations  appropriately  applied  for  atmospheric 
conditions  specified  during  the  simulation  process. 

The  final  step  in  the  verification  of  the  spectral 
calculations  is  to  compare  simulated  space  measure¬ 
ments  against  actual  satellite  radiometer  measure¬ 
ments.  Some  data  have  already  been  collected  for 
this  purpose.  Also,  software  is  under  development 
to  perform  the  necessary  calculations  and  to  analyze 
the  results. 

We  are  preparing  to  process  45  days  of  MFR 
data  to  compare  our  ozone  data  with  similar  data 
for  the  same  period  using  the  TIROS  N  and 
Nimbus  7  sensors.  The  data  period  for  comparison 
covers  January  I  to  mid-February  1979.  In  prepara¬ 
tion  for  this  comparison,  several  codes  required  for 
generating  the  simulated  radiances  and  total  ozone 
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retrieval  models  are  being  modified  to  implement 
changes  to  improve  their  accuracy  and  com¬ 
putational  efficiency.  For  example,  the  running  time 
of  the  codes  used  to  generate  simulated  radiances 
has  been  reduced  more  than  50%.  Additionally,  re¬ 
cent  ozonesonde  data  have  been  obtained  from  the 
World  Ozone  Data  Center  (Downsview,  Canada)  in 
order  to  bring  our  ozonesonde  data  base  up  to  date. 
The  ozonesonde  data  are  used  as  input  to  the 
simulated  radiance  calculations.  The  ozonesonde 
data  base  includes  nearly  6000  soundings  from 
which  representative  sets  of  data  are  selected  for  the 
season  for  which  the  radiances  are  simulated. 

Software  has  been  developed  to  perform 
various  analyses  on  retrieval  results.  Extraction  of 
various  parameters  and  data  for  specific  regions  and 
times  from  existing  archives  is  now  possible.  These 
extractions  can  now  be  printed  and  displayed,  and 
various  manipulation  options  are  being  developed 
for  performing  analyses  of  various  retrieval 
problems. 

3.3  DATA  BASING  AND 
DATA  PROCESSING 

The  major  task  for  the  past  year  has  been  the 
development  of  an  automated  data  basing,  archiv¬ 
ing,  and  processing  system.  The  computing 
capability  at  the  Atmospheric  Release  Advisory 
Capability  (ARAC)  center  at  LLL  will  be  used  for 
data  basing  and  data  processing.  We  have 
purchased  some  additional  hardware  to  supplement 
the  capabilities  of  the  DEC  VAX  1 1 /780  computer 
system  which  became  operational  in  June  1979  at 
the  ARAC  center.  A  seven-track  tape  drive  was 
added  to  enable  the  input  data  tapes  from  AFGWC 
to  be  handled,  and  a  floating  point  accelerator  was 
added  to  enhance  the  computation  speed  for  MFR 
data  processing.  A  176  megabyte  disc  system  was 
added  to  increase  the  online  data  storage  capacity. 
One  disc  can  hold  one  month  of  satellite  and 
processed  data,  which  is  a  convenient  size  for 
building  the  long-term  archival  data  base  and  for 
data  analysis  purposes. 

High  density  tape  (6250  bpi)  has  been  selected 
as  the  long-term  storage  medium  for  the  data  base. 
Disc  packs  had  initially  been  considered  as  the 
storage  medium,  but  they  were  rejected  when  in¬ 
vestigation  showed  that  there  is  a  serious  risk  of 
catastrophic  disc  failure  and  head  crashes  when  the 


disc  packs  are  handled  frequently.  The  tape  medium 
provides  considerably  less  vulnerability,  can  be  in¬ 
expensively  duplicated,  is  cost  effective  over  the 
data  collection  period  even  with  the  necessary  ad¬ 
ditional  hardware,  and  provides  only  slightly  in¬ 
ferior  accessibility  to  the  data. 

The  presently  available  1600  bpi  tape 
technology  will  be  used  through  the  software 
development  phase  of  this  effort  since  the  6250  bpi 
tape  technology  has  only  recently  become  available 
for  use  on  the  VAX  11/780  computer  system.  The 
high  density  tape  drive  has  been  ordered  and  will 
arrive  before  regular  data  processing  is  scheduled  to 
begin. 

An  LLL-developed  scientific  data  management 
system  called  FRAM1S  has  been  chosen  for  use 
with  the  SOAC/DMSP  data  base.  This  is  a 
“relational  algebra"  approach  at  the  forefront  of 
current  data  base  technology.  LLL’s  Computation 
Department  has  already  invested  in  a  multi-man- 
vear  effort  to  develop  this  capability  and  plans  con¬ 
siderable  expansion  and  enhancement  in  the  coming 
year.  The  system  is  currently  available  on  the  CDC 
7600  computers,  and  we  have  completed  most  of  the 
modifications  necessary  to  transfer  FRAMIS  to  the 
VAX  1 1/780  computer  system.  Testing  of  the  initial 
capability  is  scheduled  to  begin  in  December  1979. 
Nearly  300  computer  routines  are  being  rewritten  to 
run  on  the  new  system. 

The  codes  that  process  and  calibrate  the  raw 
DMSP  radiometer  data  are  being  converted  to  the 
VAX  11/780  system.  Part  of  the  conversion  in¬ 
volves  the  word  size  dependency  problem 
(AFGWC  =  36  bits,  LLL  7600  =  60  bits,  and  VAX 
11/780  =  32  bits),  and  part  involves  the  develop¬ 
ment  of  a  new  quality  control  package  to  evaluate 
the  sensor  data  stream  continuously  while  process¬ 
ing  the  radiometer  data  and  to  provide  a  record  of 
that  performance  for  the  data  base  archive  quality 
assurance. 

Significant  progress  has  been  made  in  how 
total  ozone  data  are  mapped  and  contoured.  We 
completed  the  development  of  a  computer  code  that 
generates  an  “intelligent  fill"  in  data-void  areas  of 
the  daily  global  ozone  analysis.  The  code  generates 
the  fill-field  by  first  mapping  and  averaging  the 
satellite  ozone  retrievals  in  2-deg  latitude-longitude 
bins.  The  bins  remaining  empty  after  24  hours  of 
satellite  data  mapping  require  filling.  Filling  is  ac¬ 
complished  by  one  of  the  follow  ing  methods,  which 
depends  on  the  data  available.  The  preference  for 


each  method  decreases  in  the  order  listed:  (1)  inter¬ 
polate  between  the  binned  data  of  the  day  preceding 
and  the  day  following,  (2)  use  the  binned  data  from 
the  day  preceding  or  the  following,  (3)  use  the  zonal 
average  values  that  are  computed  from  the  binned 
data  for  that  day,  or  (4)  use  zonal  average  values 
that  are  taken  from  known  climatology. 

The  resulting  global  fill- Held  is  not  combined 
with  the  binned  data  until  it  is  smoothed.  It  is 
smoothed  by  applying  an  elementary  filter  that  ap¬ 
proximates  a  bivariate  normal  function  three  suc¬ 
cessive  times.  After  each  pass  of  the  filter,  the 
filtered  fill-field  is  overlayed  with  the  2-deg  binned 
data  for  that  day.  The  overlaying  and  filtering 
process  blends  the  fill-field  in  the  data-void  areas 
with  the  binned  data  outside  the  data-void  areas  for 
that  day.  Also,  the  filtering  process  removes  small- 
scale  features  (e  g.,  5-deg  horizontal  wavelength  and 
smaller)  but  retains  large-scale  features  (e.g.,  a 
wavelength  of  60  deg  is  retained  at  96%  of  its 
original  amplitude)  in  the  data-void  areas.  Thus 
only  the  large-scale  features  that  change  little  in  am¬ 
plitude  and  position  from  day  to  day  are  used  as  an 
“intelligent  fill.” 

The  generated  fill-field  may  be  used  with  either 
the  spline  analysis  code  or  the  newly  developed  filter 
analysis  code,  which  are  used  to  generate  ozone 
contour  maps.  In  either  code,  the  fill-field  is  used 
only  to  fill  the  data  void  regions  of  the  global  data¬ 
set. 

The  methodology  of  the  filter  analysis  code  is 
essentially  to  compute  an  average  grid  point  value 
in  each  2-deg  latitude-longitude  bin,  fill  the  data- 
void  bins  using  the  fill-field,  and  filter  the  gridded 
data  using  a  Lanczos  filter  (Duchon,  1979)  having  a 
preselected  cutoff  frequency.  The  methodology  of 
the  spline  analysis  is  to  divide  the  globe  into  10-deg 
latitude-longitude  regions,  fit  the  data  within  a 
region  in  a  least-squares  sense  with  B-splines  assur¬ 
ing  continuity  across  region  boundaries,  and 
treating  filled  portions  of  the  region  as  data  points. 

A  study  of  the  two  methods  has  shown  that  the 
filter  method  requires  85%  less  computer  time  and 
fits  the  input  data  better  than  the  spline  method  for 
all  filter  cutoff  wavelengths  (or  cutoff  frequency) 
smaller  than  20-deg  latitude-longitude.  The 
smoothness  of  the  analysis,  generated  using  the 
filter  method,  can  be  varied  by  varying  the  cutoff 
wavelength;  increasing  the  cutoff  wavelength  in¬ 
creases  the  smoothness.  Two  examples  of  analyses 
with  cutoff  wavelengths  (Ac)  of  10  and  16  deg  are 


shown  in  Figs.  15  and  16,  respectively.  The  shaded 
areas  of  each  map  indicate  data-void  regions.  The 
analysis  in  Fig.  16  is  aesthetically  more  pleasing 
than  Fig.  15,  but  has  less  information  content  than 
the  analysis  in  Fig.  15.  We  experimented  with  a 
number  of  values  of  Ac  ranging  from  6  to  20  deg  and 
selected  10  deg  as  a  balance  between  aesthetics  and 
information  content. 

We  have  also  developed  codes  that  generate 
analyses  at  2. 5-deg  latitude-longitude  grid  intervals 
matching  the  National  Meteorological  Center 
(NMC)  global  temperature/pressure-height  grid. 
This  analysis  requires  35%  less  computer  time  than 
the  2-deg  analysis  using  the  filter  method.  The  2.5- 
deg  analysis  is  essentially  the  same  as  the  2-deg 
analysis  except  for  a  modulation  of  the  minimum 
and  maximum  (L  and  H)  values  on  the  order  of  1%, 
Figure  17  is  an  example  of  the  2. 5-deg  analysis  using 
Ac  =  10  deg. 

Our  plan  is  to  process  the  45  days  of  data 
(January  1-February  15,  1979)  for  comparison  with 
data  from  TIROS  N  and  Nimbus  7  using  the  filter 
method  at  2-deg  latitude-longitude  grid  point  inter¬ 
vals.  When  regular  data  processing  begins,  we  will 
process  the  data  at  2. 5-deg  grid  intervals  using  the 
filter  method. 


3.4  DATA  ANALYSIS 

During  the  past  year  the  number  of  Dobson 
observatories  that  have  agreed  to  take  special  ozone 
measurements  increased  from  33  to  37.  Also, 
measurements  from  the  second  Dobson  instrument 
(which  is  automated)  at  Arosa  are  being  sent  to  us. 
The  locations  of  the  participating  Dobson  obser¬ 
vatories  are  shown  in  Fig.  18.  Special  observations 
are  made  daily  at  the  time  of  DMSP  satellite  over¬ 
passes.  The  data  that  we  receive  are  entered  into  a 
Dobson  ozone  data  base  for  comparison  with  the 
MFR-derived  ozone  data.  The  Dobson  data  base 
extends  from  Spring  1977  to  the  present.  All  of  the 
Dobson  data  that  we  receive  are  also  forwarded  to 
the  NOAA  Satellite  Physics  Branch,  Silver  Spring, 
Maryland,  where  they  are  in  turn  sent  to  other 
NOAA  and  NASA  users.  Quality  assurance  codes 
were  developed  to  check  the  Dobson  data  for  any 
obvious  or  suspected  errors.  Software  was  also  writ¬ 
ten  to  provide  information  in  a  graphical  format  on 
the  availability  of  ozone  data  in  the  Dobson  data 
base  for  each  individual  station. 
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21,  1977  at  2-deg  grid  interval  ( Ac=  10° ). 
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Analysis  for  June  21, 1977  at  2-deg  grid  interval  (Ac=  16' 


leg  grid  interval  (Xc=  10°). 
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FIG.  18.  Location  of  Dobson  spectrophotometer  observatories  where  special  measurements  of  total  ozone  are 
taken  for  comparison  with  MFC-derived  total  ozone  data. 


Preliminary  analyses  have  been  made  between 
the  MFR  ozone  data  that  we  processed  during  the 
feasibility  study  and  conventional  meteorological 
data  fields.  Upper-atmospheric  data  of  tem¬ 
perature,  pressure-height,  wind  velocity  and 
humidity  were  obtained  from  the  National 
Meteorological  Center,  Boulder,  Colorado.  Com¬ 
puter  codes  were  built  to  access  the  upper- 
tropospheric  and  stratospheric  data.  Black  and 
white  and  color  graphics  routines  were  used  to  dis¬ 
play  multi-parameter  variability.  A  polar  map  pro¬ 
jection  is  currently  being  used  for  displaying  the 
data. 

A  comparison  between  the  spatial  variabilities 
of  total  ozone  and  100-mb  temperatures  is  shown  in 
Fig.  19  for  the  Southern  Hemisphere.  The  polar 
projection  extends  from  the  South  Pole  to  20°S 
latitude.  The  ozone  data  are  the  average  of  the  13 
days  of  MFR  data  that  have  been  processed  from 
June  1977.  The  Southern  Hemisphere  was  chosen 
for  comparison  because  it  is  dynamically  more  ac¬ 
tive  than  the  Northern  Hemisphere  at  this  time  of 
year. 

The  total  ozone  isolines  in  middle  latitudes  in¬ 
dicate  a  dominance  of  wave  number  five.  The  100- 
mb  temperature  data  indicate  a  dominance  of  wave 
number  five  or  six.  At  several  locations  the  total 


ozone  and  temperature  waves  are  out  of  phase  with 
each  other.  Examples  of  this  may  be  seen  over 
South  America  and  near  Africa.  At  these  locations, 
relatively  higher  (lower)  amounts  of  total  ozone  are 
associated  with  warmer  (colder)  lower-stratospheric 
temperatures.  Total  ozone  maxima  occurring  near 
longitudes  of  10°E.  170°E  and  80°W  at  50°S 
latitude  are  located  near  cores  of  warm 
stratospheric  air.  The  relationships  between  tem¬ 
perature  and  total  ozone  just  described  are  not  ap¬ 
parent  at  all  locations  in  the  map.  The  region  of 
high  total  ozone  south  of  Australia  corresponds 
well  with  the  location  of  high  ozone  values  for  June 
indicated  by  London  et  al.  (1976)  for  data  from 
1957-1967.  The  MFR  total  ozone  data  show  more 
structure  in  this  feature  than  is  indicated  by  London 
et  al.  ( 1976). 

Software  is  being  developed  that  will  be  used  to 
perform  a  power  spectral  analysis  of  the  ozone  field 
from  MFR  data.  Also,  work  has  begun  on  code 
development  for  spherical  harmonic  analyses  of  the 
MFR  total  ozone  data.  This  line  of  research  is  one 
of  several  to  be  followed  in  later  months.  A  number 
of  methods  of  discerning  relationships  among 
various  atmospheric  variables  are  being  tested  and 
analyzed.  Thus  far  the  methods  have  been  applied 
to  a  three-day  period  (in  June  1977)  in  the  Southern 
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FIG.  19.  Southern  Hemisphere  June  1977  total 
ozone  (units:  M.atm.cm;  solid  line)  from  satellite 
MFR  and  100-mb  temperature  (°C;  dashed  line) 
from  radiosonde. 


Hemisphere.  Several  promising  graphical  display  should  prove  to  be  very  useful  when  we  begin 
and  data  analysis  tools  have  been  developed  which  regular  data  processing. 


4.  WORK  IN  PROGRESS 


We  are  in  the  process  of  developing  a  two- 
dimensional  transport-kinetics  model  which  will  be 
in  the  programming  and  initial  validation  stage  dur¬ 
ing  most  of  Fiscal  Year  1980.  This  model  will  be 
used  to  study  latitudinal  and  seasonal  effects  of 
various  atmospheric  perturbations. 

A  steady  state  version  of  the  one-dimensional 
transport-kinetics  code  is  under  development  which 
will  significantly  reduce  the  computation  time 
needed  for  steady-state  calculations.  Most  of  our 
sensitivity  studies  involve  computing  the  change  in 
total  ozone  at  steady  state,  so  this  code  will  be  used 
extensively  in  the  future. 


Routine  processing  of  the  DMSP  MFR  data 
will  begin  in  the  spring  of  1980  and  will  extend 
through  Fiscal  Year  1981.  In  all,  data  from  four 
satellites  with  overlapping  coverage  over  a  three- 
year  period  will  be  processed.  Daily  maps  of  total 
ozone  and  monthly  statistics  and  comparisons  with 
Dobson  data  will  be  published  for  distribution  to 
the  scientific  community.  We  will  also  be  con¬ 
ducting  studies  of  the  spatial  and  temporal 
variability  of  ozone  related  to  natural  and  man¬ 
made  influences. 
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APPENDIX  A 

DESCRIPTION  OF  LLL  ONE-DIMENSIONAL 
TRANSPORT-KINETICS  MODEL 
(AS  PRESENTED  AT  NASA  HARPERS  FERRY  WORKSHOP  1979) 

The  local  concentrations  of  trace  species  in  an  air  parcel  are  determined  by  the  chemical  and 
photochemical  processes  and  nonchemical  sources  and  sinks  occurring  within  the  parcel,  and  by  the  transport 
and  radiative  (luxes  into  and  out  of  the  parcel.  One-dimensional  mathematical  models  of  the  chemical 
processes  in  the  stratosphere  are  governed  by  the  chemical  species  conservation  equation 


Scj/dt  +  F^zd)  =  Pj[c,J(z,t,c),  k(z,t).p)j  -  Lj(cJ(z,t,c)  k[T(z,t),p|}Cj  +  Sj(z,t) 


(A-l) 


where  c,  =  c,(z,t)  is  the  concentration  of  the  ilh  chemical  constituent;  c  is  the  general  representation  of  all  con¬ 
stituents;  P,  and  L,C|  are  the  production  and  loss  of  Cj  due  to  photochemical  interactions;  T  is  the  ambient  air 
temperature;  />  is  the  ambient  air  density;  F,  is  the  vertical  transport  flux  of  c,;  Sj  represents  any  other  possible 
sinks  or  sources  of  c,;  J  represents  photodissociation  coefficients;  k  represents  chemical  reaction  rate  coef¬ 
ficients;  and  all  of  these  variables  are  defined  at  a  given  vertical  position  z  at  time  t.  The  explicit  display  of  the 
major  interdependent  relations  of  these  variables  in  Eq.  (A-l)  illustrates  the  nonlinearity  and  general  complex¬ 
ity  of  this  mathematical  system.  Equation  (A-l )  is  defined  over  a  spatial  domain  D  with  appropriate  boundary 
conditions.  There  is  one  conservation  equation  for  each  chemical  species  treated. 


PHYSICAL  DOMAIN 

The  LLL  one-dimensional  model  extends  from  the  ground  to  56.25  km.  The  model  currently 
calculates  the  vertical  concentration  distributions  of  39  (2  of  which  are  used  only  in  sensitivity  studies)  at¬ 
mospheric  trace  constituents.  The  model  contains  134  (14  of  which  are  used  only  in  sensitivity  studies) 
chemical  or  photochemical  reactions.  Table  A- 1  lists  the  species  solved  for  in  the  model.  Of  these  species, 
O('D),  H,  and  N  are  assumed  to  be  in  instantaneous  equilibrium.  The  vertical  grid  structure  is  variable,  but 
for  the  calculations  reported  here,  we  have  a  0.5-km-thick  layer  at  the  surface,  1-km  thick  layers  extending 
from  0.5  to  34.5  km,  a  1.75-km  thick  layer  between  34.5  and  36.25  km  and  2.5-km  thick  layers  extending  to 
56.25  km. 


TRANSPORT  REPRESENTATION 

The  net  vertical  transport  flux,  Fp  of  any  minor  constituent  Cj  is  represented  through  a  diffusion  ap¬ 
proximation  in  which  Fj  is  assumed  to  be  proportional  to  the  gradient  of  the  mixing  ratio  of  that  trace  species: 

Fi  =  'Kip  a4  (ci/p)  (A'2) 

where  K,  is  the  one-dimensional  vertical  diffusion  coefficient.  The  LLL  K,  profile  (shown  in  Fig.  A- 1 )  used  in 
the  calculations  for  this  report  was  originally  based  on  an  analysis  of  N^O  and  CF^  measurements  (NAS, 
1976)  with  considerations  also  given  to  measurements  of  radionuclide  debris  transport  in  the  lower 
stratosphere. 

The  LLL  one-dimensional  model  has  been  designed  such  that  profiles  of  K,  utilized  by  other  groups 
(or  at  previous  times)  can  be  easily  incorporated.  Such  profiles  have  been  utilized  to  test  the  sensitivity  of  the 
results  to  transport  parameterization  uncertainties. 
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TABLE  A-l  •  Species  calculated  in  LLL  one-dimensional  model. 


ot3P) 

Cl 

CH4 

°3 

CIO 

HCO 

NO 

CIONOj 

CHjO 

no2 

cino2 

CH3 

n2o 

HCI 

CHjOOH 

HNO3 

OCIOa 

CHjO 

OH 

HOCI 

ch302 

ho2 

CH3CI 

CO 

h2°2 

CF2d2 

H2 

NO  3 

CFCI3 

O('d) 

N2°5 

CH3CCI3 

N 

h2o 

CIOj3 

H 

HONO 

CCI4 

HNO . 

4 

aUsed  only  in  sensitivity  studies. 


CHEMISTRY 

We  have  used  two  1979  versions  of  model  chemistry  in  this  report  (see  Tables  A-2  through  A-4). 
1979a  chemistry  was  based  primarily  on  the  rate  recommendations  in  JPL  (1979).  However,  several  reactions 
discussed  in  JPL  (1979)  are  omitted  in  the  model,  and  several  reactions  not  discussed  in  JPL  (1979)  are  in¬ 
cluded.  This  chemistry  was  used  for  many  sensitivity  calculations  carried  out  in  the  spring  of  1979.  1979b 
chemistry  was  based  almost  exclusively  on  the  draft  chapter  on  chemical  reaction  rates  prepared  at  the  NASA 
Harpers  Ferry  Workshop  (June,  1979),  It  has  a  few  comparatively  minor  differences  from  the  final  draft  of 
that  report,  and  it  includes  a  few  reactions  not  assessed  at  the  NASA  Workshop. 

Two  reactions  treated  in  the  1979a  chemistry  are  controversial  and  are  of  some  importance  to  our 
sensitivity  studies.  These  are 

M 

CIO  +  NO,  ™  CIONO,  (A3) 


and 


CIO  +  OH  -*•  HCI  +  02  .  (A4) 

In  the  1979a  chemistry,  we  used  the  slower  of  the  two  JPL-recommended  rate  constants  for  chlorine  nitrate 
formation,  and  we  adopted  an  expression  for  HCI  formation  from  OH  +  CIO  that  is  about  half  the  upper 
limit  for  that  reaction  path.  In  the  1979b  chemistry,  we  used  the  faster  of  the  two  recommendations  for  the 
chlorine  nitrate  formation  rate  coefficient.  Three  considerations  inspired  this  choice:  (I)  the  majority  of  the 
chemistry  panel  seemed  to  favor  the  faster  expression,  (2)  even  if  the  bulk  of  the  reaction  between  CIO  and 
NO2  leads  to  other  products  (as  suggested  by  those  favoring  the  slower  rate  coefficient),  the  other  products 
might  easily  have  an  effect  on  stratospheric  chemistry  similar  to  that  of  CIONOi.  and  (3)  it  improved  the  com¬ 
parison  between  calculation  and  observation  for  both  CIO  and  CIONOi. 
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In  our  1979b  chemistry  we  omitted  the  reaction  forming  HC1  from  OH  +  CIO.  The  decision  to  in¬ 
clude  this  reaction  in  the  1979a  chemistry  was  based  on  privately  communicated  preliminary  results  that 
seemed  to  suggest  that  the  reaction  probably  occurred,  but  they  have  since  been  interpreted  as  providing  only 
an  upper  limit.  Both  of  these  controversial  choices  of  rate  coefficient  have  a  significant  impact  on  model  sen¬ 
sitivities  (especially  for  C1X)  but  they  oppose  each  other.  As  a  result,  1979a  and  1979b  chemistries  yield 
qualitatively  similar  perturbational  sensitivities  for  both  NOx  and  CIX  perturbations.  The  1979b  chemistry  is 
less  controversial  than  the  1979a  chemistry  and  is  to  be  preferred  for  purposes  of  comparison  with  other 
workers.  All  primary  assessments  have  been  repeated  using  1979b  chemistry.  However,  several  sensitivity 
studies  were  not  repealed,  since  it  seemed  unlikely  that  the  qualitative  results  of  those  sensitivity  studies  would 
differ  if  they  were  repeated,  and  because  the  1979a  chemistry  is  well  within  the  limits  of  reasonable  uncertainty 
in  our  present  knowledge  of  the  atmosphere.  The  two  chemistries  are  also  useful  in  emphasizing  the  existence 
of  processes  for  which  no  clear  recommendation  is  available. 
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TABLE  A-2.  Chemical  reactions  and  rate  coefficients  where  k  =  A  eB  T  used  in  1979  model  chemistry. 


Reaction 

A* 

B 

Note 

I. 

M 

o  +o2-  o3 

See  Table  A-3 

- 

3 

2. 

0+03->  202 

1.5  X  10'11 

-2218 

1 

3. 

o3  +  no-  N02  +  02 

2.3  x  10'12 

-1450 

1 

4. 

0  +  NOj  -  NO  +  02 

9.3  x  10'12 

0 

1 

5. 

n2o  +  o(‘d)-  n2  +  o2 

4.8  x  10'11 

0 

1 

5.1  X  10'11 

0 

2 

6. 

NjO  +0(^D)-*  2NO 

6.2  x  10‘U 

0 

1 

5.9  X  10  * 1 

0 

2 

7. 

N+02  -  NO  +0 

4.4  X  10'12 

-3220 

1 

8. 

N  +  NO  -*  Nj  +  0 

3.4  X  10'11 

0 

1 

9. 

O('d)  +  HjO  -  20H 

2.3  X  1010 

0 

1 

10. 

o3  +  oh  -»  ho2  +  o2 

1.6  X  1012 

-940 

I 

11. 

O  +  OH  ->  02  +  H 

4.0  X  101 1 

0 

1 

12. 

03  +  H02  -  OH  +  202 

1.1  X  10'14 

-580 

1 

13. 

0  +  HOj  -  OH  +  02 

3.5  X  10-“ 

0 

1 

14. 

h  +  o2  ->  ho2 

See  Table  A-3 

- 

3 

IS. 

o3  +  h  -*  oh  +  o2 

1.4  X  10"10 

-470 

1 

16. 

HOj  +  HOj  -  H2Oj  +  02 

2.5  X  1012 

0 

1 

17. 

ho2  +  oh  -  H20+02 

4.0  X  1011 

0 

1 

18. 

OH  +  NOj  ”  HNOj 

See  Table  A-3 

- 

1 

19. 

OH  +  HNOj  -  HjO  +  NOj 

8.5  X  10  "14 

0 

1,10 

20. 

HjOj  +  OH  -  HjO  +  HOz 

1.0  X  1011 

-750 

1 

21. 

Nj  +  O(^D)  ”  N20 

See  Table  A-3 

- 

3 

22. 

N  +  N02  -*  N20  +  0 

2.1  X  10  11 

-800 

1 

23. 

NO  +  0  -*  NOj 

See  Table  A-3 

- 

3 

24. 

NO  +  HOj  -  N02  +  OH 

4.3  X  10'12 

200 

1 

3.4  x  1012 

250 

2 

25. 

H2+0(*D)  -  OH  +  H 

9.9  X  10  “ 

0 

1 

26. 

OH  +  OH  -»  HjO+O 

1.0  x  10"“ 

-500 

1 

27. 

N  +  Oj  -  NO  +02 

2.0  x  10  “ 

-3000 

4 

28. 

no2+o3  -  no3+o2 

1.2  X  10"13 

-2450 

1 

29. 

OH  +  OH  ^  HjOj 

See  Table  A-3 

3 

30. 

H202  +  O  -  OH  +  H02 

2.8  X  I0"12 

-2125 

1 

31. 

00  + OH  ”  H  +  COj 

See  Table  A-3 

- 

3 

•When  two  entries  ire  given,  the  lower  one  corresponds  to  1979.  chemistry  and  the  upper  one  corresponds  to  1979b  chemistry. 
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TABLE  A-2.  (Continued) 


Reaction 

A* 

B 

Note 

32. 

0('d>  +  M  -  O  +  M 

2.2  x  10‘M 

99 

5 

33. 

ci  +  o3  -  cio  +o2 

2.8  x  I0'u 

-257 

1 

34. 

a  +  no2  ”  CINOj 

See  Table  A-3 

- 

3 

35. 

CIO  +  0  -*  ci  +  o2 

7.7  X  10'11 

-130 

1 

36. 

NO  +  CIO  -  NOj  +  Cl 

7.8  X  10'12 

250 

1 

37. 

ao  +  NOj  ^  CIO  NO  2 

See  Table  A-3 

- 

3 

38. 

HCI  +  0('d)  -  Cl  +  OH 

1.4  X  10'10 

0 

1 

39. 

OH  +  HC1  -  H20  +  Cl 

2.8  X  10'12 

-425 

1 

40. 

0  +  HQ  -»  OH  +  Cl 

1.14  x  10  11 

-3370 

1 

41. 

a  +  H02  -  HCI  +  02 

4.5  X  10'11 

0 

1 

42. 

CFClj  +  O(^D)  -  3a 

2.2  X  10'10 

0 

1,6 

43. 

CF2CI2  +0(*D)  -  2 Cl 

1.4  X  10'10 

0 

1,6 

44. 

a  +  H2  -  HCI  +  H 

3.5  X  10'H 

-2290 

1 

45. 

a  +  h2o2  -  hci  +  ho2 

8.9  X  10'12 

1.7  X  10'12 

-925 

-384 

11 

2 

46. 

aONOj  +  O  -  CIO  +  NOj 

3.0  X  1012 

-808 

1,10 

47. 

CHjCl  +  OH  -»  Cl  +  HjO  +  H02 

2.2  X  10'12 

-1142 

1 

48. 

NO  +  NOj  -  2N02 

2.0  X  10')' 

8.7  X  IO'12 

0 

0 

1 

3 

49. 

N02  +  0  ”  N03 

See  Table  A-3 

- 

3 

50. 

N02  +  NOj  ”  NjOj 

See  Table  A-3 

- 

3 

51. 

NjOj  -  N02  +  NOj 

See  Table  A-3 

- 

3 

52. 

N2°5  +  H2°  2HN03 

1.0  X  10'20 

0 

7 

S3. 

0(‘d)  +  03  -  202 

1.2  X  10'10 

0 

1 

54. 

HOj  +  H02  +  HjO  -  HjOj  +  02  +  H20 

See  Table  A-3 

- 

3 

55. 

0  +  NOj  ->  Oj  +  NOj 

1.0  X  10-» 

0 

1 

56. 

ho2  +  no2  ”  hno4 

See  Table  A-3 

- 

3 

57. 

HN04  ^  HOj  +  NOj 

See  Table  A-3 

- 

3 

58. 

OH  +  HN04  -  HjO  +  NOj  +  Oj 

6.0  X  IO'12 

-750 

1.8 

59. 

a  +  HN04  -  HCI  +  NOj  +  Oj 

3.0  X  10'12 

-300 

8 

60. 

HOj  +  CIO  -  Oj  +  HOCI 

7.0  X  IO'13 

500 

1 

61. 

a  +  HOCI  -  HCI  +  CIO 

3.0  X  IO'12 

-300 

8 

62. 

OH  +  Hoa  -  HjO  +  CIO 

3.0  x  10'12 

-800 

1 

*When  two  entries  are  given,  the  lower  one  corresponds  to  1979a  chemistry  and  the  upper  one  corresponds  to  1979b  chemistry. 

TABLE  A-2.  (Continued) 


»3.  0  +  HOC1  -  OH  +  CIO 

1.0  x  I0‘“ 

>4.  OH  +  CH4  -  CH3  +  H20 

2.4  X  10'12 

>5.  0  +  CH4  -  CHj  +  OH 

3.5  X  1011 

>6.  O^D)  +  CH4  -  CH20  +  H2 

1.0  x  IO11 

»7.  0('d)+CH4  -  CHj  +  OH 

1.3  x  1010 

>8.  CH4  +  Cl  -»  HCl  +  CH3 

9.9  X  1012 

>9.  Cl  +  CHjCl  -  H02  +  CO  +  2HCI 

3.4  x  1011 

0.  CHjOj  +  NO  -  N02  +  CHjO 

7.0  x  IO12 

8.0  x  10'12 

1.  a  +  ch2o  -  HCl  +  HCO 

9.2  x  tO'11 

2.  CH302  +  ho2  -  CH3OOH  +  o2 

6.0  X  10"12 

1.0  X  IO'12 

3.  CHjO  +  02  -  CH20  +  HOz 

5.0  X  10'13 

4.  OH  +  CHjO  -  HzO  +  HCO 

1.0  X  10  11 

1.7  X  10'11 

5.  0  +  CH20  -  HCO  +  OH 

3.2  X  10  “ 

2.8  X  10'11 

6.  hco  +o2  -  co  mo2 

5.0  X  10'12 

7.  OH  +  CH3OOH  -  CH3Oj  ♦  HzO 

5.0  X  10»2 

6.2  X  IO'12 

8.  CH3  +  O  -  CHjO  +  H 

1.0  x  10’10 

9.  CHjOj  +  Oj  -  CHjO  +  20j 

1.0  x  to'14 

0.  CH302  +  O  -  CHjO  ♦  02 

3.0  X  10  11 

i.  ao  +  oh  -  ho2  +  ci 

9.2  X  IO'12 

2.0  X  1011 

2.  CHj  +  02  -1  CHj02 

See  Table  A-3 

3.  CIO  +  OH  -  HCl  +  02 

Not  used  ("VlO"22) 
2.0  x  IO'12 

4.  H2  +  OH  -  HjO  +  H 

1.2  x  IO'11 

Not  used 

5.  H  +  H02  -  H2  +02 

4.2  X  10'M 

Not  used 

6.  OH  +  CHjOOH  -  CHjO  +  HjO  +  OH 

5.0  x  IO'12 
Not  used 

7.  O  v  HN04  -*  OH  +  N02  +  02 

1.0  x  IO'12 

Not  used 

8.  OH  +  CK)N02  -  HOCI  +  NOj 

1.2  X  IO"12 

Not  used 

*When  two  entries  ire  given,  the  lower  one  corresponds  to  1979*  chemistry  and  the  upper  one  corresponds  to  1979b  chemistry. 


TABLE  A-2.  (Continued) 


Reaction 

Aa 

B 

Note 

89.  Cl  +  CIONOj  -  2C1  +  NOj 

1.7  X  I012 

Not  used 

-607 

1,10 

90.  HONO  +  OH  -  H20  +  NOj 

6.6  X  1012 

Not  used 

0 

1 

91.  OH  +  NO  ^  HONO 

92-97.  (not  used;  reactions  used  only  in  sensitivity  studies) 

See  Table  A-3 

Not  used 

- 

3 

98.  0  +  OCIO  -  CIO  +  02 

2.5  x  10  1 1 

-1166 

1 

99.  NO  +  OCIO  ^  N02  +  CIO 

2.5  x  I012 

-600 

I 

aWhen  two  entries  are  given,  the  lower  one  corresponds  to  1979a  chemistry  and  the  upper  one  corresponds  to  1979b  chemistry. 


Notes  to  Table  A-2. 

1 .  Draft  report  of  NASA  Harpers  Ferry  Workshop. 

2.  JPL  (1979).  Where  only  one  entry  is  given  for  the  rate  coefficient,  1979a  and  1979b  are  the  same.  Isually  this  means  that 
the  first  two  references  at  the  end  of  this  Appendix  give  the  same  recommendation. 

3.  The  reaction  is  pressure-dependent.  See  Table  A-3  for  discussion. 

4.  Estimate  designed  to  be  compatible  with  upper  limit  given  in  the  first  reference  at  the  end  of  this  Appendix,  and  low  enough 
to  have  no  significant  effect  on  model  performance.  Reaction  is  retained  only  to  facilitate  reintroduction  if  the  evaluated 
upper  limit  should  prove  to  be  in  error. 

5.  Weighted  average  of  the  rates  of  Of1  D)  +  N2  and  Of'  D)  +  from  the  Fust  two  references. 

6.  Product  chemistry  has  been  simplified. 

7.  Estimated  reaction  rate.  This  estimate  is  designed  to  include  a  possible  heterogeneous  contribution  to  the  overall  reaction. 
Important  only  in  the  lower  troposphere. 

8.  Estimated  reaction  rate.  This  rate  is  estimated  based  on  the  assumption  that  UNO,  and  HOCT  resemble  H,0,  (as  treated  in 
JPL,  1979)  in  reactions  with  Cl  and  OH. 

9.  Estimated  reaction  rate.  Rate  is  estimated  based  on  the  assumption  that  CHj02  closely  resembles  HO2  in  reaction  with  O  or 
°3- 

10.  Products  are  not  given  in  the  first  two  references  at  the  end  of  this  Appendix.  The  assumed  products  are  based  on  the  products 
that  seem  most  plausible  based  on  chemical  considerations. 

1 1.  Rate  based  on  a  draft  of  the  first  reference  that  trivially  differs  from  the  final  draft. 

12.  1979a  chemistry  treated  the  reactions  of  HO  and  CIO  based  on  privately  communicated  qualitative  preliminary  results.  The 
treatment  is  nearly  an  upper  limit  to  the  plausible  rate  coefficients  based  on  the  recent  results  of  Leu  and  Lin  (1979). 


Our  treatment  of  photolysis  reactions  has  also  been  modified.  There  is  evidence  for  a  moderate  tem¬ 
perature  dependence  for  many  photoabsorption  cross  sections.  With  the  exception  of  ozone.  NO,  and  Oj 
photolysis,  we  have  not  treated  this  temperature  dependence  explicitly,  but  have  used  cross  sections  measured 
at  roughly  230  K  for  all  temperatures.  As  a  result  our  calculated  trace  species  photodissociation  rates  should 
be  more  accurate  for  the  stratosphere  than  for  the  lower  troposphere. 

For  ozone  photolysis  we  use  quantum  yields  based  on  the  recommendations  of  NASA  (1979).  Our 
treatment  of  Oj  photolysis  is  based  on  Hudson  and  Mahle  (1972)  while  our  treatment  of  NO  photolysis  is 
based  on  Frederick  and  Hudson  (1979). 
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TABLE  A-3.  Rate  coefficients  used  for  pressure-dependent  reactions. 


Expression  1 

A^MIUM/T)"0 

k  =  - - - 

1  +  A0(M1(300/T)n0/Ai(300/T)ni 


/A0(M](300/T)n0\ 

2 

1  + 

g'°  \  A.(300,T)n>  J 

Reaction 

V 

n0 

A. 

i 

n. 

Note 

ho2  +  no2  ”  hno4 

2.1  x  1031 

5.0 

6.5  x  10'12 

5.0 

1 

M 

oh  +  no2  -  HNOj 

2.6  x  10'30 

2.9 

2.4  x  10'11 

1.3 

1 

ao  +  NO,  ^  CIO  NO, 

1.6  x  10'31 

3.4 

1.5  X  10  11 

1.9 

1,2 

3.5  X  10‘32 

3.8 

1.5  x  10  * 1 

1.9 

1.2 

M 

0  +  °2  -  03 

6.2  X  10'34 

2.1 

- 

0 

1 

CHj  +  02  ”  CHj02 

2.2  X  10'31 

2.2 

2.0  x  10'12 

1.7 

1 

1  M 

OCDJ  +  N-  -*  N20 

3.5  x  10‘37 

0.45 

- 

0 

1 

a  +  no2  ^  cino2 

1.6  x  10'30 

1.9 

3.0  X  10'U 

1.0 

1 

M 

h  +  o2  -»  ho2 

5.5  X  10'32 

1.4 

- 

0 

1 

M 

OH  +  NO  -  HNOz 

6.7  X  10  31 

3.3 

3.0  x  10"11 

1.0 

1 

OH  +  OH  ”  H202 

2.5  x  10'31 

0.8 

3.0  X  10'11 

1.0 

1 

NO,  +  NO,  “  N,0. 

1.4  x  10'30 

2.8 

9.0  x  1013 

-0.7 

1 

1.8  x  10-32e1316^ 

0 

9.5X1013e58^ 

0 

3 

O  +  NO  ”  NO, 

1.2  x  10‘31 

1.8 

3.0  x  10'11 

-0.3 

1 

I^XlO^e584'7 

0 

- 

4 

O  +  NOj  ”  NOj 

9.0  X  10'32 

1.0  X  10'31 

2.0 

0 

2.2  X  10  * 1 

0 

0 

1 

4 

M 

N,0,  -  NO,  +  NO, 

usxio-V11180^ 

2.8 

7.52X1014e11180^ 

-0.7 

5 

,6xioV9864^ 

0 

7.94  X  I014e',n22 

0 

4 

hno4  “  ho2  +  no2 

5.2  X  lO-8  e10015^ 

7 

1  +  4.86  x  10'12  M0,61 

OH  +  CO  ^  H  +  CO, 

1.35  x  1013  A  +  - — 

k  = 

1 

\  2.46  x  10 19  J 

HOj  +  H02  +  H20  -  H202 

+  HjO  +  02 

k  = 

i,xi0-34e373°fr 

1  +  M  •  3.5  X  10'16  e'2060/T 

®When  two  entries  ire  given,  the  lower  one  corresponds  to  1979i  chemistry  and  the  upper  one  corresponds  to  1979b  chemistry. 
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Notes  to  Table  A-3 

1.  Expression  given  in  NASA  (1979).  1979a  chemistry  differed  from  1979b  chemistry  in  that  the  parameter  0.6  in  expression  (1) 
was  set  equal  to  0.8  in  the  expression  used  for  1979a  chemistry. 

2.  Both  expressions  are  recommended  with  no  clear  preference.  The  lower  value  is  used  in  the  1979a  chemistry  and  the  upper 
one  is  used  in  the  1979b  chemistry. 

3.  Private  communication,  D.  Golden,  SRI  International,  Inc.  (1979). 

4.  Based  on  data  in  NBS  513  (1978). 

5.  Based  on  data  in  the  first  reference  and  the  equilibrium  constant  from  NBS  5 13  (1978). 

6.  Based  on  Cox  (1978). 

7.  Based  on  Graham  er  a/.  (1978). 


TABLE  A-4.  Photolysis  reactions.  Alternative  products  of  reaction  are  shown  in  parentheses,  but  they  were  not 
used  in  either  the  1979a  or  1979b  chemistry. 


Reaction 

Note 

Reaction 

Note 

1.  02  -  20 

1,2 

16. 

CF2C12  -  2CI 

4,8 

2.  Oj  -  0+02 

1,3 

17. 

CFClj  -  3C1 

4,8 

3.  03  -  0('d)  +  02 

1,3 

18. 

CC14  -  4CI 

4,8 

4.  NOj  -  NO  +  O 

1,4 

19. 

NjOj  -  N03  +  N02(2N02+0) 

4,6 

S.  N20  -*  N2  +0(lD) 

4 

20. 

N03  -  NO  +  02 

4 

6.  NO  —  N  +  O 

5 

21. 

no3  -  no2  +  O 

4 

7.  HNOj  -  OH  +  NOj 

4 

22. 

HjO  -  H  +  OH 

9 

8.  Hj02  -  20H 

4 

23. 

HN04  -  HO  +  NOj  (HOj  +  NOj) 

10,6 

9.  H02  -  OH  +0 

4 

24. 

HOC1  -  OH  +  Cl 

4 

io.  ciono2  -  a  +  no3(cio  +  no2) 

4,6 

25. 

CHjOOH  -  CHjO+OH 

4 

11.  HCI  -  H  +  Cl 

4 

26. 

CHjO  -  HCO  +  H 

11 

12.  CIO  -  Cl  +  O 

7 

27. 

CHjO  -  CO  +  H2 

11 

13.  CIO  -  CI  +  Ol'D) 

Not  used 

28. 

CHjCI  -  CH3  +  Cl 

4 

14.  aN02  -  Cl  +N02 

4 

29. 

HONO  -  OH  +  NO 

4 

15.  OCtO  -  CIO  +0 

7 

Notes  to  Table  A-4. 

1 .  Contributes  to  optical  depth  of  model  atmosphere. 

2.  Schumann-Runge  bands  are  given  special  treatment  based  on  Hudson  and  Mahle  (1972). 

3.  Quantum  yields  of  reactions  (2)  and  (3)  are  given  special  treatment  based  on  temperature-dependent  treatment  of  JPL  (1979). 

4.  Based  on  the  data  of  JPL  (1979).  Where  data  for  several  temperatures  is  given,  we  have  used  the  data  at  230  K. 

5.  Nitric  oxide  photolysis  is  based  on  treatment  of  Frederick  and  Hudson  (1979).  We  have  used  the  photolysis  rates  averaged 
over  the  sunlit  hemisphere  for  day  time  photolysis  rates.  1979a  chemistry  used  the  treatment  of  Cieslik  and  Nicolet  (1973). 
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6.  The  products  used  furXNO^  (X  =  Cl.  OH,  NO2)  changed  between  our  1979a  and  1979b  chemistries.  I  or  1979a  they  were 
based  on  the  path  of  lowest  endoergicity  (except  for  l^O^  which  was  based  on  a  recommendation  of  Johnston),  f  or  1979b 
they  arc  all  based  on  analogy  with  CIONO2  data  of  Chang  et  al.  (1979).  This  treatment  is  highly  uncertain. 

7.  Based  on  data  in  Watson  (1974). 

8.  Product  chemistry  has  been  simplified. 

9.  Treatment  based  on  Thompson  et  al.  (1963). 

10.  Treatment  based  on  Graham  et  al.  (1978). 

1 1.  Treatment  based  on  quantum  yields  of  Moortgat  and  Warneck  (1979)  and  cross  sections  of  McQuigg  and  Calvert  (1969). 


BOUNDARY  CONDITIONS 

The  model  now  allows  for  either  fixed  concentrations  or  a  flux  condition  at  the  surface  as  a  lower 
boundary  condition.  For  most  of  the  calculations  in  this  study,  six  species  were  assumed  to  have  fixed  concen¬ 
trations  (See  Table  A-5),  while  a  surface  flux  was  assigned  to  the  other  species.  Zero  flux  was  assumed  except 
for  those  species  shown  in  Table  A-5.  When  those  species  with  fixed  boundary  conditions  in  Table  A-5  were 
given  flux  boundaries,  a  flux  was  determined  to  give  an  ambient  concentration  the  same  as  those  in  Table  A-5. 

Zero  flux  was  assumed  for  all  species  except  NO  and  N02  at  the  upper  boundary.  NO  and  N02  are 
assumed  to  have  a  very  small  flux  from  the  mesosphere  into  the  stratosphere. 

Water  vapor  concentrations  are  fixed  in  the  troposphere  and  calculated  in  the  stratosphere.  All  runs 
are  made  with  fixed  boundary  conditions  unless  otherwise  noted. 


SOURCES  AND  SINKS 

In  addition  to  sources  and  sinks  from  the  chemistry  and  boundary  conditions  described  above,  there 
are  additional  sinks  due  to  dry  and/or  wet  removal  for  many  species  in  the  model.  A  source  for  nitric  oxide 
from  cosmic  ray  dissociation  of  N2  is  also  included  based  on  the  results  of  Nicolet  ( 1974). 

Wet  removal  processes  are  parameterized  by  a  first-order  loss  rate.  The  wet  removal  of  the  trace 
species  HN03,  H202,  HC1,  CIO,  C10N02,  C1N02,  HN04,  HOC1,  CH20  and  CH,OOH  is  assumed  to  vary 
with  altitude  as  shown  in  Table  A-6.  N02  is  assumed  to  have  a  loss  rate  half  the  above  rate. 

Dry  deposition  rates  at  the  surface  are  also  parameterized  by  a  first-order  loss  rate  as  shown  in 
Table  A-7. 


TABLE  A-5.  Boundary  conditions. 


Fixed  concentrations  (molecules  cm  2) 

-2  -K 

Surface  flux  (molecules  cm  s  / 

NjO 

8.25  x  lO12 

NO 

3.41  X  109 

ch4 

4.03  x  1013 

no2 

6.59  X  109 

»2 

1.42  x  1013 

HN03 

1.67  X  109 

CHjCI 

1.56  x  1010 

HCI 

3.67  x  1010 

h2o 

4.30  x  I017 

CCI4 

1.40  x  I06 

CO 

3.04  x  I012 

CF2CL2 

function  of  time  and  scenario 

CFCIj 

CH3CCI3 
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TABLE  A-6.  Wet  removal  parameterization. 


Altitude,  km 

Loss  rate,  sec  * 

Altitude,  km 

Loss  rate,  sec  * 

0 

3.86  x  I0"6 

6 

1.93  X  10"6 

1 

3.86  X  10"6 

7 

1.93  X  10'7 

2 

3.86  X  10'6 

8 

9.58  X  10'7 

J 

3.86  X  10'6 

9 

4.78  x  10'7 

4 

3.86  X  10'6 

10 

0 

5 

3.86  X  10'6 

TABLE  A-7.  Deposition  rates  in  the  lowest  layer  (z  =  0). 


Species 

Loss  rate,  sec'* 

Species 

Loss  rate,  sec'* 

0(3P) 

2.0  x  10'S 

no3 

2.0  X  10'5 

°3 

1.0  x  10'S 

N2°5 

2.0  X  10'5 

NO 

1.0  x  106 

h2o 

0 

no2 

3.0  X  106 

hno4 

2.0  X  10'S 

n2o 

0 

HOC1 

2.0  X  10'5 

hno3 

2.0  X  105 

HCO 

2.0  X  10'5 

OH 

2.0  X  10'5 

ch2o 

1.0  X  10'5 

ho2 

2.0  X  10  5 

ch3 

2.0  X  10'5 

h2°2 

2.0  x  10'5 

CHjOOH 

2.0  x  10  S 

a 

2.0  x  10'5 

ch3o 

2.0  X  10'5 

ClONOj 

2.0  x  105 

CH3°2 

2.0  X  10  5 

CIO 

2.0  x  10  s 

CO 

5.0  x  10'7 

ch4 

0 

ao3 

2.0  X  10"S 

»2 

1.0  x  10'7 

ocio 

2.0  X  10'S 

CHjCI 

0 

HONO 

2.0  X  10'5 

CINOj 

2.0  x  10  s 

CFC13 

0 

ho 

2.0  x  10  s 

CF2CI2 

0 

ca4 

0 

CH3CCI3 

0 

MULTIPLE  SCATTERING 


In  order  to  accurately  compute  photodissociation  rates,  it  is  important  to  describe  radiative 
processes,  such  as  multiple  scattering,  in  addition  to  attenuation  by  gases  such  as  CL,  0\  and  N02.  The  impor¬ 
tance  of  molecular  scattering  varies  significantly  with  wavelength  over  the  spectral  range  200-750  nm.  At 
shorter  wavelengths  gaseous  absorption  dominates  and  very  little  solar  tlux  reaches  the  lower  atmosphere. 
The  region  290-330  nm  is  a  transition  region  where  molecular  scattering  is  very  important,  especially  at  the 
longer  wavelengths  in  this  interval.  Beyond  330  nm,  surface  reflection  is  very  important  since  the  atmosphere 
is  nearly  transparent  in  this  spectral  region. 

Multiple  scattering  is  included  in  the  model  using  a  simplified  method  that  is  computationally  fast  so 
that  it  can  be  used  for  diurnal  calculations.  The  method  is  similar  to  that  of  Isaksen  ei  at.  (1977)  in  terms  of  the 
numerical  method  but  quite  different  in  terms  of  the  physical  assumptions.  The  atmosphere  is  divided  into  op¬ 
tically  thin  layers  and  each  layer  can  absorb  and  scatter  radiation.  The  layer  is  assumed  to  scatter  radiation 
isotropically  with  half  of  the  scattered  flux  going  upward  and  the  other  half  downward  at  an  average  zenith 
angle  of  ±0.  The  earth's  surface  is  also  assumed  to  scatter  isotropically,  and  a  surface  albedo  of  0.25  is  used  to 
approximate  the  ef  fect  of  clouds  on  the  upward  scattered  radiation.  Using  a  high  surface  albedo  and  no  clouds 
gives  results  that  are  nearly  identical  to  those  from  dividing  the  atmosphere  into  clear  and  cloudy  regimes  and 
averaging  the  results  (the  exception  being  the  region  below  the  cloud  layer,  which  is  not  important  for  the 
model  applications  considered  here).  For  each  atmospheric  layer  there  is  a  contribution  to  the  solar  flux  den¬ 
sity  from  the  direct  flux  and  by  the  diffuse  fluxes  incident  on  the  layer  from  above  and  from  below .  The  flux 
density  due  to  the  various  fluxes  together  can  be  much  greater  (depending  on  the  wavelength  and  altitude) 
than  the  flux  density  computed  considering  only  gaseous  absorption  (Luther  and  Gelinas,  1976;  Luther  and 
Wuebbles,  1976). 


TEMPERATURE  FEEDBACK 

The  temperature  profile  above  1 3  km  is  calculated  using  a  stratospheric  radiative  transfer  model,  and 
the  temperature  profile  is  specified  at  lower  altitudes.  The  model  includes  solar  absorption  and  long-wave  in¬ 
teraction  by  Oy,  H20,  and  C02,  along  with  solar  absorption  by  N02.  The  techniques  adopted  for  treating 
long-wave  radiative  transfer  are  the  same  as  those  described  by  Ramanathan  (1974).  This  formulation  was 
chosen  because  it  is  computationally  efficient,  and  its  accuracy  has  been  demonstrated  (Ramanathan,  1974, 
1976)  by  comparison  with  much  more  complex  models.  The  effects  of  and  justification  for  the  simplify  ing 
assumptions  used  in  the  model  are  discussed  by  Ramanathan  (1976). 

A  band  absorptance  formulation  is  used  to  treat  the  9.6-gm  band  of  Oy  and  the  fundamental  and 
several  hot  and  isotopic  bands  of  C02  in  the  15-um  region.  An  emissivity  formulation  is  used  to  treat  long¬ 
wave  radiative  transfer  by  HiO.  Solar  absorption  by  Oy  is  treated  by  using  the  empirical  formulation  given  by 
Lindzen  and  Will  (1973).  The  band  absorptance  formulation  by  Houghton  (1963)  is  adopted  for  solar  absorp¬ 
tion  by  H20,  and  the  band  absorptance  formulation  by  Ramanathan  and  Cess  (1974)  is  adopted  for  solar  ab¬ 
sorption  by  C02.  The  empirical  formulation  of  Luther  (1976)  is  used  for  solar  absorption  by  N02.  Solar  ab¬ 
sorption  by  Oy  and  NCL  are  treated  independently  because  absorption  by  these  species  is  weak  in  the  region 
where  their  absorption  bands  overlap.  Solar  radiation  scattered  from  the  troposphere  is  included  by  assuming 
an  albedo  of  0.3.  Doppler  broadening  effects  are  included  for  C02  and  Oy  as  described  in  Appendix  B  of 
Ramanathan  (1976),  The  temperature  dependences  of  the  band  absorptance  and  band  intensity  are  included 
in  the  longwave  calculations  of  COi  and  Oj. 

A  single  cloud  layer  is  included  at  6.5  km  with  42%  cloud  cover  as  was  suggested  by  Cess  (1974).  The 
lapse  rate  within  the  troposphere  is  assumed  to  be  -6.5  K/km,  and  the  temperature  at  the  earth’s  surface  is 
specified  to  be  288  K. 
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NUMERICAL  METHOD 


Each  of  the  39  species  in  the  model  has  its  concentration  calculated  at  each  of  44  vertical  levels  ex¬ 
tending  from  the  surface  to  55  km.  The  numerical  technique  used  to  solve  the  set  of  over  2000  differential 
equations  (resulting  from  a  continuity  equation  for  each  species  at  each  grid  level)  is  the  method  described  by 
Chang  et  al.  ( 1974).  The  main  advantage  of  this  method,  which  is  a  variable  order,  multistep,  implicit  method, 
is  its  ability  to  solve  sets  of  mathematically  stiff  differential  equations  for  almost  any  set  of  input  parameters, 
initial  and  boundary  conditions,  in  particular  those  resultirg  from  the  chemical  kinetics  system  described  in 
Table  A-2. 


DIURNAL  AND  DIURNAL-AVERAGED  MODELS 

We  have  developed  a  fully  diurnal-averaged  model  that  is  consistent  with  our  diurnal  model.  The 
diurnal  model  is  used  to  generate  species  profiles  for  comparison  with  measurements  and  for  perturbation 
studies  involving  short  time  integrations  (e.g.,  solar  eclipse  effects).  The  diurnal-averaged  model  is  used  for 
perturbation  and  sensitivity  studies  involving  longer  time  integrations. 

The  procedure  that  is  used  in  developing  the  fully  diurnal-averaged  model  is  also  applicable  to  two- 
dimensional  models.  If  the  continuity  equation  is  averaged  over  a  time  period  (24  hours  in  our  case)  that  is 
very  small  compared  to  the  time  scale  of  the  problem  of  interest,  then  one  obtains  averaged  terms  of  the  form 
kjjCjCj  and  JjCj  where  c(  is  the  concentration  of  species  i  at  time  t  and  altitude  z,  k,j  is  the  two-body  chemical  rate 
coefficient,  and  J,  is  the  photodissociation  rate  coefficient  for  species  i. 

We  define  the  diurnal  weighting  factors  <vit(z)  and  dj(z)  by 


and 


Vicj =  Vijcicj 


Jic.-M«c. 


(A-5) 


(A-6) 


Since  k(j  is  defined  and  is  independent  of  time,  we  have 


aij 


cicj/cicj 


(A-7) 


and 


=  J;  Cj  /  c; 


(A-8) 


The  computation  of  photodissociation  rates  can  be  an  expensive  part  of  stratospheric  model  calcula 
tions,  hence  evaluation  of  J,’s  in  the  diurnal-averaged  model  can  be  expensive.  If  we  define  d\  by 


=  JiCj/(J[u,on  ct)  ,  (A-9) 

then  the  computation  in  the  diurnal-averaged  model  is  greatly  simplified.  The  diurnal  model  is  used  to  deter¬ 
mine  c~cj,  J~c,,  c,.  cj,  and  Jjnoon  so  that  <iM  and  rf,  can  be  obtained  for  every  chemical  and  photochemical  reac¬ 
tion  in  the  model. 
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